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Preface 


This  report  is  the  result  of  an  effort  to  further 
develop  computer  performance  evaluation  (CPE)  tools  for 
personnel  at  the  Air  Force  Avionics  Laboratory, 
Wright-Patterson  AFB,  OH.  A  previously  developed  FORTRAN 
program  to  solve  multi-class  closed  queueing  networks  was 
modified  to  include  an  approximation  algorithm  to  improve 
the  accuracy  of  modeling  job  swapping. 

Not  only  was  this  an  academic  learning  experience,  I 
feel  I  have  gained  some  insight  into  the  problems  of  getting 
support  for  a  worthwhile  project  from  the  management 
structure.  Sometimes  this  can  be  the  most  difficult 
obstacle  in  a  research  project.  Also,  I  have  learned  a  few 
things  about  myself  in  the  process.  I  am  hopefull  that  I 
will  use  this  knowledge  to  take  advantage  of  my  strengths 
and  improve  on  some  of  my  weaknesses. 

First,  I  would  like  to  thank  my  thesis  advisor,  Ltc 
James  Bexfield,  for  his  support  in  this  thesis  effort.  His 
encouragement  and  guidance  were  only  exceeded  by  his 
patience. 
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Also  of  great  moral  support  were  ray  parents.  Their 
many  cards  and  letters  and  long-distance  phone  calls  were 
always  full  of  encouragement. 

I  would  like  to  thank  all  of  ray  classmates  in  both  the 
G0R-81D  and  GCS-81D  for  their  help  through  the  year  with  the 
many  problem  sets,  terra  projects,  and  labs.  A  special 
thanks  goes  to  Ken  Bauer  for  his  artwork  on  the  plots  in 
this  thesis.  More  importantly,  his  light-hearted  humor  and 
warm  smile  were  always  there  to  cheer  me  up.  May  you  never 
lose  your  tremendous  concern  for  others. 

Finally,  I  would  like  to  thank  a  small  software  company 
called  Mark  of  the  Unicorn  for  developing  Scribble,  the  best 
text  formatter  for  a  CP/M  home  computer  I've  ever  seen,  and 
for  their  patient  help  when  I  needed  it. 

Michael  H.  Cox 


-iii- 


f 

\r 


Table  of  Contents 


Preface  ii 

List  of  Figures  viii 

List  of  Tables  ix 

Abstract  x 

Chapter  1  Introduction  1 

1.1  Computer  Performance  Evaluation  (CPE) 

1 

1.1.1  Need  for  CPE  1 

1.1.2  CPE  Techniques  2 

1.1.3  Development  of  Analytical  Models 

in  CPE  5 

1.1. 3.1  Single  Queue,  Infinite 

Population  Model  5 

1.1. 3. 2  Finite  Population  Model  6 

1.1. 3. 3  Central  Server  Model  7 

1.1. 3. 4  Classical  Swapping  Model 

9 

1.1. 3. 5  Chen's  Swapping  Model  10 

1.2  Problem  Statement  11 

1.2.1  Motivation  for  Research 

Objectives  11 

1.2.2  Background  12 

1.2.3  Research  Objectives  13 

1.2. 3.1  Primary  Objective  13 

1.2. 3. 2  Specific  Objectives  13 

1.2.4  Scope  of  the  Thesis  14 

Chapter  2  Operating  Systems  and  CPE  1 

2.1  Operating  System  Objectives  2 

2.2  Resource  Manager  Functions  2 

2.3  The  DECsystera-10  Operating  System: 

TOPS- 10  Monitor  5 

2.3.1  Monitor  Programs  5 

2.3.2  Job  States  7 


i 


l 


-iv- 


i 


Table  o£  Contents  (cont.) 


2.3.3  Monitor  Queueing  Structure  11 

2. 3. 3.1  In-core  Versus  Out-core 

Chains  11 

2. 3. 3. 2  Processor  Queues  11 

2. 3. 3. 3  Long-term  Wait  Queues  12 

2. 3. 3. 4  Processor  Queue  Time  Slices 

14 

2.4  Chapter  Summary  17 

Chapter  3  Analytical  Modeling  1 

3.1  Early  Queueing  Network  Models  1 

3.2  Multi-class  Queueing  Network  Model 

3 

3.2.1  Variables  in  the  Model  4 

3.2.2  Open  Network  Outside  Arrival 

Processes  5 

3.2.3  Service  Time  Distributions  6 

3.2.4  State-dependent  Service  Rates 

7 

3.2.5  Service  Center  Types  7 

3.2.6  States  of  the  Model  9 

3. 2. 6.1  Type  1  Service  Center  9 

3. 2. 6. 2  Type  2  and  Type  3  Service 

Centers  10 

3. 2. 6. 3  Type  4  Service  Centers  10 

3.2.7  Calculating  Equilibrium  State 

Probabilities  11 

3. 2. 7.1  Balance  Equations  11 

3. 2. 7. 2  Product  Form  Solution  12 

3.2.8  Example  Problem  16 

3. 2. 8.1  The  Balance  Equations  17 

3. 2. 6. 2  The  Product  Form  Solution 

24 

3.3  Chen's  Swapping  Model  26 

3.3.1  Variables  in  the  Model  27 


Table  of  Contents  (cont.) 


3. 3. 1.1  Derivation  of  P14  and  P^2 


30 

3.3.2  Program  Swapping  Behavior  32 

3. 3. 2.1  Derivation  of  P43  32 

3. 3. 2. 2  Derivation  of  P13  34 

3. 3. 2. 3  Approximation  Algorithm  36 

3.4  Chapter  Summary  37 

Chapter  4  Computer  Implementation  1 

4.1  The  McKenzie  Program  2 

4.1.1  Program  Capabilities  2 

4.1.2  Program  Structure  3 

4.1.3  Program  Inputs  6 

4.2  The  Chen  Modification  6 

4.2.1  Program  Structure  6 

4.2.2  Program  Inputs  8 

4.3  Programming  Notes  for  Future 

Modification  11 

4.4  Chapter  Summary  13 

Chapter  5  Analytical  Modeling  Results  1 

5.1  Modeling  Interactive/Batch  Workloads 

1 

5.2  DECsystem-10  System  Configuration 

Parameters  2 

5.2.1  Hardware  Parameters  2 

5.2.2  Workload  Parameters  4 

5.2.3  DECsystem-10  Swapping  Model  6 

5.3  Model  Comparisons  8 


vi- 


Table  of  Contents  (cont.) 


5.3.1  Definition  of  Performance 

Measures  8 

5.3.2  Tabular  Results  9 

5.3.3  Chen's  Versus  the  Classical 

Swapping  Model  12 

5. 3. 3.1  Probability  Structure  12 

5. 3. 3. 2  Performance  Predictions  15 

5.3.4  Chen's  Swapping  Model: 

Multi-class  Versus 

Single-class  19 

5.4  Chapter  Summary  24 

Chapter  6  Conclusions  and  Recommendations 

1 

6.1  Conclusions  1 


6.2  Recommendations  for  Future  Research 

2 


Bibliography 


Vita 


Appendix 


A-J 


-vii 


List  of  Figures 


Figure 

Paae 

1-1 

Spectrum  of  Computer  Modeling  Techniques 

1-4 

1-2 

Single-resource  Queueing  Model 

1-6 

1-3 

Finite  Population  Model 

1-7 

1-4 

Central  Server  Model 

1-8 

1-5 

Classical  Swapping  Model 

1-10 

2-1 

Monitor  Cycle  Programs 

2-8 

2-2 

Job  State  Transitions 

2-10 

3-1 

Method  of  Stages 

3-6 

3-2 

State  Transition  Diagram 

3-18 

3-3 

Global  Balance  Equations 

3-19 

3-4 

Local  Balance  Equations 

3-21 

3-5 

Chen's  Swapping  Model 

3-29 

4-1 

McKenzie  Program  Flow  Chart 

4-5 

4-2 

New  Program  Flow  Chart 

4-9 

5-1 

DECsystem-10  Job-Swapping  Model 

5-7 

5-2 

Probability  Structure  (Two  Classes) 

5-13 

5-3 

Two  Classes 

5-16 

5-4 

Chen's  Model 

5-21 

-viii- 


List  of  Tables 


Table  Page 

3-1  Transition  Matrix  and  Service  Rates  3-16 

3-2  Example  Calculations  3-25 

3-3  Probability  Transition  Matrix  for 

Chen's  Swapping  Model  3-31 

5-1  System  Configuration  Parameters  5-3 

5-2  Classical  Swapping  Model  Results  5-9 

5-3  Chen's  Swapping  Model  Performance 

Results  5-10 

5-4  Important  Probability  Transition 

Matrix  Values  5-11 


I 


Abstract 


An  improved  model  of  the  DECsy stem-10  job-swapping 
behavior  was  developed.  This  model  combines  a  previously 
developed  closed  queueing  network  model  with  a  job-swapping 
model  developed  by  Chen/"(Ref  Chen's  swapping  model 
provides  an  approximate  solution  to  a  network  queueing  model 
with  a  state-dependent  probability  transition. 


This  combined  model  is  then  tested  on  a  hypothetical, 
though  realistic  workload  containing  both  interactive  and 
batch  jobs.  The  two  classes  of  jobs  are  treated  first  as 
separate  classes,  as  one  class  having  the  weighted  average 
job  characteristics  of  both  classes,  and  as  one  class  having 
just  the  interactive  job  characteristics.  The  results  of 
these  experiments  and  a  comparison  between  Chen's  swapping 
model  and  the  classical  are  presented. 

The  results  of  the  experiment  indicate  that  it  is 
important  to  model  multiple  classes  for  systems  which  have  a 
significant  amount  of  batch  activity.^Also,  Chen's  swapping 
model  provides  a  more  realistic  mo&^l  of  job-swapping 
behavior  for  the  DECsystem-10.  Therefore,  combining  the 


t 

i 
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multi-class  model  with  Chen's  swapping  model  improves  the 
modeling  accuracy  for  the  DECsystem-10 .  Recommendations  for 
extensions  to  this  multi-class  Chen  model  are  also 
discussed. 
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Chapter  1 


Introduction 


1.1  Computer. Performance  Evaluation  (CPE) 

1.1.1  Need  for  CPE 

In  recent  years,  science  and  engineering  have  become 
more  concerned  with  the  economic  aspects  in  their  fields  due 
to  increased  costs  of  high  technology  and  tighter  fiscal 
budgets.  Great  attention  has  been  given  to  the  development 
and  refinement  of  techniques  which  help  predict  behavior  of 
systems  and  thus  yield  insights  into  what  cost-performance 
tradeoffs  can  be  made  (Ref  23:1). 

Any  system  which  is  in  the  process  of  being  designed, 
procured,  or  modified  must  satisfy  certain  predetermined 
performance  specifications.  This  is  especially  true  in 
computer  engineering  where  system  specifications,  reference 
manuals,  and  user  guides  abound.  Designers  and  engineers 
use  design  and  performance  evaluation  prediction  techniques 
to  obtain  systems  which  meet  these  specifications. 
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Prospective  users  of  a  system  use  these  techniques  to 
determine  which  combination  of  subsystems  of  components 
comes  closest  to  matching  their  requirements,  given  certain 
cost  constraints.  Current  users  of  a  system  use  these  same 
evaluation  techniques  and  tools  to  help  make  decisions  about 
existing  systems  concerning  system  upgrades  and  additions. 
Therefore,  performance  evaluation  is  needed  at  all  stages  of 
the  life  cycle  of  a  computer  product.  (Refs  9:1,23:1-2). 

Computer  engineering  is  a  fairly  recent  development  and 
therefore  performance  evaluation  is  less  developed  than 
older  branches  of  engineering,  but  definitely  is  not  less 
important  .  The  products  of  computer  engineering  (central 
processing  units  (CPU's),  memory,  Printers,  card  readers, 
tape  drives,  disk  drives,  etc.)  are  primarily  designed  to 
perform  certain  functions  related  to  the  processing  of 
information.  How  well  these  systems  execute  their  tasks  is 
a  matter  of  tremendous  technical,  enonomic,  and  military 
importance  to  the  Air  Force  (Ref  9:3) 

1.1.2  CPE  Techniques 

Techniques  for  evaluating  a  computer  system  can  be 
divided  into  two  categories:  measurement  techniques  and 
modeling  techniques.  Measurement  techniques  involve  using 
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performance  monitors  to  detect  and  analyze  system  events. 
These  monitors  can  be  hardware  oriented  or  software 
oriented,  and  their  use  with  benchmark  programs  is  called 
benchmarking.  Modeling  techniques,  on  the  other  hand,  do 
not  directly  involve  the  system  in  question,  but  deal  with  a 
conceptual  representation  of  the  system  (Refs  1,  9,  29). 

These  modeling  techniques  can  be  divided  into  four 
categories  (Ref  1)  : 

1.  Rules  of  thumb,  e.g.  CPU  utilization  should  not 
exceed  35  percent  for  on-line  applications  or  40 
percent  for  batch  applications. 

2.  Linear  projection,  e.g.  "Computer  useage  doubled  in 
the  last  year,  so  it  will  double  again  next  year." 

3.  Analytical  queueing  models. 

4.  Simulation  models. 

Another  CPE  technique  discussed  by  Svobodova  (Ref  29:48-49) 
and  more  extensively  by  Sanabria  (Ref  25)  is  empirical 
modeling.  This  technique  combines  the  linear  projection 
technique  with  benchmarking.  Measurement  data  are  collected 
from  the  real  computer  system  and  regression  or  some  other 
curve-fitting  technique  is  used  to  develop  performance 
prediction  curves. 
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Figure  1-1  Spectrum  of  Computer  Modeling 
Techniques  (Ref  1:13) 

These  evaluation  techniques  involve  increases  in 
accuracy  as  one  moves  from  the  very  simplistic  rules  of 
thumb  to  the  very  complex  benchmarking  techniques  (see 
Figure  1-1) .  Unfortunately,  the  same  increase  in  accuracy 
also  involves  an  increase  in  cost.  The  computer  engineer 
involved  in  CPE  must  then  make  tradeoffs  as  to  the  amount  of 
accuracy  he  can  afford.  These  tradeoffs  are  extremely 
prevalent  when  comparing  the  two  most  widely  used  types: 
simulation  and  analytical  queueing  models. 

There  are  three  reasons  for  the  increasing  popularity 
of  analytical  models  over  simulation  models  for  computer 
system  modeling  (Refs  1,  4,  11,  21,  23,  26,  and  29): 

-  Analytical  models  capture  the  most  important  features 
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of  actual  systems,  i.e.  jobs  moving  from  one  queue  to 
another  waiting  for  service  from  independent  devices 
within  the  system. 

-  The  assumptions  of  the  analysis  are  realistic.  General 
device  service  time  distributions,  load-dependent 
devices  service  times,  and  multiple  classes  of  jobs  can 
be  modeled. 

-  The  algorithms  that  solve  the  equations  of  the  model 
are  available  as  highly  efficient  queueing  network 
evaluation  packages.  Because  of  their  efficiency  and 
simplicity,  these  models  are  cheaper  than  simulations 
to  develop  and  run. 


1.1.3  Development  of  Analytical  Models  in  CPE 

The  development  of  more  flexible  analytical  modeling 
techniques  for  computer  system  modeling  has  paralleled  the 
evolution  of  computer  systems  from  single-programmed,  batch 
systems  to  multi-programmed,  combination  interactive/batch 
systems. 

1.1. 3.1  Single  Queue,  Infinite  Population  Model 

The  earliest  techniques  modeled  the  entire  computer 
system  as  a  single  resource  with  a  single  queue  and  an 
infinite  population  of  jobs  (see  Figure  1-2).  This  model  was 
used  for  the  early  batch  computers,  since  the  CPU  was  the 
dominate  resource  and  only  one  job  could  be  in  the  system  at 
a  time. 
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Figure  1-2  Single-resource  Queueing  Model  (Ref  21:947). 


1.1. 3. 2  Finite  Population  Model 

With  the  advent  of  interactive  systems,  modeling  the 
system  as  a  single-resource  queue  with  an  infinite 
population  of  jobs  became  less  realistic.  Because  there  are 
a  finite  number  of  jobs  circulating  within  an  interactive 
system,  the  rate  of  arrival  of  new  requests  for  service  will 
tend  to  decrease  as  the  queue  length  grows.  The  finite 
population  model,  also  known  as  the  machine  interference 
model,  was  used  to  model  this  phenomena  (see  Figure  1-3).  In 
the  early  interactive  systems,  only  one  complete  job  could 
be  in  memory  at  a  time.  The  execution  of  the  jobs,  I/O 
activity,  and  the  swapping  activity  were  not  overlapped. 
Therefore,  program  execution,  I/O  activity,  and  the  swapping 
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time  were  summed  and  used  as  the  CPU  service  time  (Ref 
21:946-947)  . 


TERMINALS 


< 


PREEMPTED  JOBS 


Figure  1-3  Finite  Populaton  Model  (Ref  21:947). 
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1.1. 3. 3  Central  Server  Model 

With  the  development  of  true  multi-programming,  i.e. 
more  than  one  entire  job  within  memory,  the  finite 
population  models  became  less  realistic.  The  most 
fundamental  characteristic  of  jobs  in  a  computer  system  is 
that  they  alternate  between  CPU  execution  and  being  blocked 
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from  further  CPU  execution  waiting  for  access  to  secondary 
storage.  With  more  than  one  job  in  memory,  the  I/O  activity 
of  one  job  could  be  overlapped  with  the  execution  of  another 
job.  A  model  which  considers  this  type  of  behavior  is  shown 
in  Figure  1-4.  Note  that  this  model  still  does  not  consider 
swapping  times  separately. 


TERMINALS 
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SECONDARY  STORAGE 
DEVICES 


PREEMPTED  JOBS 
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Figure  1-4  Central  Server  Model  (Ref  3) . 


1-8 


Introduction 


1.1. 3. 4  Classical  Swapping  Model 

As  more  and  more  computing  systems  used  swapping  as 
their  memory  management  strategy,  it  became  necessary  to 
more  accurately  model  the  job  swapping  activities.  The 
classical  method  to  model  job  swapping  explicitly  was  to  add 
a  swapping  device  node  after  the  terminal  interaction  node 
(see  Figure  1-5)  .  Thus,  each  job  is  swapped  out  after  each 
interaction  with  the  user  at  the  terminal  (Ref  20) .  To 
improve  the  overall  accuracy,  two  other  improvements  to  the 
interactive  model  were  developed  (Ref  2:249) 

1.  Multiple  classes  of  jobs  with  each  class  having  its 
unique  transition  probabilities  and  node  service 
times. 

2.  Non-exponential  service  times  using  the  method  of 
stages. 

A  more  detailed  account  of  these  improvements  and  how  they 
were  integrated  into  one  model  is  given  in  Chapter  3. 

1.1. 3. 5  Chen's  Swapping  Model 

Chen  (Ref  5)  proposed  a  new  approach  to  modeling 
job-swapping  behavior:  that  it  be  modeled  by  state-dependent 
transition  probabilities  in  a  closed  queueing  network  with  a 
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single  class  of  jobs.  Since  an  exact  solution  to  this 
problem  does  not  currently  exist,  Chen  developed  an 
approximation  algorithm  that  iteratively  solves  a  closed 
queueing  network.  Chen's  model  is  described  in  more  detail 
in  Chapter  3. 


TERMINALS 


Pigure  1-5  Classical  Swapping  Model  (Ref  5) . 
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1.2  Problem  Statement 

1.2.1  Motivation  for  Research  Objectives 

The  Avionics  Laboratory  at  Wright-Patterson  AFB  is 
responsible  for  research,  development,  and  validation  of 
avionics  systems  for  present  and  future  aircraft  in  the  Air 
Force  inventory.  This  includes  the  development  and 
validation  of  both  the  hardware  (instrument  displays, 
control  consoles,  etc.)  and  software  (the  computer  programs 
that  control  the  hardware)  for  navigational  aids,  weapons 
delivery  systems,  and  electronic  warfare  systems.  The 
laboratory  is  one  of  the  most  advanced  facilities  for  the 
development  and  validation  of  avionics  systems. 

One  of  the  primary  test  beds  for  avionics  system 
development  is  the  Digital  Electronics  Corporation's 
DECsystem-10  mainframe  computer.  The  DECsystem-10  is  a 
multi-programming,  multi-processing,  time-sharing  computer 
with  real-time  processing  capabilities.  It  provides 
high-level  language  support  (FORTRAN,  PASCAL,  COBOL,  BLISS, 
JOVIAL,  and  in  the  future  ADA);  assembly  language  support; 
document  and  manual  preparation  programs  such  as  editors, 
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text  formatters,  and  spelling  checkers;  and  graphics  and 
plotting  programs.  Several  mini-computers  (PDP-ll's)  are 
directly  interfaced  to  the  DECsystem-10.  This  allows 
real-time  communication  between  the  DECsystem-10  and  the 
PDP-11  which  are  used  to  control  testing  and  simulation  of 
current  and  experimental  avionics  systems. 

Since  the  DECsystem-10  is  such  a  vital  Air  Force 
resource  in  the  research  and  development  of  state-of-the-art 
avionic  systems,  it  is  important  that  this  resource  is  used 
as  efficiently  and  effectively  as  possible.  To  insure  the 
proper  management  of  this  resource,  the  system  managers  must 
have  the  proper  tools  to  predict  the  impact  of  future 
workload  changes,  potential  operating  system  changes,  and 
reconfiguration  of  hardware  resources.  With  these  tools, 
the  managers  of  the  DECsystem-10  will  gain  a  better 
understanding  of  the  factors  impacting  computer  performance 
and  will  be  better  able  to  make  decisions  to  properly  manage 
the  DECsystem-10. 

1.2.2  Background 

In  1977  McKenzie  developed  a  computer  program  to  solve 
multi-class  closed  queueing  networks  with  service  centers 
having  different  service  disciplines.  This  program  was  used 
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to  solve  simple  finite  population  and  central  server  type 
queueing  models  of  the  DECsystem-10 ,  the  results  of  which 
helped  the  Avionics  Lab  to  justify  the  purchase  of 
additional  memory.  Both  the  results  from  McKenzie's  limited 
validation  and  Saxton's  entire  thesis  (Ref  24),  indicated 
that  more  detailed  models  of  the  DECsystem-10  job-swapping 
behavior  were  needed  (Refs  19:141  and  24:107). 

1.2.3  Research  Objectives 

1.2. 3.1  Primary  Objective 

The  primary  objective  of  this  thesis  effort  was  to 
improve  the  models  of  the  DECsystem-10  by  better  modeling 
the  job-swapping  behavior.  This  will  enable  the  system 
administrators  to  make  more  informed  decisions  concerning 
the  present  use  of  the  computer  system  resources,  as  well  as 
more  accurate  and  responsive  planning  for  future  computer 
system  acquisitions. 

1.2. 3. 2  Specific  Objectives 

In  order  to  meet  the  above  primary  objective,  the 

following  tasks  were  accomplished  in  this  thesis  effort: 

1.  The  McKenzie  program  to  compute  performance  measures 
for  closed,  multi-class  queueing  networks  was  brought 
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up  on  the  CDC  computer.  The  program  had  to  be  typed 
in  from  a  listing  in  the  thesis  and  subsequently 

debugged  in  order  to  make  it  runnable. 

2.  The  program  was  modified  to  include  an  approximation 

algorithm  developed  by  Chen  (Ref  5)  to  better  model 

the  DECsystem-10  swapping  behavior.  By  combining  the 
multi-class  modeling  capability  of  the  McKenzie 
program  with  Chen's  algorithm,  we  now  have  a  new,  more 
powerful  tool  to  accurately  model  the  performance  of 
an  interactive  computer  system  with  job-swapping. 

3.  The  modified  program  was  run  using  a  workload 
consisting  of  interactive  and  batch  jobs  and  the 
results  used  to  answer  the  following  questions: 

-  Does  Chen's  swapping  model  improve  the  accuracy 
of  the  performance  predictions  when  compared  to 
the  classical  swapping  model? 

-  Does  combining  Chen's  swapping  model  and  multiple 
classes  improve  the  accuracy  of  the  performance 
predictions?  Specifically,  is  it  important  to 
model  an  interactive  computer  system  containing 
both  interactive  and  batch  jobs  with  a 
multi-class  model  versus  using  a  single  class 
model  where 

*  The  single  class  has  the  weighted  average  of 
the  characteristics  of  the  batch  and 
interactive  jobs  in  the  multi-class  model; 
or 

*  The  single  class  has  the  characteristics  of 
the  largest  class  of  jobs  in  the  system? 


1.2.4  Scope  of  the  Thesis 


The  the  next  chapter  will  discuss  the  objectives  and 
basic  functions  of  computer  operating  systems.  This  will  be 
followed  by  a  description  of  the  DECsystem-10 ' s  operating 
system  in  enough  detail  to  determine  how  it  affects  job 
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swapping.  Chapter  3  provides  a  detailed  discussion  of  the 
state-of-the-art  analytical  queueing  model  with  a  product 
form  solution  developed  by  Basket,  Chandy,  Muntz,  and 
Palacios.  Chapter  3  concludes  with  a  detailed  description  of 
Chen' 3  swapping  model  and  its  underlying  assumptions. 
Chapter  4  contains  a  description  of  McKenzie's  computer 
implementation  and  how  it  was  modified  to  include  Chen's 
algorithm.  Chapter  5  presents  the  results  of  modeling  a 
hypothetical  workload  using  the  modified  model  to  answer  the 
above  questions.  Finally,  Chapter  6  will  present  a  summary 
of  the  conclusions  as  well  as  recommendations  for  future 
work  in  this  area. 
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The  primary  method  of  improving  the  performance  of  an 
existing  computer  system  without  buying  new  and  improved 
hardware  is  to  modify  the  behavior  of  the  program  that 
controls  the  computer  system.  This  program  is  called  the 
operating  system.  Operating  systems  usually  have  built-in 
software  parameters  which  may  be  set  to  allow  the  system 
administrator  to  implement  various  scheduling  policies  and 
priority  schemes.  Poorly  chosen  settings  of  these  operating 
system  parameters  can  cause  inefficient  use  of  system 
resources.  Therefore,  it  is  important  that  the  system 
administrator  has  the  tools  to  properly  determine  these 
operating  system  parameters.  Before  discussing  this  problem 
in  detail  and  how  it  relates  to  the  specific  case  dealt  with 
in  this  paper,  we  will  first  review  some  of  the  basic 
objectives  and  functions  of  operating  systems. 
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2.1  Operating  System  Objectives 


The  operating  system  is  a  collection  of  system  programs 
(algorithms)  designed  to  meet  the  following  objectives  (Ref 
6:2)  : 


Provide  the  programmers  with  an  efficient  environment 
for  program  development,  debugging,  and  execution. 

Provide  a  range  of  problem-solving  facilities 

(application  programs). 

Provide  all  this  at  the  lowest  cost  by  sharing 
resources  and  information. 


2.2  Resource  Manager  Functions 

In  order  to  meet  the  above  objectives,  the  operating 
system  must  efficiently  manage  the  systems's  resources. 
Therefore,  its  primary  role  is  that  of  resource  manager.  It 
must  accomplish  the  following  (Ref  18:8): 

1.  Keep  track  of  the  resources. 

2.  Enforce  policy  that  determines  who  gets  what,  when, 
and  how  much. 

3.  Allocate  the  resources. 

4.  Reclaim  the  resources. 
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The  operating  system's  primary  role  can  be  broken  down 
by  resource  type:  memory,  processors,  devices  (disk  drives, 
tape  drives,  printers,  etc.),  and  information  (programs  and 
data  such  as  editors,  compilers,  file  directories  and  other 
software  resources) .  Below  are  listed  the  resource 
management  functions  by  type  and  typical  names  given  to  some 
of  the  routines  that  perform  these  functions  (Ref  18:9-10) : 


Memory  Management  Functions 

1.  Keep  track  of  the  resource  (memory),  what  parts 
are  in  use  and  by  whom?  What  parts  are  not  in 
use  (called  free) ? 

2.  If  multiprogramming,  i.e.  more  than  one  program 
in  execution  at  one  time,  decide  which  job  gets 
memory,  when  it  gets  it,  and  how  much. 

3.  Allocate  the  resource  (memory)  when  the  jobs 
request  it  and  the  policy  of  2  above  allows  it. 

4.  Reclaim  the  resource  (memory)  when  the  job  no 
longer  needs  it  or  has  been  terminated. 

CPU  Processor  Management  Functions 

1.  Keep  track  of  the  resource  (processors  and  the 
status  of  jobs) .  The  system  program  that  does 
this  has  been  called  the  traffic  controller. 

2.  Decide  who  will  have  a  chance  to  use  the 
processor;  the  iob  scheduler  chooses  from  all  the 
jobs  submitted  to  the  system  and  decides  which 
one  will  be  allowed  into  the  system,  i.e.  have 
resources  assigned  to  it.  If  multiprogramming, 
decide  which  job  gets  the  processor,  when,  and 
how  much. 

3.  Allocate  the  resource  (processor)  to  a  job  by 
setting  up  necessary  hardware  registers;  this 
system  program  is  often  called  the  dispatcher. 
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4.  Reclaim  the  resource  (processor)  when  the  job 
relinquishes  processor  usage,  exits  from  the 
system  after  job  completion,  or  exceeds  allowed 
amount  of  usage  and  aborts. 

Pevice  Management  Functions 

1.  Keep  track  of  the  resource  (disk  and  tape  drives, 
channels,  control  units) ;  this  is  typically 
called  the  I/£  fcmffis  controller. 

2.  Decide  what  is  an  efficient  way  to  allocate  the 
resource  (device) .  If  it  is  to  be  shared,  then 
decide  who  gets  it,  and  how  much  he  is  to  get; 
this  is  called  I/O  scheduling. 

3.  Allocate  the  resource  (device)  and  initiate  the 
I/O  operation. 

4.  Reclaim  the  resource  (device).  In  most  cases  when 
the  I/O  terminates,  the  device  is  released 
automatically,  but  the  operating  system  must  be 
informed  of  the  new  status  of  the  device. 

Information  Management  Functions 

1.  Keep  track  of  the  resource  (information),  its 
location,  use,  status,  etc.  These  collective 
facilities  are  often  called  the  file  system. 

2.  Decide  who  gets  use  of  the  resources,  enforce 
protection  requirements,  and  provide  acessing 
routines. 

3.  Allocate  the  resource  (information),  e.g.,  open  a 
file. 

4.  Deallocate  the  resource,  e.g.,  close  a  file. 
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2.2  Ihe-DECSYStemrlfl  Operating  Systems  TQPS-10  Monitor 

The  TOPS-IO  Monitor  performs  the  accounting, 
scheduling,  resource  allocation,  and  service  routines 
necessary  to  operate  in  a  multiprogramming,  time-sharing 
environment.  It  both  controls  user  jobs  and  provides 
services  to  them.  The  monitor  gives  the  appearance  of  a 
single-user  machine  to  all  the  users  on  the  DECsystem-10  by 
rapidly  switching  control  to  each  user.  It  manages  all  I/O 
operations,  according  to  requests  from  user  programs  and 
from  device  interrupts.  It  attempts  to  allocate  all  system 
resources  in  such  a  way  as  to  give  the  best  overall  system 
performance  (Ref  8:INTRO-10). 

2.3.1  Monitor  Programs 

The  monitor  consists  of  many  separate  and  more  or  less 
independent  programs  which  are  called  according  to  events 
which  occur  within  the  system.  The  system  programs  are 
divided  into  two  types  depending  on  whether  the  program 
executes  synchronously  or  asynchronously.  Synchronous 


1.  The  clock  interrupt  is  an  AC  line  generated  interrupt  and 
therefore  occurs  every  1/60  of  a  second.  This  time  interval 
has  been  given  the  name  "jiffy"  and  the  term  occurs 
frequently  in  DECsystem-10  literature. 
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system  programs  execute  once  between  system  clock 
interrupts1  and  complete  their  function  before  the  next 
clock  interrupt.  After  each  clock  interrupt,  the  monitor 
cycle  is  executed  followed  by  the  running  of  a  user  job 
selected  by  the  Job  Scheduler  (Ref  8) .  The  synchronous 
system  programs  are  the 

CLOCK1  The  system  program  which  acts  as  the  "main" 

calling  program  for  all  other  system 

programs.  This  program  contains  the  code 
used  by  the  Job  Scheduler  to  make  the  actual 
resource  assignments. 

Job  Scheduler  The  system  program  which  selects  the  next  job 
to  run  in  the  remainder  of  the  jiffy  and 
controls  the  allocation  of  system  resources. 
This  program  is  also  known  as  SCHED. 

Command  Processor 

The  system  program  which  reads  in  a  typed 
command  and  interprets  it. 

Context  Switching  Routine 

The  system  program  which  transfers  control  to 
the  user  job  selected  by  the  Job  Scheduler. 


Asynchronous  system  programs  execute  on  an  as  needed 
basis  by  jobs  within  the  system  and  may  have  program  cycles 
lasting  longer  than  one  jiffy.  The  asynchronous  system 
programs  are  described  below 

Job  Swapper  The  system  program  called  by  the  Job 
Scheduler  to  keep  as  many  runnable  jobs  in 
memory  as  possible.  Even  though  it  is  called 
each  time  during  the  execution  of  the  Job 
Scheduler  (a  synchronous  program) ,  the  Job 
Swapper  is  not  synchronous  since  the  swapping 
in  and  out  of  a  job  usually  takes  several 
jiffies  to  complete.  The  Job  Swapper 
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"remembers"  what  it  was  doing  when  called  by 
the  Job  Scheduler.  Each  time  the  Job  Swapper 
is  called,  it  continues  processing  where  it 
last  terminated  during  the  last  jiffy, 
o 

UUO  Processor  Routine 

The  system  program  which  executes  whenever  a 
job  requests  I/O  service  such  as  disk  reads 
or  writes. 

Core  Management  Routine 

The  system  program  which  dynamically 
allocates  memory  as  needed  by  other  system 
programs  or  user  jobs. 


The  system  programs  described  above  and  how  they 
interact  with  one  another  during  a  clock  interrupt  is  shown 
in  Figure  2-1.  All  dotted  arrows  represent  asynchronous 
events  which  may  or  may  not  occur  within  a  monitor  cycle, 
depending  on  the  state  of  the  system  (Ref  8) . 


2.3.2  Job  States 

The  next  state  transition  that  a  job  can  enter  depends 
on  the  current  state  of  the  job  (Ref  6:31-35).  The  job  state 
is  represented  by  many  status  variables  maintained  within 
the  monitor  system  program  and  job  status  tables.  A  job  can 
be  in  one  of  four  states  (Ref  8:5-1): 


2.  Unimplemented  User  Operations  (UUO)  are  machine  language 
operation  codes  which  directly  call  the  monitor  to  execute 
supervisor  level  tasks  such  as  I/O  requests  to  system 
peripherals. 
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1.  The  job  is  in  a  processor  queue,  has  memory,  and  is 
running  on  a  CPU. 

2.  The  job  is  in  a  processor  queue,  ready  to  execute  on  a 
CPU.  As  defined  by  DECsystem-10  literature,  a  job  in  a 
processor  queue  is  defined  to  be  ready  to  execute  if 
the  job  is  not  waiting  for  a  shareable  resource-3, even 
though  the  job  may  have  no  memory  and  must  first  be 
swapped  in. 

3.  The  job  is  in  a  short-term  wait  state.  Jobs  in  this 
state  were  previously  in  the  running  state,  but 
requested  the  use  of  a  shareable  resource.  Jobs  in 
the  short-term  wait  state  can  not  be  selected  to 
execute  until  they  are  assigned  the  shareable  resource 
and  are  once  again  runnable. 

4.  The  job  is  in  a  long-term  wait  state.  Jobs  in  this 
state  were  previously  in  the  running  state,  but  then 
requested  the  use  of  a  non-shareable  resource. 
Non-shareable  resources  are  distinguished  from 
shareable  resources  by  the  length  of  time  the  job  may 
have  to  wait  to  be  assigned  the  resource. 
Non-shareable  resources  include  such  things  as  line 
printers  or  card  readers,  as  well  as  less  obvious 
non-shareable  resource  such  as  a  response  from  a 
terminal. 


The  possible  states  of  a  DECsystem-10  are  summarized  in 
Figure  2-2.  The  states  labeled  HPQ  1-15,  PQ1,  and  PQ2  are 
the  processor  queues  and  are  discussed  in  the  following 
section. 


3.  A  shareable  resource  is  some  part  of  the  system,  either 
hardware  or  software,  which  can  be  used  by  only  one  job  at  a 
time,  but  is  shared  among  different  jobs  over  relatively 
short  periods  of  time.  Examples  of  shareable  resources  are 
I/O  channels  and  disk  drives  used  to  satisfy  disk  I/O 
requests 
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2.3.3  Monitor  Queueing  Structure 

Scheduling  in  the  DECsystem-10  is  based  on  the  use  of 
queues  and  wait  state  codes.  The  jobs  in  the  system  are 
maintained  in  a  master  set  of  queues,  which  are  logically 
divided  into  long-term  wait  queues  and  processor  queues. 

2. 3. 3.1  In-core  Versus  Out-core  Chains 

Each  queue  in  the  long-term  wait  queues  and  processor 
queues  are  separated  into  two  chains:  one  contains  all  the 
jobs  in  the  queue  that  are  in-core  and  the  other  contains 
all  the  jobs  that  are  out-core  <have  no  memory).  This 
breakdown  enables  efficient  scanning  of  the  queues  by  the 
Job  Scheduler  and  the  Job  Swapper.  The  scheduler  only  has  to 
scan  jobs  in  the  in-core  chains  to  find  a  runnable  job  with 
memory.  The  Job  Swapper  only  has  to  scan  jobs  in  the 
in-core  chains  to  find  a  job  to  swap  out  to  make  room  for  a 
higher  priority  job  or  scan  the  out-core  chains  to  find  a 
job  that  needs  to  be  swapped  in  (Ref  8:SCH-8). 

2. 3. 3. 2  Processor  Queues 

Jobs  in  the  processor  queues  can  either  be  running, 
ready  to  run  (with  or  without  memory),  or  in  a  short-term 
wait  for  a  shareable  resource.  The  processor  queues  are  the 


2-11 


Operating  Systems  and  CPE 


high-priority  queues  (HPQ's),  PQ1,  and  PQ2.  Each  are 
described  below  (Ref  28:2-1): 

HPQ's  (Up  to  15  levels,  called  HPQ1  through  HPQ15) 

contain  jobs  that  require  real-time  response, 
such  as  the  line-printer  spooler  and  the 
card-reader  spooler  system  programs.  These 
queues  are  scanned  first  by  the  Job  Scheduler 
to  find  an  executable  job. 

PQ1  Contains  jobs  that  require  fast  response, 

such  as  time-sharing  jobs.  This  queue  is  the 
next  most  often  scanned  processor  queue. 

PQ2  Contains  jobs  that  require  long-term 

computing,  such  as  those  that  compile 
programs.  Jobs  in  this  queue  are  scanned 
least  often  by  the  Job  Scheduler,  since  fast 
response  is  not  necessary. 


In  addition  to  providing  a  priority  structure  for  the  Job 
Scheduler  to  scan  for  an  executable  job,  it  establishes  a 
priority  structure  for  the  Job  Swapper.  The  priority 
structure,  though,  is  exactly  reversed.  The  Job  Swapper 
first  scans  PQ2,  PQ1,  and  then  the  HPQ's  to  find  a  job  it 
can  swap  out.  The  Job  Swapper  can  not  swap-out  jobs 
performing  disk  I/O's  and  tries  to  avoid  swapping  out  an 
executable  job.  Therefore,  jobs  in  long-term  wait  states 
and  other  short-term  wait  states  besides  disk  I/O's  are  the 
Job  Swapper's  primary  target. 


2. 3. 3. 3  Long-term  Wait  Queues 


The  long-term  wait  queues  hold  jobs  that  are  in  a 
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long-term  wait  state.  The  queues  and  their  purpose  are 
described  below  (Ref  28:2-3): 


CMQ 


TIOWQ 

JDCQ 


EWQ 


Command  Wait  Queue.  The  user  has  typed  a 
monitor  command  that  cannot  be  executed  until 
the  command  program  is  in  memory,  and  the 
program  is  not  in  memory. 

Teletype  I/O  Wait  Queue.  Waiting  for  the  user 
to  type  a  response  or  waiting  for  the  device 
to  print  output  already  sent  to  it. 

DAEMON  Wait  Queue.  The  job  is  waiting  for 
service  by  DAEMON.  The  DAEMON  is  a  system 
program  which  runs  as  a  user  job  and  performs 
various  functions  such  as  recording 
accounting  data  or  error  logging  required  by 
other  user  jobs.  It  is,  in  effect,  a 
non-resident  portion  of  the  monitor. 

Event  Wait  Queue.  This  queue  encompasses  many 
types  of  resource  allocation  waiting  lines 
such  as  waiting  for  a  magnetic  tape 
controller,  etc. 


Jobs  in  a  long-term  wait  queue  are  requeued  to  the  rear  of 
the  PQ1  when  their  long-term  wait  is  satisfied. 


Describing  the  above  as  queues  is  actually  a  misnomer. 
Having  a  queue  implies  that  some  order  exists  in  the  queue, 
i.e.  f irst-come-f irst  serve  (FCFS) ,  last-in-first-out 
(LIFO) ,  etc.  In  the  long-term  wait  queues  the  order  in 
which  they  leave  the  queue  is  not  dependent  on  their  order 
in  the  queue.  Instead,  the  order  in  which  they  leave  is 
totally  dependent  on  the  job  characteristics.  For  example, 
even  though  user  A  may  arrive  before  user  B  to  the  terminal 
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I/O  wait  queue  (TIOWQ) ,  if  user  B  responds  first,  he  leaves 
the  queue  before  user  A. 

The  main  purpose  of  having  long-term  wait  queues  is  two 

fold: 


-  Remove  jobs  from  the  processor  queues  so  that  the  Job 
Scheduler  does  not  waste  time  scanning  jobs  that  are 
not  expected  to  become  runnable  for  long  periods  of 
time. 

-  Provide  a  set  of  queues  that  the  Job  Swapper  can  scan 
first  to  find  a  job  to  swap-out  to  make  room  for  other 
jobs. 

The  important  thing  to  note  here  is  that  when  there  are  more 
jobs  than  can  fit  in  memory,  jobs  in  the  long-term  wait 
queues  will  be  swapped  from  the  in-core  chain  to  the 
out-core  chain  before  jobs  in  the  processor  queues. 

2. 3. 3. 4  Processor  Queue  Time  Slices 

Other  factors,  besides  which  queue  the  job  is  in, 
affect  the  likelihood  that  a  job  will  be  swapped.  These 
include  how  long  it  can  be  in  memory  without  becoming 
eligible  for  swap-out  and  its  position  in  the  queue.  These 
two  factors  are  determined  by  the  time  slice  assigned  to  a 
job. 


When  a  job  enters  one  of  the  processor  queues,  it  is 
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assigned  a  time  slice.  The  time  slice  consists  of  two 
components:  the  in-core  protect  time  (ICPT)  and  the  quantum 
run  time  (QRT) .  The  ICPT  and  QRT  are  actually  counts.  The 
ICPT  is  the  number  of  times  the  Scheduler  can  try  to  execute 
the  job  before  the  job  becomes  eligible  for  swap-out.  The 
QRT  is  the  number  of  times  the  Scheduler  can  execute  the  job 
before  it  is  requeued  to  the  rear  of  a  processor  queue. 

Each  time  the  Job  Scheduler  scans  a  job  to  see  if  it  is 
executable,  the  ICPT  is  decremented.  Note  that  the  Job 
Scheduler  may  not  execute  a  job  because  it  is  in  some 
short-term  wait  state,  but  the  ICPT  is  still  decremented. 
When  this  count  reaches  zero,  the  job  becomes  eligible  for 
swap  out.  The  ICPT  provides  a  mechanism  wherein  the  swapper 
is  prevented  from  immediately  swapping  out  a  job  which  has 
just  been  swapped  in,  and  in  effect  "locks"  the  job  in 
memory  until  its  ICPT  is  zero  (Ref  8) . 

Each  time  the  Job  Scheduler  chooses  a  job  to  execute, 
that  job's  QRT  is  decremented.  When  the  QRT  reaches  zero, 
the  job  either  requeues  to  a  lower  priority  queue,  e.g.  PQ1 
jobs  requeue  to  the  end  of  the  PQ2  queue,  or  the  job 
requeues  to  the  rear  of  the  queue  it  is  in,  e.g.  PQ2  and 
HPQ  jobs  requeue  to  the  rear  of  their  respective  queue. 
When  treated  in  this  way,  the  QRT  limits  the  amount  of  CPU 
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time  the  job  receives  in  a  particular  queue,  thus  providing 
a  fairness  consideration  in  assigning  CPU  time.  Also,  since 
the  Job  Swapper  usually  scans  the  processor  queues  from  back 
to  front,  jobs  that  have  most  recently  expired  their  QRT  and 
were  requeued  to  the  rear  of  the  queue  are  more  likely  to  be 
swapped  out  (Ref  8)  .  This  provides  a  fairness  consideration 
in  assigning  the  use  of  main  memory. 


PQ1  Time Slice 

The  PQ1  time  slice  is  the  amount  of  time  that  a  job 
receives  fast  interactive  response  after  being  swapped  in. 
The  QRT  of  the  PQ1  queue  is  usually  much  smaller  than  other 
processor  queues  in  order  to  provide  quick  response  to 
interactive  jobs.  The  short  QRT  will  be  sufficient  for 
interactive  jobs  before  they  become  blocked  to  a  long-term 
wait  state,  e.g.  waiting  for  terminal  response.  The  more 
CPU-intensive  jobs  will  expire  their  QRT  and  be  requeued  to 
the  lower  priority  PQ2  queue.  There,  they  will  be  assigned 
a  new,  much  larger  QRT,  but  now  have  a  lower  priority  for 
execution  and  a  greater  likelihood  of  being  swapped  out. 
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PQ2  Time  Slice 

For  PQ2  jobs,  the  parameters  for  ICPT  and  QRT  control 
the  bias  of  the  scheduler  for  throughput  versus  response  and 
for  I/O  versus  CPU.  Throughput  versus  response  is  controlled 
by  increasing  or  decreasing  the  magnitude  of  both 
parameters.  As  the  parameters  are  increased,  jobs  expire 
their  time  slices  more  slowly,  swapping  rate  decreases,  and 
less  core  is  allocated  for  swapping.  These  effects  improve 
throughput,  but  average  response  is  correspondingly  degraded 
because  interactive  jobs  wait  longer  to  swap  in.  When  you 
decrease  both  parameters  the  effect  is  reversed.  I/O  versus 
CPU  response  is  controlled  by  changing  the  ratio  of  ICPT  to 
QRT.  Increasing  only  QRT  favors  CPU-bound  jobs.  Increasing 
only  ICPT  favors  I/O-bound  jobs,  while  reducing  it  tends  to 
favor  CPU-bound  jobs  (Ref  28:5-3). 


2.4  Chapter  Summary 


The  following  are  the  important  aspects  to  remember 
about  the  TOPS-IO  Monitor 


The  Monitor  executes  every  jiffy  at  which  time  the  Job 
Scheduler  selects  a  job  to  run  the  remainder  of  the 

jiffy. 
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Jobs  in  the  system  can  be  in  one  of  four  states 

1.  Executing  on  the  CPU. 

2.  Ready  to  run  but  waiting  for  its  turn  to  execute 
on  the  CPU. 


3. 

In  a  short-term 
shareable  resource. 

wait 

state 

waiting 

for 

a 

4. 

In 

a  long-term 

wait 

state 

waiting 

for 

a 

non-shareable  resource. 


Jobs  in  the  system  wait  in  processor  queues  for  the  CPU 
or  for  shareable  resources  while  jobs  waiting  for 
non-shareable  resources  wait  in  the  long-term  wait 
queues.  Interactive  jobs  wait  for  terminal  response  in 
a  long-term  wait  queue. 

The  processor  queues  and  the  long-term  wait  queues  are 
divided  into  in-core  chains  (jobs  that  have  memory)  and 
out-core  chains  (jobs  without  memory  that  will  have  to 
be  swapped  in)  . 

Jobs  in  the  long-term  wait  queues  will  be  swapped  from 
the  in-core  chain  to  the  out-core  chain  before  jobs  in 
the  processor  queues. 

The  time  slices  assigned  to  jobs  by  the  Job  Scheduler 
are  made  up  of  two  components:  the  incore  protect  time 
(ICPT)  and  the  quantum  run  time  (QRT) . 

*  The  ICPT  provides  a  mechanism  where  in  the  swapper 
is  prevented  from  immediately  swapping  out  a  job 
which  has  just  been  swapped  in,  and  in  effect 
"locks"  the  job  in  memory  until  its  ICPT  is  zero. 

*  The  QRT  limits  the  amount  of  CPU  time  the  job 
receives  in  a  particular  queue,  thus  providing  a 
fairness  consideration  in  assigning  CPU  time. 

-  The  short  QRT  in  PQ1  will  be  sufficient  for 
interactive  jobs  before  they  become  blocked 
to  a  long-term  wait  state  while  more 
CPU-intensive  jobs  will  expire  their  QRT  and 
be  requeued  to  PQ2  where  they  will  be 
assigned  a  new,  much  larger  QRT. 


2-18 


Operating  Systems  and  CPE 


The  longer  QRT  in  PQ2  provides  the  extensive 
processing  needs  of  jobs  needing  large 
amounts  of  CPU  processing.  The  ratio  of  the 
PQ2  QRT  and  the  ICPT  determines  the 
schedulers  bias  towards  CPU-intensive  jobs 
versus  I/O  intensive  jobs. 
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3.1  Earlv  Queueing  Network  Models 


Jackson  (Ref  14)  and  Gordon  and  Newell  (Ref  10)  give 
solutions  for  the  case  of  networks  of  queues  with  the 
following  assumptions 

1.  There  are  a  finite  number  of  nodes,  M. 

2.  There  are  a  fixed  finite  number  of  customers,  N. 

3.  The  manner  in  which  customers  visit  the  various 
resources  is  governed  by  a  transition  matrix  P  = 
[p..J,  where  p-.  is  the  probability  that  a  customer 
departing  from  node  i  will  next  visit  node  j. 

4.  The  service  time  distributions  at  nodes  are 
exponentially  distributed. 

5.  The  service  rate  may  be  a  function  of  the  number  of 
customers. 

6.  All  customers  are  identical  in  that  their  routing 
probabilities  and  service  times  at  the  M  nodes  are  the 
same. 


For  this  class  of  models,  a  state  of  the  model  is  given 
by  the  number  of  customers  at  each  node.  Thus  the  set  of 
states  is  just 


i 

i 
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M 


'  n .  f  /  •  *  •  /  n>.)  I  n  ■ 

12  M  i=l  1 


N 


The  equilibrium  state  probabilities  are  given  by 


P(nlf  n2,...nM) 


M 

c  n 

i»l 


n. 


ni 


n  u. 


i(j> 


where 

X  The  mean  arrival  rate  to  the  ith  node. 

u,(j)  The  instantaneous  departure  rate  from  the  ith 

1  resource  when  there  are  j  customers  queued  at 

this  node. 

C  The  normalization  constant  is  chosen  so  that 

all  the  equilibrium  state  probabilities  sum 
to  one. 


The  calculation  of  C  can  be  very  time  consuming  since 
the  number  of  states  is 


However,  methods  of  calculating  C  have  been  found  which 

2 

increase  only  as  MN  (Ref  3}  allowing  larger  queueing 
networks  to  be  solved. 
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.2  Multi-class  Queueing  Network  Model 

The  earliest  queueing  models  developed  by  Jackson  and 
Gordon  and  Newell,  and  the  computational  algorithm  derived 
by  Buzen  suffered  from  two  limitations: 

1.  Only  one  class  of  customers  was  allowed  in  the 
network. 

2.  All  the  service  time  distributions  were  exponential. 

Recently  developed  queueing  models  addressed  many  of 
the  important  aspects  of  applying  queueing  models  to 
time-sharing  computer  systems  such  as  multi-class  jobs  and 
non-exponential  service  time  distributions.  Still,  no  one 
unified  model  existed  to  describe  all  these  aspects  of 
time-sharing  systems. 

The  model  developed  in  reference  2  by  Basket,  Chandy, 
Muntz,  and  Palacios  and  described  below  accomplished  this 
unification.  It  combined  recent  results  on  the  networks  of 
queues  of  several  different  service  disciplines  and  a  broad 
class  of  service  time  distributions  with  earlier  results  on 
networks  of  queues  containing  different  classes  of 
customers.  This  model  thus  more  accurately  describes  a 
time-shared  computer  system. 
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3.2.1  Variables  in  the  Model 

The  model  describes  a  system  with  an  arbitrary  but 
finite  number  (N)  of  service  centers  with  an  arbitrary  but 
finite  number  (R)  of  job  classes.  The  routes  through  the 
network  of  service  centers  for  the  R  classes  of  jobs  is 
determined  by  a  probability  transition  matrix 
P  *  [P4  r.+  «J  *  where  P.  .  is  the  probability  of  a  class 
r  job  at  service  center  i  going  to  service  center  j  as  a 
class  s  job.  The  transition  matrix  defines  a  Markov  chain 
assumed  to  be  decomposable  into  m  subchains,  E.,E,,...E  . 

Let  n^r  be  the  number  of  jobs  of  class  r  at  service 
center  i  in  state  Sj  of  the  network  model.  Also  let 

M(S/E . )  -  2  n. 

2  (i,r)  €  E j  ir 

be  the  total  number  of  jobs  within  the  subchain.  Then 


M(S) 


m 


E 

i«l 


M (S/E j ) 
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is  the  total  number  of  jobs  in  the  network.  A  system  is 
closed  when  M(S/Ej)  *  constant,  1  <  j  <  m. 

3.2.2  Open  Network  Outside  Arrival  Processes 

In  an  open  network  the  arrivals  to  node  in  a  queueing 
network  can  be  external  from  the  network.  The  fixed 
probability  of  an  outside  arrival  of  a  class  r  customer  at 
service  center  i  is  <3ir*  The  probability  of  a  class  r  job 
leaving  the  system  from  service  enter  i  is 


1  - 


E 

1<  j<N 
1<S<R 


pi , r ;  j ,  s 


The  external  arrival  process  can  be  state  dependent  and 
is  of  two  general  types 


1.  The  arrival  rate,  X(M(s))  is  dependent  on  the  total 
number  of  customers  in  the  system,  M(S) . 

2.  The  arrival  rate,  \j(M(S/E.)  is  decomposed  into  m 
Poisson  arrival  streams  corresponding  to  the  m 
subchains  described  above.  The  m  arrival  rates  are 
dependent  on  the  number  of  customers  in  the  subchain, 
M(S/E..) . 
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3.2.3  Service  Time  Distributions 

Exponential,  hyperexponential,  and  hypoexponential 
distributions  are  eligible  service  time  distribution  for 
this  model.  Also,  any  service  time  distribution  that  can  be 
represented  as  a  network  of  states  can  be  used  as  a  service 
time  distribution.  In  this  method  of  stages  technique, 
represents  the  service  rate  of  a  job  in  service  stage  1  at 
service  center  i  who  is  in  class  r.  There  are  s^r  service 
stages  for  a  job  in  class  r  at  node  i.  The  probability  of 
going  to  the  next  stage  is  a^^,  while  the  probability  of 
completing  service  is  birl<  These  concepts  are  illustrated 
in  Figure  3-1. 


Figure  3-1  Method  of  Stages. 
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3.2.4  State-dependent  Service  Rates 


Three  types  of  state-dependent  service  rates  can  be 
incorporated  into  the  Basket,  et.  al.  model 

1.  The  service  rate  at  a  service  center  can  depend  on  the 
total  number  of  jobs  at  that  service  center.  This 
form  of  state-dependent  service  rate  is  useful  to 
model  multiple,  identical  servers.  If  there  are  k. 
identical  servers  at  service  center  i,  then  by  taking 
the  service  rate  with  one  job  present  and  multiplying 
it  by  the  following  function 


ni#  1  £  ni  <  k, 

k,  n^^  >,  k 

we  obtain  the  service  rate  for  service  center  i. 

2.  The  service  rate  of  a  job  class  can  depend  at  a 
service  center  on  the  number  n.  of  class  r  jobs  at 
service  center  i. 

3.  The  service  rate  of  a  service  center  can  depend  on  the 
number  of  jobs  at  other  service  centers  in  the 
network. 


Note  that  these  various  forms  of  state-dependent  service 
rates  can  be  mixed. 


xi(ni) 


3.2.5  Service  Center  Types 

Service  centers  can  be  one  of  four  types: 

1.  The  service  discipline  is  f irst-come-f irst-served 
(FCPS) . 
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-  All  customers  have  the  same  service  time 
distribution  at  this  service  center. 

-  The  service  time  distribution  must  be  a  negative 
exponential. 

-  Only  state-dependent  service  rate  types  1  and  3 
can  be  used. 

This  type  of  service  center  is  most  often  used  to 
model  the  I/O  devices  such  as  disks,  drums,  and  tape 
drives  in  a  computer  system. 

2.  The  service  discipline  is  processor  sharing  (PS), 
i.e.  when  there  are  n.  customers  in  the  service 
center,  each  is  receiving  service  at  a  rate  of  1/n. 
times  the  service  center's  rate  for  its  job  class. 

Each  class  of  customer  can  have  a  distinct 
service  time  distribution. 

-  All  three  service  time  distributions  can  be 
used. 

-  All  three  state-dependent  service  rate  types  can 
be  used. 

This  type  of  service  center  is  most  often  used  to 
model  a  central  processor  that  uses  a  round-robin 
scheduling  algorithm.  This  is  appropriate  since 
round-robin  scheduling  approaches  the  processor 
sharing  service  time  as  the  round-robin  time  quantum 
approaches  zero. 

3.  The  number  of  servers  in  the  service  center  is  greater 
than  or  equal  to  the  maximum  number  of  customers  that 
can  be  at  this  center  in  a  feasible  state  (this  is  the 
infinite  server  (IS)  case). 

-  Each  class  of  customer  can  have  a  distinct 
service  time  distribution. 

-  All  three  service  time  distributions  may  be 
used. 

-  All  three  state-dependent  service  rate  types  can 
be  used. 
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This  type  of  service  center  is  most  often  used  to 
model  the  terminal  users  in  a  time-sharing  computer 
system. 

4.  The  service  discipline  is  preemptive-resume  last-come- 
first-served  (LCFS) 

-  Each  class  of  customer  can  have  a  distinct 
service  time  distribution. 

-  All  three  service  time  distributions  can  be 
used. 

-  All  three  state-dependent  service  rate  types  can 
be  used. 

This  type  of  service  center  is  most  often  used  to 
model  a  central  processor  in  which  jobs  can  preempt 
other  jobs  using  the  processor. 


3.2.6  States  of  the  Model 

The  state  of  the  model  is  represented  by  a  vector 
(x-^  ,  x2  ,...xN>  where  xi  represents  the  state  of  service 
center  i.  The  representation  of  the  state  of  the  service 
center  is  dependent  on  the  type  of  service  center. 


3. 2. 6.1  Type  1  Service  Center 

If  service  center  i  is  of  type  1,  then  x^  * 
(xu  ,xi2  ^..x^),  where  k  *  n^  is  the  number  of  customers 
at  center  i  and  x^j  (1  <  j  <  n^f  1  <  x^j  <  R)  is  the  class 
of  customer  who  is  jth  in  FCFS  order.  The  first  customer  is 
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served  while  the  remainder  are  waiting  for  service.  This 
type  of  service  center  state  space  is  very  large  since  one 
must  not  only  account  for  the  number  of  customers  of  each 
type,  but  the  order  of  the  queue  must  be  properly 
represented. 

3. 2. 6. 2  Type  2  and  Type  3  Service  Centers 

If  service  center  i  is  of  type  2  or  3,  then 

x^  *  (vil  ,  vi2  r...viR),  where  vir  is  a  vector 

th 

(mlr  ,  m2r  ,  ...  mkr;)  where  k  *  s^r.  The  1  component  of 

v.  is  the  number  of  customers  of  class  r  in  center  i  and  in 
lr 

the  1th  stage  of  service.  The  number  of  stages  for  a  class 
r  customer  at  service  center  i  is  sir.  Note  here  that  a 
state  is  distinguished  from  another  state  by  the  number  of 
customers,  the  class  of  the  customer,  and  the  stage  of 
service  of  the  customer,  but  not  the  order  of  the  customers 
at  the  service  center  since  there  is  no  waiting  line  for 
these  service  disciplines. 

3. 2. 6. 3  Type  4  Service  Centers 

If  the  service  center  is  of  type  4,  then 

xi  “  **ri'  mi*  2 '  m2* '  *•*  ^rn  fIBn^  where  the  ordered 
pair  (rj,  m..)  describes  the  customer  in  LCPS  order, 


I 
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i.e.  r j  is  the  class  of  the  customer  and  nij  is  the  stage  of 
service  of  the  customer.  Note  again  that  like  the  type  1 
service  centers,  a  queue  exists  and  one  must  take  into 
account  the  order  of  customers  in  the  queue  if  there  is  more 
than  one  class. 

3.2.7  Calculating  Equilibrium  State  Probabilities 

3. 2. 7.1  Balance  Equations 

The  balance  equation  technique  for  the  solution  of 
equilibrium  state  probabilities  is  based  on  the  concept  that 
the  rate  of  customers  transitioning  into  a  state  is  equal  to 
the  rate  of  customers  transitioning  out  of  the  state,  or 
more  formally,  for  all  states,  Sj,F 

P(Sj)  Irate  of  flow  — >S {  ]  » 
all  states  J  J 

Si 

J  P(S^) Trate  of  flow  out  of  1 

These  are  the  global  balance  equations  of  the  queueing 
network.  The  balance  equations  establish  a  set  of  linear 
equations  which  can  then  be  solved  for  the  equilibrium  state 
probabilities. 
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A  second  type  of  balance  equation  for  a  queueing 
network  exists,  that  of  independent  (local)  balance.  This 
concept  equates  the  rate  of  flow  into  a  state  by  a  customer 
entering  a  stage  of  service  to  the  flow  out  of  that  state 
due  to  the  customer  leaving  that  stage  of  service.  A 
customer  can  be  associated  with  a  stage  of  service  in  the 
following  ways 

1.  If  the  customer  is  in  service  at  a  service  center 

(always  the  case  in  type  2  or  3  service  centers) ,  then 
he  is  in  one  of  the  stages  of  his  service  time 
distribution  at  that  service  center. 

2.  If  the  customer  is  queued  at  a  service  center  (only 

possible  for  type  1  or  4  service  centers) ,  then  he  is 

in  the  stage  of  his  service  time  distribution  he  will 

enter  when  he  is  next  served.  For  FCFS  this  is  stage 
1;  for  LCFS  this  is  the  stage  of  service  when  last 
preempted. 


3. 2. 7. 2  Product  Form  Solution 


Before  presenting  the  solution  to  this  type  of  queueing 
network,  we  must  first  define  a  few  additional  terms. 


Additional  Notation 


Bach  of  the  subchains  defined  above  has  associated  with 
it  a  set  of  linear  equations  of  the  form 
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£  eir^i/r;j,s  +  9js  ‘  ejs<  ^'s)  Ek 

(i,r)  e  Ek 

where  q.  _  is  the  rate  of  exogenous  arrivals  of  class  s 
j  s 

customers  to  service  center  j.  If  <ljs*0  (j,s)  Ek,  then  the 
network  is  closed  with  respect  to  E^.  In  the  case  of  closed 
networks,  the  above  linear  equations  do  not  provide  a  unique 
solution.  By  setting  the  value  of  one  of  the  e^^  a  unique 
solution  can  be  obtained.  For  convenience,  one  of  the  eir 
is  set  to  1,  which  allows  the  rest  of  the  e^r  to  be 
interpreted  as  relative  arrival  rates  of  class  r  customers 
to  service  center  i  (relative  in  that  the  arrival  rates  eir 
times  the  rate  of  the  eir  that  was  set  to  1) .  Note  that  the 
network  may  be  closed  with  respect  to  less  than  m  of  the 
subchains  E^. 

One  more  term  that  appears  in  the  product  form  solution 
needs  to  be  defined.  This  is  the  probability  that  the  rth 
class  at  service  center  i  is  in  the  1th  stage  of  service, 
denoted  as 

1 

Airl  “  airj 

where  a^rj  is  defined  as  in  Figure  3-1. 
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The  Theorem 


For  a  network  of  service  stations  which  is  open, 
closed,  or  mixed  in  which  each  service  center  is  of  type  1, 
2,  3,  or  4,  the  equilibrium  state  probabilities  are  given  by 

P(S  *  xl,x2, *  *  * ,XN>  *  ^  d(S)  f^(Xj)  f 2 ^x2^ *  *  * f n (XN^ 

where  C  is  a  normalizing  constant  chosen  to  make  the 
equilibrium  state  probabilities  sum  to  1,  d(S)  is  a  function 
of  the  number  of  customers  in  the  system,  and  each  fi  is  a 
function  that  depends  on  the  type  of  service  center  i. 


If  service  center  i  is  of  type  1,  then 


n. 


ni 


fiui>  *  (i/ui> 1  n  teixii1 


j-i 


If  service  center  i  is  of  type  2,  then 


R  s 


ir 


ID  • 

'  V  n  n  Ieir*irl/ulrl) 

r«l  1-1 
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If  service  center  i  is  of  type  3,  then 


If  the  arrivals  to  the  system  depend  on  the  total  number  of 
jobs  in  the  system,  M(S) ,  and  the  arrivals  are  of  class  r 
and  for  center  i  according  to  fixed  probabilities  pir,  then 

M(S)-1 

d(S)  -  n  X (i) . 
i-0 

If  we  have  the  second  type  of  state-dependent  arrival 
process,  then 

m  M(S/E j ) -1 

d(S)  -  FI  II  A(i). 

j-1  i-0 

If  the  network  is  closed,  d(S)  >1. 
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3.2.8  Example  Problem 

To  illustrate  the  concepts  defined  above/  a  very  simple 
example  will  be  formulated.  First,  the  global  and 
independent  balance  equations  for  the  problem  will  be 
defined  using  the  "rate-in  »  rate-out"  concept.  Then  the 

Table  3-1  Transition  Matrix  and  Service 

Rates. 


Probability  Transition  Matrices 

by  Job 

Class 

Job  Class  1 

Job  Class  2 

0.6  0.4 

0.8  0.2 

1.0  0.0 

1.0  0.0 

Service  Rates 

by  Job  Class 

Node 

Job  Class  1 

Job  Class  2 

1 

1.0 

1.0  (processor 

sharing) 

2 

2.0 

1.0  (infinite 

server) 

N1  - 

2  N-  -  2  S.  * 

1  R  -  2 

2  lr 
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equilibrium  state  probabilities  will  be  solved  for  using  the 
product  form  solution.  These  solutions  will  be  verified  by 
substituting  the  values  in  some  of  the  balance  equations  to 
demonstrate  that  they  do  indeed  solve  the  balance 
equations. 

The  example  will  consist  of  a  closed  queueing  network 
with  two  service  centers.  There  will  be  two  classes  of  jobs 
with  two  jobs  from  each  class  in  the  network.  Service 
center  1  will  use  the  processor-sharing  service  discipline, 
while  service  center  2  will  use  the  infinite  server  service 
discipline.  This  network  is  a  very  simple  model  of  an 
interactive  time-sharing  computer  system  where  service 
center  1  represents  the  computer  and  service  center  2 
represents  the  terminal  users.  The  values  of  the 
probability  transition  matrix  and  service  rates  are  shown  in 
Table  3-1  along  with  a  summary  of  the  problem  statement. 

3. 2. 8.1  The  Balance  Equations 

The  state  of  the  system  can  be  represented  as  an 
ordered  4-tuple  *nn'n2i'ni2'n22*  where  the  first  subscript 
represents  the  service  center  and  the  second  represents  the 
class.  Therefore,  n^  is  the  number  of  jobs  at  service 
center  1  from  job  class  1. 
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Given  the  transition  matrix  and  service  rates  in  Table 
3-1,  the  state  transition  diagram  is  as  shown  in  Figure  3-2. 
Using  the  diagram  in  Figure  3-2  the  global  and  local  balance 
equations  can  be  derived.  The  global  balance  equations  of 
this  example  are  shown  in  Figure  3-3.  The  local  balance 
equations  are  shown  in  Figure  3-4. 
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P(2,0,0,2) lu22(2) (1)  +  u,,(2/ (0+2) ) (.4+. 6) ]  =P(2, 1,0,1)  u, , (1/1+2) ) ( . 2) 


Figure  3-3  Global  Balance  Equations 
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(Local)  Balance  Equations 


(1/ (2+D)  (.2+. 8)1  *  P<2,1,0,1)u12(1/<2+1>>  (.8) 


Figure  3-4  Independent  (Local)  Balance  Equations  (cont. 
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3. 2. 8. 2  The  Product  Form  Solution 


Since  there  are  no  stages  of  service  in  either  service 
center,  the  function  f^x^)  and  f2(x2>  can  be  simplified  to 
the  following 


f 1 C  xi )  ■  n1I 


2 

Jl  lelr/ulr]"lr<l/mitn 


f2(x2)  * 


2 

n 

r=l 


Ie2r/u2r”2r‘1/m2c!) 


Solving  for  the  eir  above  for  this  example,  we  get  the 
following  linear  equations  for  class  1  jobs 


ell  "  *6ell  +  e21 
*4ell  *  e21 

and  for  class  2  jobs 


e12  ‘  -8e12  +  e22 
•2e12  "  e22 

which  yield  the  following  solutions  for  the  e^r's 


i 

i 


3-24 


Analytical  Modeling 


ell 


1.0 


e21 

e12 


0.4 

1.0 


e22 


0.2 


Table  3-2  Example  Calculations 


1  State 

£i(*i 

)  f2(x2) 

P(state) 

1  (2/2/0, 0) 

3/2 

1 

1.5C  ■  .670 

1  (1/2/1, 0) 

3/4 

.4 

. 3C  »  .134 

1  (2, 1,0,1) 

3/2 

.2 

. 3C  -  .0356 

1  (1,1, 1,1) 

1 

,08 

. 08C  «  .0356 

1  (0,1, 2,1) 

1/2 

.016 

. 008C  *  .00356 

1  (1,0, 1,2) 

1 

.008 

. 008C  ■  .00356 

1  (0,0, 2, 2) 

1 

.0016 

. 0016C  «  .000715 

1  (0,2, 2,0) 

1/4 

.08 

. 02C  »  .00894 

1  (2, 0,0, 2) 

1 

• 

o 

to 

.02C  -.00894 

2.2376C  -  1  -> 

C  -  .4469 

The  results 

of  the 

computation 

are  shown  in  Table 

3-2. 

The  solutions  can  be  verified  by  substituting  them  into  some 
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of  the  balance  equations.  For  example,  substituting  the 
equilibrium  probabilities  into  the  24th  independent  balance 
equation,  we  get 


(.00894) (1) (2)  (1)  *  (.134)  (1)  (1/3) (.4) 
.01787  *  .01787 

and  for  the  17th  balance  equations 


(.003575)  (1)  (1)  (1)  =»  (.00894)  (2)  (1)  (.2) 

.003575  *  .003575 

Therefore  the  product  form  solution  of  the  above 
theorem  provides  solutions  to  the  balance  equations  of  this 
problem. 


1*1  Chen '  s_-SMflpp.ina,  Haflal 

The  classical  swapping  model  described  in  Chapter  1, 
along  with  the  inclusion  of  multi-class  jobs  and 
non-exponential  service  times,  provides  a  very  flexible 
model  of  a  computer  system  and  its  job  swapping  behavior. 
But  the  model  still  is  not  completely  realistic;  job 
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swapping  behavior  is  more  complex  and  should  depend  on  the 
main  memory  size,  the  number  of  jobs  competing  for  memory, 
and  the  job  sizes. 

Chen  (Ref  5)  proposed  a  new  approach:  that  job  swapping 
behavior  be  modeled  by  state-dependent  routing  probabilities 
in  a  closed  queueing  network  with  a  single  class  of  jobs. 
The  probability  that  a  program  needs  to  be  swapped  in  (or 
out)  is  expressed  as  a  function  of  the  system  state,  the 
individual  job  characteristics,  and  main  memory  size. 
Unfortunately,  this  type  of  queueing  model  has  no  known 
exact  solution.  Therefore,  Chen  developed  a 
successive-approximation  algorithm  that  iteratively  solves  a 
closed  queueing  network.  When  two  successive  iterations 
yield  results  that  are  within  specified  tolerances  of  each 
other,  the  swapping  model  is  solved.  Chen's  model  and 
algorithm  are  described  below. 

3.3.1  Variables  in  the  Model 

First  we  will  consider  a  simple  model  of  an  interactive 
computer  system  described  by  Chen  (Ref  5)  and  shown  in 
Figure  3-5.  This  system  consists  of  one  CPU,  one  disk,  one 
swapping  drum,  and  a  set  of  terminals.  Queues  exist  for  the 
CPU,  swapping  drum,  and  disk,  but  no  queue  is  necessary  for 
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the  terminal  node  since  there  is  a  "server"  (terminal  user) 
at  each  active  terminal. 

The  variables  in  the  model  are 

-  Number  of  jobs  in  the  system  (N) . 

-  Main  memory  size  (user  area)  (M) . 

-  Average  job  size  (J) . 

-  Average  CPU  service  time  (1/u^). 

-  Average  disk  I/O  time  (l/u2) . 

-  Average  time  to  swap  a  job  in  and  out  of  memory  (I/U3) • 

-  Average  user  think  time  (l/u^) . 

-  Average  CPU  time  needed  per  interaction  (Tcpo) . 

-  Average  number  of  disk  I/O  requests  per  interaction 

'W  • 

-  Average  job  size  (J) . 


These 

variables 

above  are 

used 

to  calculate 

the 

probability 

transition 

matrix  p 

«  p.  . 

13 

where  P^j  is 

the 

probability  a  job  being  serviced  at  node  i  will  request  the 
service  of  the  node  j.  In  terras  of  the  context  of  the 
model,  the  meaning  of  some  of  the  important  P^j's  ate  given 
below 

P.3  The  probability  that  the  job  is  not  in  main 

memory  and  needs  to  be  swapped  in. 

P41  The  probability  that  the  job  is  already  in 

main  memory  and  does  not  need  to  be  swapped 
in. 
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P,,  The  probability  that  the  job's  time  quantum 

xx  has  expired 

P12  The  probability  that  the  job  needs 

XA  information  on  a  disk  file  or  needs  to  write 

its  results  to  a  disk  file. 

P13  The  probability  that  the  main  memory  size  is 

not  sufficient  and  the  job  must  be  swapped 
out. 

The  state-dependent  probabilities  in  Chen's  algorithm  are 

Pll'  P43'  P41r  and  P13*  Note  that  P4i  and  pi3  are  always 
zero  in  the  simpler  classical  swapping  model  contained  in 
Ch -pi* i  1.  A  general  form  of  the  probability  transition 
matrix  in  Chen's  model  is  shown  in  Table  3-3. 

3. 3. 1.1  Derivation  of  P14  and  P12 

Consider  a  typical  sequence  of  events  where  'i' 
represents  the  job  leaving  the  CPU  for  an  interaction  at  the 
terminal  node  and  'c'  represents  the  job  requesting  another 
CPU  service  of  1/u^^  seconds  (on  the  average).  A  typical 
sequence  would  then  look  like 

cc* • • cicc. • • cicc* • • 

The  number  of  CPU  requests  before  an  interaction  would  then 
be  geometrically  distributed  and  the  average  number  of  CPU 
service  completions  per  interaction,  Ncp would  be 
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N 


CPU 


E  W14U  -  -  1/P14 


(1) 


Table  3-3  Probability  Transition  Matrix 
for  Chen's  Swapping  Model 


P11 

P12 

P13 

P14 

1.0 

o 

• 

o 

o 

. 

o 

o 

• 

o 

1.0 

0.0 

0.0 

0.0 

P41 

O 

• 

o 

P43 

o 

• 

o 

The  average  CPU  time  needed  per  interaction  is  equal  to  the 
product  of  the  average  number  of  times  the  job  requests  CPU 
service  times  the  service  rate  or 


TCPU  “  1/,Qi  NCPU  =  d/u1)  (1/P^4) 


therefore  we  have 


P14  “  1/(ui  TCPU} 


(2) 


It  follows  immediately  that  since  there  are  Ncpu  CPU  service 
completions  and  NDIQ  of  them  must  have  been  due  to  an  I/O 
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request,  then 

P12  *  NDIO/NCPU  "  NDIO  P14 


(3) 


3.3.2  Program  Swapping  Behavior 

The  job  swapping  behavior  is  represented  by  the 
state-dependent  routing  probabilities  P43  and  pi3*  Since 
there  are  several  memory  management  schemes  that  can  be 
implemented  in  an  operating  system,  the  expressions  for  P43 
and  P13  will  vary  from  system  to  system.  Chen  derived 
expressions  for  these  probabilities  that  are  reasonable  for 
the  TOPS-IO  Monitor  on  the  DECsystem-10.  Their  derivations 
and  the  assumptions  for  their  derivations  are  discussed 
below. 

3. 3. 2.1  Derivation  of  P43 

The  assumptions  made  for  deriving  P43  are  as  follows 
(Ref  5:955) 

1.  Each  job  in  the  main  memory  is  allocated  the  same 
amount  of  memory.  Therefore,  the  maximum  number  of 
jobs  that  can  be  allocated  in  the  main  memory 
simultaneously  is  a  constant  and  denoted  by  A. 

A  »  >M/JI 
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This  assumption  is  justified  since  main  memory  is 
partitioned  into  a  fixed  number  of  parts  and 
interactive  jobs  usually  will  fit  in  the  same  number 
of  partitions. 

2.  Among  all  jobs  in  the  main  memory,  the  jobs  in  the 
"think"  mode  have  the  highest  priority  to  be  swapped 
out  if  memory  space  is  needed  for  the  jobs  to  be 
swapped  in.  This  assumption  is  justified  since  many 
computer  systems  use  this  memory  management  strategy 
(Ref  5).  Jobs  in  the  "think"  mode  are  very  unlikely  to 
need  CPU  service  in  the  near  future.  This  is  true 
since  the  average  think  time  is  measured  in  seconds, 
while  most  other  activities  in  the  system  occur  in  the 
millisecond  or  even  microsecond  range.  In  the 

DECsy stem-10 ,  jobs  contained  in  the  TIOWQ 

correspond  to  jobs  in  the  think  mode. 


Let  n^  denote  the  number  of  jobs  at  node  i.  If  the 
total  number  of  jobs  in  the  system  is  less  than  the  number 
of  jobs  that  can  fit  in  memory  (N  <  A),  then  no  swapping 
activity  is  necessary.  In  most  cases,  though,  N  is  greater 
than  A,  and  so  N  -  A  jobs  cannot  be  in  memory.  Since  the  n4 
jobs  in  the  think  mode  have  the  highest  priority  to  be 
swapped  out  (Assumption  2) ,  these  will  be  swapped  out  until 
N-A  jobs  are  swapped  or  until  all  jobs  in  the  think  mode  are 
swapped  out.  Therefore,  the  probability  that  a  job  needs  to 
be  swapped  in  at  the  beginning  of  an  interaction,  P^,  is 
equal  to  the  probability  that  a  job  in  the  think  mode  is  not 
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in  main  memory.  Chen  uses  the  approximation 


P 


43  * 


(N-A)/n4  n4  >  N  -  A  and  N  >  A 

.  0  N  ^  A 

1  n4  ^  N  -  A 

(4) 


and  therefore  the  probability  that 
swapped  in  at  the  beginning  of  the 


a  job  does  not  need  to  be 
interaction  is 


(5) 


3. 3. 2. 2  Derivation  of  Pi3 


To  derive  P13,  two  more  assumptions  must  be  made  (Ref 
5:955) 


1.  P.3  is  linearly  related  to  the  difference  between 

number  of  active  jobs  (n.+n2+n3)  and  the  maximum 
number  of  jobs  (A)  that  cah  fit-3  in  memory,  if  that 
difference  is  positive.  The  difference,  n,+n2+n3-A, 
is  the  number  of  jobs  which  need  memory  but  can  not 
have  it.  When  this  difference  is  positive,  then  some 
of  these  jobs  as  well,  as  the  jobs  in  the  think  mode, 
must  be  swapped  out.  When  this  difference  is  less 
than  or  equal  to  zero,  P.3  is  zero.  This  models  the 
TOPS-IO  Monitors  swapping1behavior  of  choosing  jobs  in 
long-term  wait  queues,  e.g.  interactive  jobs  waiting 
form  terminal  response  or  batch  jobs  waiting  for  a 
tape  to  be  mounted,  over  jobs  that  are  in  the 
processing  queues  and  more  likely  to  run. 
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2.  In  the  worst  case, the  maximum  number  of  times  that  a 
job  will  be  swapped  out  before  the  end  of  the 
interaction  is  one.  This  assumption  represents  the 
ICPT's  effect  of  locking  the  job  in  memory 


Assumption  1  tells  that  P^  reaches  a  maximum  when 
nl+n2+n3  reaches  its  maximum  value  N.  This  maximum  value  can 
be  determined  using  similar  reasoning  used  in  the  derivation 
of  Equation  (3) .  By  assumption  2  the  maximum  number  of  times 
a  job  can  be  swapped  out  is  one.  So  the  maximum  probability 
of  being  swapped  out  is 

P13  =  NSWAP/NCPU  *  1/(1/P14)  =  P14 

Therefore,  from  assumption  1,  P13  varies  linearly  from  0  to 
P14  as  n^n^^  varies  from  A  to  N.  Noting  that  P^,  P12, 
P13,  and  P14  must  sum  to  1,  thus  being  bounded  by  1-P12-P14, 
and  combining  the  above  observations,  we  have 


P 


13 


r 

min[l-P12-P14,  P14(  (n^nj+nj-AJ/tN-A) )  ] 
nl+n2+n3  >  A  and  N  >  A 

(6) 

0  otherwise 


and  therefore 

P11  *  1"P13"P12"P14 
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to  insure  the  probabilities  sum  to  one. 


3. 3. 2. 3  Approximation  Algorithm 

The  above  formulation  makes  the  probability  transition 
matrix  of  Table  3-3  dependent  on  n^,  i.e.  the  probability 
transition  matrix  is  state-dependent.  Unfortunately,  this 
queueing  network  model  has  no  exact  solution  known,  so  Chen 
proposes  using  the  following  approximation  algorithm  (Ref 


5:956) 

Step 

1 

Assume  an  initial  value  for  n4 

Step 

2 

Use  n.  to  calculate  PA,  and  by  Equations 

(4)  and  (6) . 

Step 

3 

Solve  the  model  as  a  closed  queueing  network 
with  fixed  probabilities  using  Buzen's 
method. 

Step 

4 

If  the  new  value  of  n.  is  very  close  to  its 
previous  value,  the  algorithm  stops. 
Otherwise,  Steps  2,  3,  and  4  are  repeated. 

According  to  Chen,  the  algorithm  will  always  converge 
since  the  sequence  of  computed  values  of  n^  is  either 
monotonically  increasing  or  decreasing  and  i s  bounded 
between  0  and  N.  Also,  from  Chen's  computational  experience, 
the  algorithm  converges  to  the  same  point,  independent  of 
its  starting  value. 
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3.4  Chapter  Summary 


The  following  items  are  the  important  points  brought 
out  in  this  chapter 

-  Early  models  developed  by  Gordon  and  Newell,  Jackson, 

and  Buzen  where  restricted  to  single  class  networks 
with  exponential  servers.  Baskett,  et.  al. 

consolidated  later  models  which  included  multi-class 
networks  and  the  method  of  stages  to  form  one  unified 
model. 

-  The  Baskett,  et.  al.  model  provided  a  more  accurate 
model  of  time-sharing  systems,  but  the  model  still  was 
not  completely  realistic;  job  swapping  behavior  is  more 
complex  and  should  depend  on  the  main  memory  size,  the 
number  of  jobs  competing  for  memory,  and  the  job 
sizes. 

-  Chen  developed  a  model  to  make  job  swapping  behavior 
depend  on  memory  size,  the  number  of  jobs  competing  for 
memory,  and  the  jobsizes,  but  failed  to  include  the 
more  general  multi-class  queueing  networks.  Also,  his 
development  only  considered  the  single  CPU  and  single 
disk  system,  and  left  the  multiple  CPU  and  multiple 
disk  system  as  an  extension. 


In  the  next  chapter  we  will  look  at  the  structure  of 
the  McKenzie  FORTRAN  program  and  how  it  was  modified  to 
include  Chen’s  algorithm.  Combining  McKenzie's  program  with 
Chen's  algorithm  extended  the  Chen  algorithm  to  multiple 
classes. 
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Computer  Implementation 


The  computational  effort  for  the  product  form  solution 
of  a  queueing  network  can  be  quite  large  for  anything  but 
very  simple  networks  (see  example  in  chapter  3) .  Therefore, 
to  compute  the  solutions  to  queueing  networks  containing 
many  nodes  and  job  classes,  one  must  perform  the 
computations  for  the  product  form  solution  on  a  computer. 

The  first  section  in  this  chapter  describes  the 
McKenzie  computer  program.  This  is  followed  by  a  section 
describing  the  modifications  made  to  McKenzie's  program  to 
implement  Chen's  swapping  model.  After  the  Chen 
modifications  are  described,  the  next  section  discusses  some 
of  the  pitfalls  and  problems  encountered  in  this  current 
effort  to  modify  the  program.  The  last  section  summarizes 
the  chapter. 
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4.1  The  McKenzifi_Prflg.r.aro 

In  1977  McKenzie  developed  a  computer  program  to  solve 
multi-class  closed  queueing  networks  with  service  centers 
having  different  service  disciplines.  This  program  provides 
a  flexible  tool  to  model  interactive  computer  systems,  but 
is  limited  to  product  form  solutions  of  closed  queueing 
networks.  The  program  was  modified  to  include  Chen's  model 
(see  Chapter  3)  to  more  accurately  model  job-swapping 
behavior. 

4.1.1  Program  Capabilities 

The  computer  program  is  based  on  the  model  of  Baskett, 
et.  al.  (Ref  2)  and  the  computational  algorithms  of  Wong 
(Ref  30).  Wong's  computational  algorithms  implement  a 
subset  of  the  entire  Basket  et.  al.  theory.  The  limitations 
of  the  model  are 

-  No  outside  arrivals  can  enter  the  network,  i.e.  the 
model  is  a  closed  queueing  network. 

-  Jobs  cannot  transition  from  one  job  class  to  another  at 
a  service  center. 

-  The  third  type  of  state-dependent  service  rate  is  not 
contained  in  the  model.  (see  Chapter  3  for  the  types 
of  state-dependent  service  rates) 
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4.1.2  Program  Structure 


McKenzie's  computer  program  was  written  in  FORTRAN  IV. 

The  program  is  modular  in  nature.  It  is  divided  into  five 

subroutines  called  from  the  main  program.  The  names  of  the 

subroutines  and  their  functions  are  described  below 

INDATA  Reads  the  data  necessary  to  run  the  program. 

The  characteristics  of  the  time-sharing 
computer  system  being  modeled  are  input  in 
this  routine. 

FUNCT  Calculates  the  values  for  the  functions  in 

the  product  form  solution  for  each  node  in 
the  network.  These  are  used  to  calculate  the 
joint  and  marginal  probability  density 
functions  for  each  node. 

NORMAL  Calculates  the  normalization  constant  for  the 

™odel.  The  normalization  constant  assures 
that  the  system  state  probabilities  sum  to 
one. 

MARGIN  Calculates  the  marginal  probability  density 

functions  for  each  node  in  the  network. 

EXPECT  Calculates  the  probability  that  there  are  k 

users  of  class  j  at  a  node  and  the  total 
probability  over  all  classes  of  jobs  that 
there  are  k  users  at  a  node.  Using  these 
probabilities,  it  then  calculates  node 
utilizations  by  job  class,  expected  queue 
length  probability  distribution  and  their 
expected  values  by  job  class,  mean  time  spent 
at  each  node  by  each  job  class,  and  system 
response  times  by  job  class. 


In  addition  to  these  programs  written  by  McKenzie,  five 

other  routines  peculiar  to  the  CDC  computer  at  AFIT  are  used 

WRITMS  A  random  access  routine  to  write  records  to 

mass  storage. 
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READMS 

A  random  access  routine  to  read  records  from 
mass  storage. 

OPENMS 

Initializes  a  file  for  random  record 
and  writes. 

reads 

CLOSEMS 

Updates  the  master  index  in  the  file 
all  records  have  been  read  or  written. 

after 

LEQT1F 

A  subroutine  that  solves  simultaneous 
equations. 

linear 

The  documentation  on  the  random  record  routines  is  contained 
in  reference  7.  LEQ1TF  documentation  is  contained  in 
reference  13. 

The  WRITMS  is  used  at  the  end  Of  INDATA,  FUNCT,  NORMAL, 
and  MARGIN  to  save  the  computed  values  for  each  node  on  a 
random  access  file  for  further  processing  by  the  next 
subroutine.  The  READMS  is  used  at  the  beginning  of  FUNCT, 
NORMAL,  MARGIN,  and  EXPECT  to  retrieve  the  results  of  the 
calculations  for  each  node  from  the  preceding  routine  in  the 
program  flow  (see  Figure  4-1) .  OPENMS  and  CLOSEMS  are  used 
at  the  very  beginning  of  the  program  and  the  very  end  of  the 
program.  The  calculations  are  done  node  by  node  since  the 
memory  requirements  necessary  to  store  the  entire  network's 
state  space  would  exceed  the  CDC  computers  memory  resources 
for  all  except  trivially  small  networks. 

The  subroutine  LEQ1TF  is  used  in  INDATA  to  compute  the 
node  departure  rates  using  the  probability  transition 
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matrices  and  node  service  times  for  each  class. 
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Figure  4-1  McKenzie  Program  Flow  Chart 
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4.1.3  Program  Inputs 


The  program  inputs  and  formats  are  described  in  detail 
in  Appendix  C  of  McKenzie's  thesis  and  will  not  be  discussed 
here. 


4.2  The,  Chen  Modification 


4.2.1  Program  Structure 


To  implement  Chen's  algorithm,  McKenzie's  program  had 
to  be  modified  to 

1.  Read  in  the  system  configuration  data  necessary  to  use 
Chen's  algorithm. 

2.  Compute  the  initial  probability  transition  matrix  for 
each  class  of  job  using  the  input  data  and  Chen's 
formulae  discussed  in  Chapter  3.  The  initial 
probability  transition  matrix  is  computed  using  Chen's 
formulae  and  then  converted  to  the  transition  matrix 
for  the  classical  swapping  model. 

3.  Recompute  the  probability  transition  matrix  for  each 
class  based  on  results  from  the  previous  iteration. 
Use  this  new  probability  transition  matrix  to  compute 
the  product  form  solution  to  the  corresponding 
queueing  model. 

4.  Stop  iterating  if  the  solution  has  converged,  i.e. 
the  absolute  difference  between  the  number  of  jobs  in 
the  think  mode  calculated  in  the  previous  iteration 
compared  to  the  current  iteration  is  within  a 
specified  tolerance,  or  if  the  maximum  number  of 
iterations  allowed  has  been  exceeded. 


4-6 


Computer  Implementation 


These  modifications  were  localized  to  the  main  program 
and  two  subroutines  in  order  to  retain  as  much  of  McKenzie's 
structure  as  possible.  INDATA  was  almost  completely 
rewritten  in  order  to  read  in  the  variables  in  Chen's 
model.  One  major  benefit  in  this  modification  is  that  input 
to  the  model  is  now  read  in  freefield  format,  i.e.  data  does 
not  have  to  appear  in  certain  columns  in  the  data  file. 
This  will  decrease  the  chance  of  human  error  when  typing  in 
the  input  to  the  model. 

EXPECT  was  slightly  modified  in  order  to  stop  it  from 
printing  the  results  of  every  iteration.  Now  the  results 
are  printed  for  the  first  iteration,  every  fifth  iteration, 
and  then  the  final  iteration. 

The  major  changes  to  the  McKenzie  program  to  include 
Chen's  algorithm  were  incorporated  in  the  subroutines  CHEN75 
and  COMPUTE.  At  the  beginning  of  each  iteration,  CHEN75  is 
called,  which  in  turn  calls  COMPUTE  to  recompute  the 
probability  transition  matrix  for  each  class  of  job  in  the 
network  using  the  formulae  from  Chen's  model. 

Chen's  model  was  also  extended  to  include  multiple  disk 
drives.  Given  that  a  job  from  a  particular  class  needs  a 
disk  I/O,  one  must  define  the  disk  access  probabilities  for 
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each  disk  in  the  system.  Defining  ndISKS  to  be  the  number 
of  disks  in  the  system,  and  the  probability  of 
performing  the  disk  I/O  on  the  ifc^  disk,  then  the 
probability  transition  matrix  is  modified  as  shown  below 
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To  conform  to  normal  probability  transition  matrix  notation, 
the  subscripts  are  renumbered.  For  example,  in  the  case  of 
a  system  with  six  disks,  the  subscripts  are  changed  so  that 
P13  becomes  P18,  P14  becomes  P^g,  P41  becomes  Pg^,  etc. 

The  changes  described  lbove  are  illustrated  in  the  new 
flow  chart  for  the  program  shown  in  Figure  4-2. 

4.2.2  Program  Inputs 

The  input  to  the  new  program  is  simpler  now,  since  all 
input  data  values  are  read  in  freefield,  instead  of 
requiring  the  values  to  be  in  particular  columns.  All  that 
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Figure  4-2  New  Program  Flow  Chart 
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is  necessary  is  that  the  input  data  values  be  separated  by 
at  least  one  space.  The  data  could  be  entered  on  one  line 
until  that  line  is  filled  up  and  continue  in  that  fashion 
for  subsequent  lines.  Instead/  it  is  recommended  that  the 
data  be  entered  in  the  following  manner  (the  names  in 
capital  letters  are  the  FORTRAN  variable  names) 

-  Network  description.  Enter  the  number  of  nodes  and 
number  of  job  classes  in  the  network. 

NODES,  I TYPES 

-  Computer  system  configuration.  Enter  the  number  of 
disk  drives  in  the  I/O  subsystem  and  the  size  of  main 
memory  available  for  user  jobs. 


NDISKS,  SIZEMEM 

Workload  configuration.  If  possible/  enter  the 
following  workload  configuration  parameters  on  separate 
lines 

1.  Disk  access  probabilities  by  class  for  each 
disk.  If  possible,  put  the  probabilities  for 
each  disk  on  separate  lines.  There  would  then  be 
NDISKS  lines. 


DISKPRB( I/J)  (probability  that  a  member 
of  the  j  class 
accesses  the  i  1  disk) 


2.  Job  class  characteristics.  Put  the 

characteristics  for  each  class  on  a  separate 
line.  This  includes  the  number  of  jobs  in  the 

class,  the  CPU  time  per  interaction  (T^p^) ,  the 

number  of  disk  I/O's  per  interaction  (N^Iq) ,  and 

average  job  size  of  the  class.  There  should  be 
ITYPES  lines  of  job  class  characteristics. 
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(NUSERS ( J) ,  TCPU(J),  DIO ( J) ,  SIZEJOB(J)) 

Node  characteristics.  The  node  characteristics  include 
the  following 

1.  Node  type  (FCFS  (NODETYP(I) *1)  ,  PS 

(NODETYP ( I) “2) ,  IS  (NODETYP(I) »3) ,  and  LCFS 

(NODETYP( I) =4) ) . 

2.  Types  of  state-dependent  service  rates  can  be 
None  ( IDEP (I) *1) ,  number  of  total  jobs  at  node 
(IDEP(I)=2),  number  of  jobs  in  a  particular  class 
at  the  node (IDEP (I) *3) . 

3.  Service  rate  for  job  class  J  and  node  I 
(SERVICE ( I, J) ) . 

4.  State-dependent  service  rates.  If  IDEP(I)  was 
not  1,  then  the  state-dependent  service  rate  data 
must  be  entered.  The  entered  values  represent 
the  ratio  of  the  service  time  when  there  is  a 
total  of  one  job  at  the  node  (IDEP(I)=2)  or  when 
there  is  one  job  of  a  particular  class 
(IDEP (I) =3)  .  For  example,  to  model  two  identical 
servers,  one  would  set  IDEP (I) *2  and  enter 


1.0  2.0  2.0  ...  2.0 

into  array  (DEP(I,1,K) ,K*1,MAXUSER)  The  number  of 
entries  is  equal  to  the  maximum  number  of  users 
of  all  types  in  the  network  (MAXUSER) . 


.3  Programming  Notes  for  Future  Modification 

If  future  thesis  students  decide  to  work  with  the 
program  used  in  this  thesis,  the  student  should  be  aware  of 
some  possible  pitfalls.  The  problems  encountered  when 
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making  changes  to  the  program  are  described  in  the  following 
paragraphs. 

The  original  McKenzie  program  depended  on  the  operating 
system  to  initialize  the  variables  in  the  model  to  zero. 
Unfortunately,  the  CDC  NOS/BE  operating  system  no  longer 
zeros  out  the  memory  locations  before  loading  a  program,  but 
instead  sets  them  to  a  special  value  signifying  an 
indefinite  value.  In  the  process  of  modifying  the  program, 
all  uninitialized  variables  were  initialized  to  zero  for 
each  iteration.  Some  variables  may  have  slipped  by,  though, 
so  if  any  "Error  Mode  4—  Indefinite  Value"  messages  appear 
after  modifying  the  program,  check  for  uninitialized 
variables. 

In  order  to  have  a  variably  dimensioned  array  passed  to 
a  FORTRAN  IV  subroutine,  one  must  not  only  pass  the  variable 
dimension  value,  but  also  the  memory  space  (in  the  form  of 
an  array  from  the  main  program) .  That  is  the  reason  the 
subroutine  parameter  lists  are  so  large.  Unfortunately,  if 
one  passes  incorrect  dimensions  of  a  passed  array  to  the 
subroutine,  no  syntax  error  is  generated.  Instead,  when  the 
program  is  run,  it  will  produce  one  of  possibly  several  mode 
errors.  Therefore,  when  making  any  changes  to  the 
dimensions  of  the  arrays  in  the  main  program,  be  sure  to 
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change  the  appropriate  variables  (MAXNODE,  MAXTYPE,  MAXUSER, 
MAX ST AT,  MAXUSE1 ,  and  ISIZE3)  to  match  the  true  dimensions. 

Finally,  one  must  be  careful  that  the  inputs  for  the 
model  are  consistent.  Specifically,  a  common  and  somewhat 
obscure  error  is  to  define  a  probability  transition  matrix 
that  unintentionally  isolates  a  node  from  the  network.  This 
causes  the  departure  rate  from  that  node  to  be  zero.  When 
the  subroutine  FUNCT  is  called  to  compute  the  functions  used 
in  the  product  form  solution  for  that  transition  matrix, 
that  zero  value  is  used  as  a  base  in  an  exponential 
expression  (line  34)  which  can  have  a  zero  exponent.  Since 
zero  to  the  zero  power  is  undefined,  another  mode  error  is 
produced;  this  time  "Error  Mode  2 — infinite  operand". 
Therefore,  always  make  sure  the  probability  transition 
matrix  used  in  the  computations  is  well-defined. 


A«.4-  Chap&ej  Summary 

In  this  chapter,  we  have  reviewed  the  basic  structure 
of  the  McKenzie  program  and  the  functions  of  its  five 
subroutines.  The  modifications  made  to  this  program  were 
identified,  and  the  new  subroutines  CHEN75  and  COMPUTE 
discussed.  The  modifications  were  made  with  the  objective 
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of  retaining  as  much  of  the  original  McKenzie  model  as 
possible.  All  of  the  changes  to  the  original  program  were 
in  the  main  program,  INDATA,  and  EXPECT.  The  addition  of 
Chen's  algorithm  was  accomplished  by  including  CHEN75  and 
COMPUTE.  Finally,  some  programming  notes  on  this  author's 
experience  with  modifying  the  program  were  presented. 
Briefly,  these  included 

1.  If  any  "Error  Mode  4 —  Indefinite  Value"  messages 
appear  after  modifying  the  program,  check  for 
uninitialized  variables. 

2.  When  making  any  changes  to  the  dimensions  of  the 
arrays  in  the  main  program,  be  sure  to  change  the 
appropriate  variables  (MAXNODE,  MAXTYPE,  MAXUSER, 
MAXSTAT,  MAXUSE1 ,  and  ISIZE3)  to  match  the  true 
dimensions. 

3.  Always  make  sure  that  the  probability  transition 
matrix  used  in  the  computations  is  well-defined. 

With  this  modification  to  McKenzie's  original  program, 
we  now  have  a  new  tool  to  model  interactive  computer  systems 
which  considers  multi-class  job-swapping  behavior  as  a 
function  of  main  memory  size,  the  number  of  jobs  competing 
for  memory,  and  the  job  size.  Chapter  5  presents  the 
results  of  using  the  multi-class  Chen  swapping  model  and  the 
classical  swapping  model  on  three  hypothetical  workload 
models. 


Chapter  5 


Analytical  Modeling  Results 


This  chapter  contains  the  results  of  modeling  the 
DECsystem-10  using  the  single-class  and  multi-class  Chen 
swapping  model  on  a  realistic,  though  hypothetical 
interactive/batch  workload  configuration.  The  workload  was 
modeled  using  three  alternate  approaches. 

The  first  section  discusses  the  three  alternate 
approaches  used  to  model  the  workload  configuration.  The 
next  section  describes  the  system  configuration  used  to 
model  the  DECsystem-10.  The  third  section  presents  the 
results  of  two  model  comparisons:  Chen's  swapping  model 
versus  the  classical  model  and  the  multi-class  versus  the 
single-class  Chen  swapping  model. 


5.1  Modeling  Interactive/Batch  Workloads 

When  a  CPE  analyst  is  faced  with  the  task  of  analyzing 
the  performance  of  an  interactive  computer  system  containing 
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batch  jobs  with  an  analytical  model,  he  has  the  following 
alternatives 

1.  Characterize  the  interactive  and  batch  job  classes 
separately  and  use  a  multi-class  queueing  model. 

2.  Use  a  single  class  model  and  make  one  of  following 
commonly  used  approximations 

-  Create  a  single,  average  job  class  by  aggregating 
the  characteristics  of  the  individual  classes. 

-  Assume  the  impact  of  batch  jobs  is  minimal,  and 
use  a  single  class  with  the  characteristics  of 
the  interactive  jobs.  This  approach  was  used  by 
Chen  in  his  original  work. 

5,2  DSCgystem-lQ  System. Configuration.  Parameters 

Because  of  the  large  number  of  system  configuration 
parameters  which  can  be  varied  within  the  models,  an 
infinite  variety  of  system  configurations  could  be 
analyzed.  A  subset  of  them  were  arbitrarily  fixed  in  order 
to  compare  modeling  approaches.  These  configuration 
parameters  can  be  divided  into  two  types:  hardware  and 
workload. 

5.2.1  Hardware  Parameters 

The  hardware  configuration  parameter  values  were  chosen 
to  model  as  closely  as  possible  the  Avionics  lab 
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DECsystem-10  hardware  configuration.  These  hardware 
parameters  are  given  in  Table  5-1.  The  hardware 
configuration  includes  two  CPU's,  six  disk  drives  within  the 
I/O  subsystem  and  one  swapping  device. 

Table  5-1  System  Configuration  Parameters 


Number  of  Disks  *  6 
Number  of  CPU's  =  2 
Number  of  Swapping  Devices  »  1 


u 

iteractlve 

Number  of  Jobs 

30 

10 

Job  Size  (K) 

5 

10 

CPU  quantum 
(sec) 

Disk  Access 

.0125 

.015 

Probabilities 

Equal 

Equal 

Disk  Service 
Times  (sec) 
Long-term  wait 

.075 

.075 

Service  (sec) 
Job  Swap  Time 

8.0 

1.0 

(sec) 

.1 

.2 

Tcpo  <8.0 

DIO 

.1 

4.0 

.5 

5.0 
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5.2.2  Workload  Parameters 


The  workload  configuration  parameter  values  shown  in 
Table  5-1  were  selected  to  model  typical  interactive  and 
batch  jobs.  These  values  are  within  the  range  of  values 
used  by  Chen  to  assure  their  reasonableness. 


Additional  knowledge  of  the  DECsystem-10  and  of 
interactive  and  batch  jobs  in  general  were  used  to  set  the 
relationships  between  the  magnitudes  of  the  workload 
parameters.  These  are  described  below  for  each  workload 
parameter. 

Job  Size  The  above  job  sizes  are  within  the  range  of 

job  sizes  used  by  Chen.  The  interactive  job 
size  was  chosen  smaller  than  the  batch  jobs 
since  this  is  true  on  the  average. 
Interactive  jobs  include  many  editing , 
compiling,  and  other  utility  programs  which 
have  small  memory  requirements.  Batch  jobs, 
on  the  other  hand,  include  the  running  of 
large  application  packages  such  as  SPSS, 
SLAM,  etc. 

CPU  Quantum  This  value  corresponds  to  the  jiffy  described 
in  Chapter  2.  The  maximum  CPU  quantum  is  one 
jiffy  (.01667  seconds),  disregarding  overhead 
of  the  operating  system.  When  overhead  rates 
of  10-25%  are  included,  the  CPU  quantums  have 
the  values  of  those  above.  Interactive  jobs 
usually  have  more  overhead,  since  they  deal 
more  with  system  utilities.  Batch  jobs 
usually  are  large,  number-crunching 
applications,  and  utilize  little  overhead. 

Disk  Access  Probabilities 

Since  no  other  information  was  available 
these  probabilities  were  arbitrarily  chosen 
to  be  equal  for  all  disks  and  for  all  job 


i 
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classes  to  simplify  the  analysis. 

Disk  Service  Times 

These  times  were  again  arbitrarily  chosen  to 
be  equal  for  all  disks  and  for  all  job 
classes.  The  value  of  the  service  time  is 
consistent  with  the  times  used  by  Chen. 

Long-term  Wait  Service  Time 

The  value  chosen  for  the  interactive  job 
class  represents  the  "think”  time  for  the 
user  (the  time  to  make  a  response  to  the 
computer)  as  well  as  the  other  long-term  wait 
service  times  common  to  both  interactive  and 
batch  jobs,  e.g.  command  wait  or  DAEMON  wait 
described  in  Chapter  2.  The  value  selected  is 
consistent  with  the  values  used  by  Chen. 

Since  jobs  in  the  batch  job  class  do  not  have 
a  "think"  time,  but  do  enter  other  long-term 
wait  states  common  to  both  interactive  and 
batch  jobs,  the  batch  job  class  was  assigned 
a  long-term  wait  service.  No  information  was 
available  to  determine  typical  values  for 
long-term  wait  service  times  for  batch  jobs. 
The  batch  long-term  wait  service  time  was 
chosen  to  be  significantly  less  than  the 
interactive  long-term  wait  service  time  since 
the  long-term  wait  states  possible  for  batch 
jobs  require  considerably  less  time  to 
satisfy,  and  are  therefore,  much  faster  than 
the  "think"  time  of  an  interactive  user. 

Job  Swap  Time  The  values  selected  are  consistent  with  the 
values  used  by  Chen.  The  swap  time  for  the 
batch  job  class  should  be  larger  since  there 
is  more  batch  job  memory  to  swap.  The  batch 
job  swapping  time  was  arbitrarily  doubled  to 
take  this  into  account. 

The  values  selected  for  the  CPU  time  per 
interaction  are  consistent  with  the  values 
used  by  Chen.  The  CPU  time  per  interaction 
for  interactive  jobs  is  lower  than  the  value 
for  batch  jobs  because  interactive  jobs  are 
generally  less  CPU-intensive  then  batch 
jobs. 

The  values  selected  for  the  number  of  disk 


CPU 


N 


DIO 
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I/O's  per  interaction  are  reasonable 
considering  the  values  used  by  Chen.  Since 
the  number  of  disk  I/O's  per  interaction  for 
interactive  jobs  is  lower  than  the  number  of 
disk  I/O's  per  interaction  for  batch  jobs, 
the  natural  conclusion  to  draw  is  that  the 
batch  jobs  are  more  I/O-intensive  than  the 
interactive  jobs  (usually  the  opposite  is 
true) .  A  more  insightful  approach  shows  that 
the  interactive  jobs  are  the  more 
I/O- intensive.  Equation  3-1  gives  the  formula 
for  the  number  of  CPU  quantums  received  per 
interaction,  Nr,c,„.  This  formula  and  Equation 
3-3  yield  a  \&lue  of  8  and  33  1/3  for 
interactive  and  batch  jobs,  respectively. 
Now  the  values  of  NDIQ  are  more  meaningful. 
In  the  case  of  the  iiixeractive  jobs,  half  of 
the  CPU  quantums  were  terminated  due  to  disk 
I/O  requests,  while  in  the  case  of  the  batch 
jobs,  only  15  percent  of  the  CPU  quantums 
were  terminated  due  to  disk  I/O  requests. 
Therefore  in  this  workload  configuration,  the 
interactive  jobs  are  the  more  I/O-intensive. 


5.2.3  DECsy stem-10  Swapping  Model 

The  multi-class  Chen  swapping  model  derived  from  the 
above  system  configuration  is  shown  in  Figure  5-1.  Note  that 
the  probabilities  and  service  times  for  the  network  are 
class  dependent  (r).  Node  1  corresponds  to  the  two  CPU's; 
nodes  2  through  7  represent  the  I/O  subsystem  containing  the 
6  disk  drives;  node  8  models  the  swapping  device;  and  node  9 
represents  the  long-term  wait  service  node. 

The  following  are  definitions  of  the  important 
probabilities  in  the  network 

P91r  Pgir  is  the  probability  that  the  job  in  class 
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Figure  5-1  DECsystem-10  Swapping  Model  (For  Class  r) 
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p98r 


P 


18r 


r  leaving  a  long-term  wait  state  (the  "think" 
mode  for  interactive  jobs)  is  already  in 
memory  and  can  proceed  directly  back  to  the 
CPU  (analogous  to  P.,  in  Chen's  original 
model)  .  41 


Pqgr  is  the  probability  that  a  job  in  class  r 
leaving  a  long-term  wait  state  (the  "think" 
mode  for  interactive  jobs)  must  be  swapped 
into  memory  (analogous  to  P.,  in  Chen's 
model) .  43 


P,g  is  the  probability  that  a  job  in  class  r 
finishing  its  CPU  quantum,  must  be  swapped 
out  to  make  room  for  another  job  which  is 
ready  to  execute,  but  has  been  waiting  for 
memory  (analogous  to  P^  in  Chen's  model). 

P..  ,  for  j  *  2, 3... 7,  are  the  probabilities 
or-3^  job  in  class  r  going  to  the  various 
disks  in  the  I/O  subsystem  from  the  CPU 
(analogous  to  P12  in  Chen's  model). 


5.«.3L  MfiAel-C<?mpaiiS9ns 


5.3.1  Definition  of  Performance  Measures 


The  following  are  definitions  of  the  performance 
measures  considered  in  these  comparisons  (Ref  29:10-12) 

CPU  utilization 

The  percent  of  time  the  CPU  is  processing  a 
job. 

Swapper  utilization 

The  percent  of  time  the  swapping  device  is 
swapping  in  or  out  a  job. 
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Response  Time  Elapsed  time  between  entering  the  last 
character  of  a  request  at  a  terminal  and 
receiving  the  first  character  of  the 

response. 


5.3.2  Tabular  Results 

The  results  of  the  models  are  shown  in  Tables  5-2 
through  5-4.  Table  5-2  contains  the  performance  measure 
results  for  the  classical  swapping  model  with  the  three 
workload  cases.  Table  5-3  contains  the  results  for  the  Chen 
swapping  model  with  the  three  workload  cases.  Table  5-4 
shows  the  behavior  of  the  probability  transition  matrix  for 
different  main  memory 

Table  5-2  Classical  Swapping  Model  Results 


Performance  Measures 

Workload 

CPU 

Swapper 

Response 

Cases 

Utilization 

Utilization 

Time 

(«) 

(«) 

(sec) 

Two  Class 

100.00 

33.46 

6.750 

Average 

99.98 

33.23 

8.756 

No  Batch 

81.05 

26.97 

3.122 

Important  Prob ability  Matrix 

Values 

Pllr 

P  P 

*18r  *91 r 

Two  Class 

.466/. 933 

0  0 

Average 

.699 

0  0 

No  Batch 

.466 

0  0 
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Table  5-3  Chen  Swapping  Model  Performance  Results 


Memory 

Workload 

CPU 

Utilization 

Swapper 

Utilization 

Response 

Time 

(K) 

(%) 

(%) 

(sec) 

Two  Class 

98.78 

6.44 

.828 

240 

Average 

70.95 

4.00 

1.083 

No  Batch 

— 

— 

Two  Class 

98.58 

25.72 

.847 

200 

Average 

70.80 

16.00 

1.107 

No  Batch 

— 

— 

— 

Two  Class 

97.94 

51.14 

.895 

150 

Average 

70.51 

32.00 

1.154 

No  Batch 

— 

— 

- 

Two  Class 

95.84 

75.24 

1.022 

100 

Average 

70.05 

48.00 

1.227 

No  Batch 

30.11 

12.50 

.504 

Two  Class 

89.15 

93.06 

1.400 

50 

Average 

69.24 

64.01 

1.350 

No  Batch 

29.95 

25.00 

.557 

Two  Class 

87.07 

95.10 

1.521 

40 

Average 

68.95 

67.89 

1.407 

No  Batch 

29.91 

27.50 

.570 

Two  Class 

85.12 

96.48 

1.640 

30 

Average 

68.65 

71.28 

1.458 

No  Batch 

29.87 

30.00 

.583 

Two  Class 

84.21 

96.99 

1.697 

20 

Average 

68.49 

72.81 

1.484 

No  Batch 

29.83 

32.50 

.598 

Two  Class 

82.52 

97.76 

1.804 

10 

Average 

68.18 

75.54 

1.536 

No  Batch 

29.78 

34.99 

.613 
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Table  5-4  Important  Probability  Matrix  Values 


Memory 

(K) 

Workload 

Pllr 

P18r 

P91r 

Two  Class 

.374/. 820 

0 

.934 

240 

Average 

.521 

0 

.941 

No  Batch 

— 

— 

— 

Two  Class 

.374/. 820 

0 

.735 

200 

Average 

.521 

0 

.764 

No  Batch 

— 

— 

Two  Class 

.374/. 820 

0 

.468 

150 

Average 

.521 

0 

.526 

No  Batch 

— 

— 

Two  Class 

.374/. 820 

0 

.188 

100 

Average 

.521 

0 

.282 

No  Batch 

.374 

0 

.641 

Two  Class 

.359/. 816 

.015/. 003 

0 

50 

Average 

.521 

0 

.025 

No  Batch 

.374 

0 

.286 

Two  Class 

.351/. 814 

.023/. 005 

0 

40 

Average 

.519 

.002 

0 

No  Batch 

.374 

0 

.214 

Two  Class 

.344/. 812 

.030/. 007 

0 

30 

Average 

.515 

.006 

0 

No  Batch 

.374 

0 

.141 

Two  Class 

.341/. 811 

.033/. 008 

0 

20 

Average 

.513 

.008 

0 

No  Batch 

.374 

0 

.068 

Two  Class 

.335/. 810 

.039/. 009 

0 

10 

Average 

.510 

.011 

0 

No  Batch 

.374 

.000 

0 
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sizes  and  workload  cases.  The  Chen  versus  classical 
swapping  model  is  discussed  first  to  verify  the  Chen 
swapping  algorithm  results.  Then  a  comparison  between  the 
two  class  model  versus  the  single-class  models  is 
presented. 

5.3.3  Chen's  Versus  the  Classical  Swapping  Model 

The  following  paragraphs  consider  the  impact  of  using 
Chen's  swapping  model  versus  using  the  classical  swapping 
model.  The  performance  measures  will  be  utilization  of  the 
CPU  and  the  swapping  device,  and  the  response  time  received 
by  the  interactive  users.  The  results  shown  are  only  for 
the  two  class  model.  Similar  results  were  achieved  with  the 
average  class  model  and  the  no  batch  model. 

5. 3. 3.1  Probability  Structure 

The  primary  difference  between  the  classical  swapping 
model  and  Chen's  swapping  model  is  that  Chen's  probability 
transition  matrix  is  dynamic,  i.e.  the  probabilities  in  the 
matrix  are  dependent  on  the  main  memory  size,  the  number  of 
jobs  competing  for  memory,  and  the  job  sizes. 

The  Figure  5-2  is  a  plot  of  the  data  in  Table  5-4  for 
the  two  class  workload  model.  It  illustrates  the  dynamic 
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Figure  5-Z  Probability  Structure  (Two 


Analytical  Modeling  Results 


relationship  between  the  state-dependent  probabilities  and 
main  memory  size  produced  by  the  two  class  workload  model. 
In  the  classical  model  P^Sr  an<*  P91r  wou^  be  zero  and 
would  be  constant.  With  Chen's  model,  the  probabilities 
have  some  interesting  properties 

-  pglr  behaves  "classically"  (Pg.  =0)  from  zero  to 
approximately  70K  of  main  memory.  At  that  point  the 
Chen  approximation  takes  effect  and  PQ1  steadily 

increases  to  a  probability  of  one  at  23U TC  of  main 
memory. 


-  Pllr  and  pi8n  behavior  is  reversed.  For  ranges  of  main 
memory  between  70K  and  250K,  the  probabilities  behave 
"classically"  (P. ,  -constant ,  P,a  “0).  Between  10K  and 
70S 


approximately  76ft  of  main 
approximation  takes  effect.  In 
inversely  proportional  to  P. 

P18r's  loss* 


memory,  the 
this  range,  P 


1 .  e . 


llr 


s  g 


Chen 
is 
is 


Similar  phenomena  occur  in  the  average  class  model  and  the 
no  batch  model. 


The  rationale  for  the  behavior  of  Pglr  is  as  follows. 
For  small  amounts  of  main  memory,  all  jobs  in  long-term  wait 
will  be  swapped  out  to  make  room  for  other  jobs  waiting  for 
memory.  Hence,  whenever  a  job  gets  blocked  at  the  CPU  and 
goes  to  a  long-term  wait  state,  it  will  be  swapped  out  and 
must  therefore  always  be  swapped  back  in  when  it  finishes 
its  long-term  wait  service.  For  large  amounts  of  main 
memory,  there  is  enough  memory  to  satisfy  all  memory  demands 
as  well  as  keep  all  jobs  in  in  a  long-term  wait  state  in 
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memory.  As  main  memory  size  approaches  250K,  an  increasing 
number  of  jobs  in  a  long-term  wait  state  may  stay  in 
memory.  Therefore,  this  increases  the  probability  that  a 
job  leaving  a  long-term  wait  state  will  already  be  in  main 
memory  and  need  not  be  swapped  in. 

The  rationale  for  the  behavior  of  Pllr  and  P18r  is  as 
follows.  For  large  amounts  of  main  memory,  there  is  enough 
for  all  executable  jobs  to  remain  in  memory.  As  main  memory 
size  gets  smaller,  a  larger  number  of  jobs  in  long-term  wait 
are  swapped  out.  This  continues  until  all  jobs  in  the 
long-term  wait  states  swapped  out.  After  all  the  inactive 
jobs  have  been  swapped  out,  the  Job  Swapper  must  start 
swapping  out  jobs  that  are  executable  in  order  to  make  room 
for  other  jobs  in  the  system.  Therefore,  jobs  expiring 
their  CPU  quantum  now  have  higher  probabilities  of  being 
swapped  out,  instead  of  returning  to  the  end  of  the  CPU 
queue  for  more  service  (P18r  becomes  larger) . 

5. 3. 3. 2  Performance  Predictions 

The  figures  5-3  a)  through  5-3  c)  are  plots  of  the  data 
contained  in  tables  5-2  and  5-3  for  the  two  class  model. 
The  horizontal  lines  on  each  plot  represent  the  values 
predicted  by  the  classical  model,  while  the  data  points 
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represent  the  values  predicted  by  the  Chen  model  at  each 
memory  size.  Note  the  sensitivity  of  the  performance 
predictions  due  to  the  effects  of  main  memory  size  in  Chen's 
swapping  model.  Specifically,  the  following  effects  are 
modeled  in  the  Chen  swapping  model 

-  In  the  range  from  zero  main  memory  to  the  main  memory 
size  at  the  intersection  of  the  Chen  swapping  model 
curve  and  the  classical  swapping  model  curve 
(approximately  70K  for  the  two  class  model) ,  the  Chen 
model  includes  the  effect  of  swapping  out  executable 
jobs  on  the  performance  measures  (Pig) •  Compared  to  the 
classical  swapping  model,  the  Chen  swapping  model  shows 

*  Increased  response  time  and  swapping  utilization 
caused  by  increased  swapping  activity. 

*  Decreased  CPU  utilization  caused  by  less  jobs  in 
memory  competing  for  service  from  the  CPU's. 

-  For  the  range  of  memory  values  between  the  point  of 
intersection  and  enough  memory  to  contain  all  jobs 
(250K  for  the  two  class  model) ,  the  Chen  model  includes 
the  effect  of  gradually  having  to  swap  in  and  out  fewer 
jobs  in  long-term  wait  as  opposed  to  the  classical 
assumption  that  jobs  in  long-term  wait  are  always 
swapped  in.  Compared  to  the  classical  swapping  model, 
the  Chen  swapping  model  shows 

*  Decreased  response  time  and  swapping  utilization 
caused  by  decreased  swapping  of  the  jobs  in  a 
long-term  wait  state. 

*  Increased  CPU  utilization  caused  by  more  jobs  in 
memory  competing  for  service  from  the  CPU's. 


5.3.4  Chen's  Swapping  Model:  Multi-class  Versus  Single-class 

The  following  paragraphs  compare  the  results  of  the  two 
class  model  to  the  average  class  and  no  batch  class  workload 
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models  when  using  Chen's  swapping  model.  The  pecformance 
measures  considered  will  be  the  same  as  those  used  in 
comparing  Chen's  swapping  model  versus  the  classical 
swapping  model:  utilization  of  the  CPU  and  the  swapping 
device,  and  the  response  time  received  by  the  interactive 
users. 


The  figures  5-4  a)  through  5-4  c)  are  plots  of  the  data 
contained  in  Table  5-3  for  the  two  class,  average  class,  and 
no  batch  workload  models.  These  plots  reveal  the  wide 
disparity  between  the  results  obtained  from  a  single  class 
approximation  when  actually  two  classes  of  jobs  exist.  In 
particular,  one  can  note  the  following  results 

-  The  no  batch  workload  model  seriously  underestimates 
the  response  time  of  interactive  jobs.  For  the  average 
workload  model,  the  response  time  is  overestimated  for 
main  memory  sizes  which  allow  all  active  jobs  to  reside 
in  memory.  For  main  memory  sizes  below  this  value, 
response  times  are  underestimated. 

-  The  single-class  workload  models  underestimate  the  CPU 
and  swapper  utilizations.  This  is  especially  true  for 
the  no  batch  workload  model.  As  main  memory  sizes 
increases,  the  average  workload  model  approaches  the 
two  class  model. 

-  For  all  performance  measures  considered,  the  average 
workload  model  was  closer  to  the  two  class  workload 
model  than  the  no  batch  workload  model. 
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5.4  Chapter  Summary 


This  chapter  presented  the  results  of  modeling  a 
computer  system  using  the  Chen  model  extended  to  multi-class 
models  for  three  different  workload  cases:  two  class, 
average  class,  and  no  batch.  The  hardware  and  workload 
configuration  parameters  were  chosen  to  as  closely  model  the 
Avionics  Lab  system  configuration  as  possible  and  to  remain 
consistent  with  the  workloads  considered  by  Chen. 


The  results  of  the  comparison  between  the  Chen  swapping 
model  and  the  classical  swapping  model  are  summarized  below 

-  For  small  amounts  of  main  memory 

*  Longer  response  time  and  higher  swapping 

utilization  caused  by  increased  swapping 

activity. 

*  Lower  CPU  utilization  caused  by  more  jobs  waiting 

for  swapping  service  and  less  jobs  in  memory 

competing  for  CPU  srvice. 

-  For  larger  amounts  of  main  memory 

*  Shorter  response  time  and  lower  swapping 

utilization  caused  by  decreased  swapping  of  the 
jobs  in  a  long-term  wait  state. 

*  Higher  CPU  utilization  caused  by  less  jobs  waiting 

for  swapping  service  and  more  jobs  in  memory 

competing  for  CPU  service. 


The  following  results  were  obtained  from  the 

comparisons  of  the  two  class  workload  model  versus  the 
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single  class  workload  models  (average  class  and  no  batch 
class) 

-  The  no  batch  workload  model  seriously  underestimates 
the  response  time  of  interactive  jobs  for  all  memory 
sizes,  while  the  average  workload  model  overestimates 
response  times  for  large  main  memory  sizes  and 
underestimates  response  times  fo-  small  main  memory 
sizes. 

-  The  single-class  workload  models  underestimate  the  CPU 
and  swapper  utilizations. 

-  For  all  performance  measures  considered,  the  average 
workload  model  was  closer  to  the  two  class  workload 
model  than  the  no  batch  workload  model. 
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6.1  Conclusions 

From  the  results  obtained  in  Chapter  5  for  the  workload 
described,  the  following  conclusions  can  be  made 

1.  Chen's  swapping  model  is  a  more  realistic  model  of  a 
computer  system  than  the  classical  swapping  model  due 
to  the  fact  that  Chen's  model  more  accurately  models 
the  effects  of  the  scarcity  of  main  memory  on 
job-swapping  behavior. 

2.  Adding  the  batch  class  to  Chen's  swapping  model 
improves  the  accuracy  of  the  performance  measure 
predictions  when  the  system  contains  significant  batch 
activity. 
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$.2  Pecominendations  for  Future  Research 


The  following  are  topics  of  studies  which  could  be 
conducted  as  an  extension  of  this  thesis 

1.  The  current  model  of  the  DECsystem-10  shown  in  Figure 
5-1  is  a  highly  abstract  model  that  does  not  consider 
many  of  the  complexities  of  the  TOPS-10  monitor  or  the 
hardware  configuration.  As  an  extension  to  the  model , 
the  following  future  thesis  efforts  could  be  conducted 

-  Modeling  the  CPU's  in  the  DECsystem-10  as  a 
single  processor-sharing  node  may  oversimplify 
the  complex  nature  of  the  HPQ's,  PQ1,  and  PQ2 
processor  queues.  Including  multiple  queues  with 
job  class  priorities  may  provide  a  more  accurate 
model  of  the  overall  system  behavior.  These 
multiple  queues  and  job  class  priorities  may  be 
implemented  using  the  shortest  elapsed  time  or 
shortest  remaining  processing  time  queueing 
disciplines  (Ref  6:178-186). 

-  Modeling  the  I/O  subsystem  by  just  modeling  the 
disk  drives  may  be  overlooking  long  delays  due  to 
contention  for  the  I/O  channel.  An  approximation 
method  described  in  reference  16  called  the 
method  of  surrogate  delays  appears  to  be  a  prime 
candidate  for  inclusion  to  the  present  computer 
program  to  more  accurately  model  I/O  channel 
contention. 

-  An  important  extension  to  the  present  model  would 
be  to  determine  how  the  numerous  scheduling 
parameters  within  the  DECsystem-10  affect  the 
workload  configuration  parameters,  service  times, 
and  the  probability  transition  matrices  in  the 
model.  This  information  would  provide  guidance 
to  systems  analysts  in  setting  optimum  values  for 
the  scheduling  parameters. 

2.  Although  the  current  model  does  provide  reasonable 
results,  the  model  has  not  been  validated.  To 
accomplish  this  validation,  one  of  the  following  two 
approaches  could  be  taken 
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-  Take  measurements  on  the  actual  system,  and 

compare  the  measurement  results  to  the  results 
predicted  by  the  model. 

-  If  measurements  on  the  actual  system  are 

impossible,  the  development  of  an  extremely 
detailed  simulation  model  which  takes  into 

account  most  of  the  intricacies  of  the  TOPS-IO 
monitor  would  be  a  step  toward  validation.  This 
would  not  completely  validate  the  analytical 
model  since  the  simulation  itself  is  not 
validated,  but  the  modeling  accuracy  possible 
with  a  simulation  should  provide  a  more  accurate, 
but  time  consuming  model. 

3.  The  impact  of  changes  in  the  workload  configuration 
using  the  current  model  would  provide  better  insights 
as  to  how  the  performance  measurements  are  affected  by 
the  workload  configuration.  Some  of  the  most 
interesting  variations  might  be  changes  in  the 

-  Job  Size 

-  Number  of  batch  jobs  in  the  system 

-  Number  of  competing  jobs  in  the  system,  i.e.  not 
in  a  long-term  wait  state 

-  Number  of  CPU  versus  I/O  intensive  jobs. 

4.  Another  type  of  workload  variation  entails  increasing 
the  number  of  job  classes  in  the  system.  To  model  a 
large  number  of  job  classes  (greater  than  4)  in  the 
system  will  require  a  restructuring  of  the  current 
program.  The  program  currently  computes  the  necessary 
probability  distributions  one  node  at  a  time,  since 
the  state  space  for  the  whole  network  would  be  too 
large  to  run  on  the  CDC  computer.  With  a  larger 
number  of  job  classes,  the  state  space  of  a  node 
becomes  unmanageable.  Therefore,  the  computations 
will  have  to  be  broken  down  even  further,  in  order  to 
fit  in  memory.  If  the  program  rewrite  is  undertaken, 
one  of  the  two  alternatives  listed  below  should  be 
chosen 

-  Rewrite  the  program  in  FORTRAN77  to  take 
advantage  of  the  language's  increased  structure 
and  more  importantly,  to  incorporate  the  random 
record  reads  and  writes  that  are  now  part  of  the 
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language.  This  will  eliminate  the  installation 
dependent  implementation  of  random  record 
processing  and  make  the  program  more  portable. 

Rewrite  the  program  in  PASCAL,  PL/I  or  some  other 
language  that  allows  recursion.  Since  the 
computational  algorithms  are  defined  recursively, 
this  may  provide  more  efficient  code  and  will 
definitely  be  a  more  straight  forward 
implementation. 
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PROGRAM  P€RF0RH(INPUT<0UTPUT,TAPE5*INPUT»TAPE6*0UTPUTiTAPE3) 


C  THE  PROGRAM  PERFORM  tS  A  CLOSED  QUEUEING  tCTWGRK 

C  MODEL  WHICH  CM  BE  USED  IN  PERFORMANCE  EVALUATION  OF 
C  TDE-SHWINC  COTUT®  SYSTEMS.  THE  PROGRAM  PERFORM  IS 
C  BASED  ON  A  aOSED  QUEUEING  NETWORK  MODEL  DEVELOPED  BV 
C  F.  BASKETT.K.  CHANDY,  R.  MJNTZ,  MD  F.  PALACIOS. 

C  TTCIR  MODEL  HAS  PRESENTED  IN  Tt€  ARTICLE  "OPEN,  CLOSED, 

C  MD  NIXED  ICTWCRK  OF  QUEUES  WITH  DIFFERENT  CLASSES  OF 

C  CUSTOMERS* >  WHICH  APPEARED  IN  THE  JQJRNAL  OF  THE  ACM, 

C  VOL.  22,  MO.  2.  PP.  248-260,  APRIL,  1975. 

C 

C  THE  COMPUTATIONAL  ALGORITHMS  WHICH  ARE  USED  TO 

C  IMPLEMENT  THE  MODEL  OF  BASKETT,  ET.  AL.  ARE  BASED  OH  TO 
C  COPUTATIONALLY  EFFICIENT  ALGORITHMS  PRESENTED  IN  THE 
C  DISSERATION  OF  J.  WONG  ENTITLED  QUEUEING  NETWORK  MODELS 
C  FOR  COFUTER  SYSTEMS,  UCLA-ENG-7579,  UNIVERSITY  OF 
C  CALIFORNIA  AT  LOS  ANGELES,  OCTOBER.  1975. 

DIMENSION  E(13,2),  EKPVAL(13,2>,  DEPART(13,2),  FNT(IOOO), 
llFACTOR(Al),  INDEX3(40),  ISKIP<2),  ISTATEK2),  ISTATE212), 
2MMETYPU3),  NORNO]N(lOOO,2),  NUSERS<2>.  PR0B(13,I3), 
SPRODCUHdOOO),  PROBHARdOOOI.  PR0BPAR(13,2,41>, 
4PROBTOTU3.41I,  SERVICE!  13,21.  070.(13,2).  UAREA<13)> 
51(13),  IDEPI13),  OEPI 13,2,40) i  TCPUI2),  010(2), 
6DISKPRB(13,2).SIZEJ0B(2) 

DOUCLE  PRECISION  IF  ACTOR 
LOGICAL  PRINTON 
INTEGER  IUSERS.  NAVGJOB 

REAL  DELTA,  TOLERAN,  PEXPVAL.  SIZEtEH,  AVGJOBS 
MAX  NODE  -  13 
MAXUSER  *  40 
MAXUSE1  *  MAXUSER  +  1 
HAXTYPE  *  2 
MAXSTAT  «  1000 
ISIZE3  *  40 

HAXNOOE  —  THE  HAX1HUM  NUMBER  OF  MOOES  WHICH  CAN  BE 
REPRESENTED  IN  THE  PROERAM. 

MAXUSER  —  THE  MAXIMUM  NUMBER  OF  TOTAL  USERS  THAT  CAN 
BE  REPRESENTED  IN  THE  PROGRAM. 


MAXUSER!  —  THE  MAXIMUM  NUMBER  OF  TOTAL  USERS  THAT  CAN 
BE  REPRESENTED  IN  THE  PROGRAM  ♦  i. 

HAXTYPE  —  THE  MAXIMUM  XIMBER  OF  TYPES  OF  USERS  THAT 
AN  E  REPRESENTED  IN  THE  PROGRAM. 

MAXSTAT  —  THE  MAXIMUM  STATE  SPACE  FOR  EACH  NOE. 


PERFORM  2 
PERFORM  3 
PERFORM  4 
PERFORM  5 
PERFORM  6 
PERFORM  7 
PERFORM  8 
PERFORM  9 
PERFORM  10 
PERFORM  11 
PERFORM  12 
PERFORM  13 
PERFORM  14 
PERFORM  15 
PERFORM  16 
PERFORM  17 
PERFORM  18 
PERFORM  19 
COXMOO  1 
COXMOD  2 
COXHOD  3 
COXMOD  4 
COXMOD  5 
COXMOD  6 
COXMOO  7 
COXHOD  8 
COXMOD  9 
COXMOD  10 
COXHOD  11 
COXMOD  12 
COXMOD  13 
COXMOO  14 
COXHOD  15 
COXMOD  16 
COXMOO  17 
PERFORM  33 
PERFORM  34 
PERFORM  35 
PERFORM  36 
PERFORM  37 
PERFORM  38 
PERFORM  39 
PERFORM  40 
PERFORM  41 
PERFORM  42 
PERFORM  43 
PERFORM  44 
PERFORM  45 
PERFORM  46 
PERFORM  47 
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c  isize3  —  he  maximum  mem  of  random  records  <30) 

C  WHICH  CAN  BE  WRITTEN  ON  TAPE3  ♦  l.  ISIZE3  IS 

C  CALCULATED  BY  HAXNODE  •  3  ♦  i. 

E(I* J)  --  THE  MEAN  ARRIVAL  RATE  OF  TYPE  J  USERS  AT 
SERVICE  CENTER  I.  THE  DIMENSIONS  OF  E  ARE 
(HAXNODE. HAXTYPE). 

EXPVALll.J)  —  TT€  EXPECTED  NUMBER  OF  TYPE  J  USERS 

AT  NODE  I  AT  EQUILIBRIUM.  TT€  DIMENSIONS 
OF  EXPVAL  ARE  I HAXNODE.  HAXTYPE). 

DEPART(I.J)  —  THE  DffARTDRE  RATE  OF  TYPE  J  USERS  FROM 
NODE  I  AT  EQUILIBRIUM.  THE  DIMENSIONS 
OF  DEPART  ARE  (HAXNODE,  HAXTYPE). 

FNT(I)  —  THE  FUNCTION  REPRESENTING  THE  RBJTTIVE 

WEIGHTING  OF  STATE  I  IN  THE  CALCULATIONS  OF 
THE  EQUILIBRIUM  STATE  PROBABILITY  DISTRIBUTIONS. 

THE  DIMENSION  OF  FNT  IS  (HAXSTAT). 

IFACTOR(I)  —  VALUE  OF  I  FATORIAL.  DIMENSION  OF 
IFACTOR  IS  (HAXUSER1). 

IMJEX3CI )  —  THE  LOCATION  OF  RECORD  I  ON  THE  RANDOH 
FILE  TAPE3.  THE  DIMENSION  OF  INOEI3  IS 
(ISIZE3). 

ISKIP  —  WORKING  ARRAY  USED  TO  SKIP  OVER  (NEEDED 
STATES  IN  THE  STATE  SPACE  FOR  A  GIVEN  NODE. 

TIE  DI TENSION  OF  ISKIP  IS  (HAXTYPE). 

ISTATE1  —  WORKING  ARRAY  USED  TO  SPECIFY  A  STATE  OF  A 

GIVEN  NODE.  BE  D  WENS  ION  OF  ISTATE1  IS  (HAXTYPE). 

ISTATE2  —  WORKING  ARRAY  USED  TO  SPECIFY  A  STATE  OF  A 
GIVEN  NODE.  THE  DIIENSION  OF  ISTATE2  IS 
(HAXTYPE). 

NOOETYPII)  —  SPECIFIES  THE  TYPE  OF  NODE  I .  THE 
DIMENSION  OF  NOOETYP  IS  (HAXNODE). 

NORTCON  —  WORKING  ARRAY  USED  TO  CALCULATE  THE  NORMAL¬ 
IZATION  CONSTANT.  THE  DIMENSIONS  OF  NQRHCON 
ARE  (HAXSTAT. 2). 

NUSERS(I)  —  NUMBER  OF  TYPE  I  USERS  REPRESENTED  IN  THE 

HOOEL.  THE  DIMENSION  OF  HUSERS  IS  (HAXTYPE). 

PROB(I.J)  —  THE  PROBABILITY  THAT  A  TASK  OF  A  GIVEN  TYPE 
WIU.  PROCEED  TO  NODE  IWOi  IT  DEPARTS 
WGte  J.  THE  D PENSIONS  OF  PROS  ARE 
(NMNK.HAXMOOE). 


PERFORM  48 
PERFORM  49 
PERFORM  50 
PERFORM  51 
PERFORM  52 
PERFORM  53 
PERFORM  54 
PERFORM  55 
PERFORM  5 6 
PERFORM  57 
PERFORM  58 
PERFORM  59 
PERFORM  SO 
PERFORM  61 
PERFORM  62 
PERFORM  63 
PERFORM  64 
PERFORM  65 
PERFORM  66 
PERFORM  67 
PERFORM  68 
PERFORM  69 
PERFORM  70 
PERFORM  71 
PERFORM  72 
PERFORM  73 
PERFORM  74 
PERFORM  75 
PERFORM  76 
PERFORM  77 
PERFORM  78 
PERFORM  79 
PERFORM  80 
PERFORM  81 
PERFORM  82 
PERFORM  S3 
PERFORM  84 
PERFORM  85 
PERFORM  86 
PERFORM  87 
PERFORM  88 
PERFORM  89 
PERFORM  90 
PERFORM  91 
PERFORM  92 
PERFORM  93 
PERFORM  94 
PERFORM  95 
PERFORM  96 
PERFORM  97 
PERFORM  98 
PERFORM  99 
PERFORHIOO 
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c 

C  PRODCUMdl  —  THE  PROBABILITY  THAT  A  SUBNETWORK  OF 
C  NOCES  IN  THE  NETWORK  10DEL  IS  IN  A  GIVEN 

C  STATE.  THE  STATES  IN  THE  STATE  SPACE  ARE 

INDEXED  BY  TT£  SUBSCRIPT  I.  TIE  DIMENSION 
OF  PROBCUN  IS  (HAISTAD. 

FROBHAR(I)  —  TIE  PROBABILITY  THAT  A  SINGLE  NODE  IN 
THE  NETWORK  KIEL  IS  IN  A  GIVEN  STATE. 

TIE  STATES  IN  TIE  STATE  SPACE  ARE  INDEXED 
BY  THE  SUBSRIPT  I.  TIE  DIMENSION  OF 
PROBMAR  IS  (HAXSTAT). 

PROfPAR(IiJ.K)  —  THE  PROBABILITY  THAT  THERE  ARE  K 
USSTS  OF  TYFE  J  AT  NODE.  THE 
DIIEK3IONB  OF  PR3BPAR  ARE  (HAXNODE. 
MAXTYPE.MAXUSE1). 

PROBTOTtliJI  —  TIE  PROBABILITY  THAT  TIERE  ARE  J  USERS 
AT  NODE  (.  THE  DIMENSIONS  OF  PROBTOT 
ARE  (HAXNODE, MAXTYPE). 

SERVICE(I,J>  —  TIE  SERVICE  TIME  OF  A  TYPE  J  USER  AT 
NODE  I.  THE  DIlENStmS  OF  SERVICE 
ARE  (HAXNODE, MAXTYPE). 

UTILdiJ)  —  TIE  UTILIZATION  OF  NODE  I  BY  USERS  OF  TYPE 
J.  TIE  DIMENSIONS  OF  UTIL  ARE  (HAXNODE* 
MAKTYPE). 

UAREA  —  WORKING  ARRAY  USED  TO  CALCULATE  TIE  MEAN  ARRIVAL 
RATE  OF  USERS  AT  THE  NODES  IN  THE  IETWORK. 

THE  DIHENSION  OF  WAREA  IS  (HAXNODE). 

X  —  NIRKINO  ARRAY  USED  TO  CALCULATE  TIE  MEAN  ARRIVAL 
RATE  OF  USERS  AT  THE  NODES  IN  THE  NETWORK. 

DIHENSION  OF  X  IS  (HAXNODE). 

IDEP(I)  —  INDICATES  STATE  DEPENDENT  SERVICE 
RATES  AT  NOOE  I. 

1  —  NO  STATE  DEPENDENT  SERVICE  RATES 

2  —  SERVICE  RATES  DEPENDENT  ON  THE  TOTAL 
NUMBER  OF  USERS  AT  THE  NODE. 

3  —  SERVICE  RATES  DEPENDENT  ON  THE  NUHBER 
OF  USERS  OF  A  TYPE  AT  THE  NODE. 

THE  DIMENSION  OP  IDEP  IS  (HAXNODE). 

C  DEP(I,J,KI  —  THE  RELATIVE  SERVICE  RATE  OF  TYPE  J  USERS 
C  AT  NODE  I  WHEN  TIERE  ARE  K  USERS  PRESENT 

C  AT  NODE  I.  THE  DIIENSIONS  OF  DEP  ARE 

C  (NAUWDE,  MAXTYPE,  HAXU3ERI. 

NRITE(M) 

1  FORMAT! IHl.PX.’THIS  IS  A  CLOSED  QUEUEING  NETWORK  MODEL  *, 
1*  FOR  THE-SHARINO  COMPUTER  SYSTEMS.*) 


PERFQRNIOI 
COXHOD  18 
CQXHOD  19 
COXHOD  20 
COXHOD  21 
COXHOD  22 
PERF0RH108 
PERF0RM109 
PERFORMliO 
PERFORMlll 
PERF0RMU2 
PERF0RM113 
PERFORM U 4 
PERF0RM115 
PERF0RM11& 
PERF0RM117 
PERF0RM118 
PERF0RI1119 
PERF0B1120 
PERF0RM121 
PERF0RM122 
PERF0RM123 
PERF0RM124 
PERF0RM125 
PERFQRH12L 
PERF0RNI27 
PERF0RM128 
PERF0RM129 
PERF0RM130 
PERF0RM131 
PERF0RM132 
PERF0RM133 
PERF0RN134 
PERF0RN13S 
PERF0RM136 
PERF0RM137 
PERFORM  38 
PERFORM l 39 
PERF0RM140 
PERF0RM141 
PERF0RM142 
PERF0RM143 
PERF0RM144 
PERF0RM145 
PERFORMS 
PERF0RH147 
PERF0RM148 
PERF0RM149 
PERF0RM15O 
PERF0RM151 
PERFORM! 52 
PERF0RM153 
PERF0RM154 
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C  OPEN  THE  MUSS  STORAGE  FILE 

C 

CM1  0PENNSI3,I>UEX3.ISIZE3,0I 
C 

C  CALL  INPUT  TO  INPUT  THE  DATA  NECESSARY  TO  EXECUTE  THE 
C  PROGRAM 

C 

CALL  INDATA(Pfl38,  SERVICE. E,NODETYP,NUSERS, HAXNODE, 
mAXUSER,l«XTYPE,HAISTAT>«I)ES.  IUSERS,  ITYPES, X.INDEX3, 1SIZE3. 
2MAREA,  ISTATES,  IDEP,OEP,EXim,NDlSKS,NAXITER,TOLERAN, 
3DISKPRB, TCPU> DIO> SI2EJOB, AVGJ08S, NAVGJOB) . RETURNS! 10 i 
PRINTON  =  .TRUE. 

PEXPWAL  =  0.0 
NITER  =  0 
DO  4  J«1,ITYFE$ 

00  4  t=L. NODES 
IF  (J  .HE.  1)  GO  TO  2 
EXPVAL(I,J)  *  IUSERS/NQDES 
GO  TO  4 

2  EXPVAL(1,.J)  *  0.0 

4  CONTINUE 
C 

C  CALL  CHEN75  TO  COMPUTE  PROBABILITY  TRANSITION  MATRIX  PROS 
C 

5  CALL  CHEM75IPR0B,  NDIBKS, D1SKPRB,  TCPUi  QIQ.NAWGJQB, 

1EXPVAL, SERVICE,  DELTA,  fEXPVAL, IUSERS. HAXNODE,  MAXTYPE, 

2ITYPES, NODES,  NITER.  E,  X,  IftfiEA) 

C 

C  CALL  FUNCT  TO  CALCULATE  THE  FUNCTION  FNT  FOR  ALL  NODES. 

C 

CALL  FUNCTIE, SERVICE,  IFACTOR,  ISTATEL , FNT ,  MUSERS . NOOETYP, 
INODES, IUSERS,  I  TYPES. MAXNODE.MAXUSER, MAXUSE1 .MAXTYPE, 
2MMSTAT,  INDEX3,  ISIZE3,  ISTATES,  IDEP,  DEP)  .RETUWISI 10) 

C 

C  CALL  NORMAL  TO  CALCULATE  THE  NORMALIZATION  CONSTANT. 

C 

CALL  NGRMAL(R1T,N0RNC0N,NU^RS,ISTATE1,ISTATE2,ISK1P, 
INAXNOOE.MAXSTAT,  NODES,  ITYPES,  MAXTYPE,  INLETS,  ISIZE3,  ISTATES) 


CALL  MWGIN  TO  CALCULATE  (Hi  MARGINAL  PROBABILITY 
DISTRIBUTION  FOR  T*  STATE  OF  THE  NODES. 


CAUL  MAROCN(FNT,NORMCCN,PROSHAR,PflOBCUM,  (USERS,  ISTATE1,  ISTATE2, 

1 ISKIP, HAXNODE, HAXSTA7 , MAXTYPE, NCCES , ITYPES, INDEX3, ISIIE3, ISTATES) 

CAU.  EXPECT  TO  CA1XU.ATE  VARIOUS  PERFORMANCE  MEASURES. 

CALL  EXPSCr(PftDBMW>Pfi)3BT0TiPR0BPAR,NUSERS,  ISTATE1, 
lNOOETYP,EXPVAL,DEPART,uriL, SERVICED, HAXNODE, 

2MAXUSE1 ,  MAXTYPE, NODES,  [TYPES,  IUSERS,  NAXSTAT,  INDEX3,  ISIZE3,  ISTATES, 
3  IDEP,  OBCNAXUSER.  PRINTON,  NITER) 

DO  7  L  *  1,  NODES 


PERF0RM155 
PERF0RM56 
PERFQRN157 
PERF0RN158 
PERFORM! 59 
PERF0RM160 
PERFORMS  1 
PERF0RM62 
PERF0RM63 
COXHQO  23 
CQXNOO  24 
COXNOD  25 
COXHOD  26 
COXMOO  27 
COXNOD  28 
COXMOO  29 
COXHOD  30 
COXMOO  31 
COXMOO  32 
COXNOD  33 
COXMOO  34 
COXHOD  35 
COXHOD  36 
COXHOD  37 
COXNOD  38 
COXMOO  39 
COXHOD  40 
COXHOD  41 
COXHOD  42 
PERFORMS? 
PERF0RM168 
PERF0RM169 
PERF0RM170 
PERFQRM171 
PERF0RM172 
PERF0RM173 
PERFORM  74 
PERFORM  75 
PERF0RN176 
PERF0RM77 
PERF0RM178 
PERF0RM79 
PERFORM80 
PERFORMS 
PERFORM  32 
PERF0RH183 
PERF0RH184 
PERF0RM135 
PERF0RM186 
PERF0RM37 
PERF0RM88 
COXHOD  43 
COXHOD  44 
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oooo—ooooo  ooo 


00  7  H  *  2,  ITYPES 

EXPVFUL.l)  =  EXPVAL(l.l)  4  EXPVAL(L,H) 

7  CONTINUE 

USLTA  »  A0S(P€XPVAL-£XPVAL‘,NDISKS+3i  1)1 
PEXPVAL  >  EXPVflL(NDISKS+3,l> 

NITER  »  NITER  4  t 
FftINTm  a  .FALSE. 

WITE(6,1210>  DELTA 
1210  FORMAT ( A///10X, 'DELTA  =  ",F10.5) 

IF  ((NITER  .LT.  HAXITER)  .AM).  (DELTA  .GT.  T0L5WD1  GO  TO  5 
PRINTON  a  .TRUE. 

CALL  EXPECT(PR08HAR>PR08TOT,PR08PAR,NUS81S»ISTATEli 
INODETYPi EH>VflL. DEPART,  UTIL, SERVICE, E, HAXNCBE, 

2HAXUSEI ,  NAXTYPE,  HIDES,  ITYPES,  IUSERS,  MAXSTAT.  IMJEX3,  ISIZE3,  ISTATES, 
3IDEP,  DEP, HAXUSER,  FRINTON,NITER) 

GO  TO  20 

10  NR1TE(6,15) 

15  FORMAT ( //LOX, ’ERRORS  HAVE  OCCURRED  IN  THE  11PUT  DESCRIPTION.*, 

1//2X, "CHECK  VOUR  CtPUT  DATA  FOR  CORRECTNESS.*) 

20  CONTINUE 

CLOSE  THE  MASS  STORAGE  FILE. 

CALL  CL0SHS(3I 
STOP 
EM) 

SUBROUTINE  IM)ATA(PROS, SERVICE. E.NOOETYP.NUSERS.MAXNODE, 
lHAXUSER.MAXrrPE, MAXSTAT, NODES, IUSBTS, ITVPES, X, INDEX3, ISIZE3, 
2UAREA, ISTATES, IDEP,DB>,EXPVAL,NDISKS<MAX1TER,T0LERAN, 
3DISKFRB,TCPU,OIO,SIZEJOB,AVGJOBS,NAVGJOB),RETURMS(ERROR) 

DIMENSION  FROBIMAWODE.MAXNODE),  9ERVICE(HAXN0CE,HAXTYPE), 
lE(MAXNOOE.MAXTYPE),  NOOETYP(NAXNOOE),  NUSERS(HAXTYPE), 

2X(HAXN00E),  MREA(MAXNODE),  INDEX3USIZE3),  IDEP(HAXNOOE>, 
3DEP(MAXN0DE,mXTYPE,MAXUSBl>,EXPVAL(HAXN0DE,MAXTYPE), 
4DISKPRB(HAXHOOE>HAXTYPE),TCFU(HAXTYPE)iDIO<HAXTYPE), 
SSIZEJOB(HAXTYPE) 

SUEROUTHE  DPUT  READS  THE  DATA  tCCESSARY  TO  RUN  THE 
PROGRAM.  THE  CHARACTERISTICS  OF  THE  TIME-SHARING  COMPUTBT 
SYSTEM  BEING  MODELED  ARE  INPUT  IN  THIS  ROUTINE. 

WRITE (fell) 

FORMAT! /7//10X,  *#*#«THE  INPUT  TO  THE  MODEL  FOLLOWS******) 

READ  THE  WIDER  OF  NODES  IN  THE  1ETN0RK  AND  THE  NUMBER 
OF  TYPES  OF  TASKS. 

READ(5,*)  NODES,  ITYPES 
IF  (  BJF<5)  )  1000,15 
15  CONTINUE 

IF  (  (NODES. GT.  MAXNODE)  .OR.  (ITYPES  .GT.  NAXTYPE)  )  GO  TO  1100 
MRITEI6,20>  NODES,  ITYPES 

20  FORMAT(//10X,"THE  NUMBER  OF  MOOES  REPRESENTED  IN  TIC  HOOEL  ■  *, 


COXHOD 

45 

COX  MOD 

46 

COXHOD 

47 

COXHOD 

48 

COXHOD 

49 

COXHOD 

50 

COXHOD 

51 

COXHOD 

52 

COXHOO 

53 

COXHOD 

54 

COXHOD 

55 

COXHOD 

56 

COXHOD 

57 

COXHOD 

58 

COXHOD 

59 

COXHOD 

60 

PERFQRH19I 

PERF0H1192 

PERF0RH193 

PERF0RH194 

PERF0RH195 

PERFQW1196 

PERF0RM197 

PERF0RM198 

PERFORMS 

PERFORH200 

INDATA 

2 

INDATA 

3 

COXHOD 

61 

COXNGO 

62 

IMJATA 

5 

INDATA 

6 

IMMTA 

7 

COXHOD 

63 

COXHOD 

64 

COXHOD 

65 

INDATA 

9 

IMMTA 

10 

IMMTA 

11 

IMMTA 

12 

IMMTA 

13 

IMMTA 

14 

IMMTA 

15 

IMMTA 

16 

INDATA 

17 

IMMTA 

18 

IMMTA 

19 

COXHOO 

66 

INDATA 

22 

IMMTA 

23 

IMMTA 

24 

IMMTA 

25 

IMMTA 

26 

A-6 


lI3//10XiirTHE  NUMBER  OF  TASK  TYPES  REPRESENTED  IN  TTC  V 

INDATA 

27 

2‘HOCEL  *  M3) 

IMMTA 

28 

READ!5,t)  M)ISKS,SIZE»1,NAXITER,T01ERAN 

COXHOD 

67 

WRITE(6,26>  (BISKS,  SIZHHEH 

COXHOO 

68 

26 

FORMAT! //10X,  "THE  NUMBER  OF  DISKS  REPRESENTED  IN  THE  MODEL  *  V 

COX MOO 

69 

1I3//10X,*THE  SIZE  OF  MEMORY  IN  K  BYTES  =  ",F10.5) 

COXHOO 

70 

IF  (  EOF (5)  .EQ.  1  )  GO  TO  1000 

COXHOO 

71 

HRITE!6,27)  MAXI  TER,  TGLERflN 

COXHOD 

72 

27 

FGRHAT(//10X,"THE  MAXIMUM  NUMBER  OF  ITERATIONS  -  M3//10X, 

COXHOO 

73 

l'THE  TOLERANCE  FOR  CODERS! ON  *  ’,F10.5) 

COXHOD 

74 

IF  (  EDF<5)  .EQ.  1)  9)  TO  1000 

COXHOO 

75 

DO  29  I  *  1,  NDISKS 

COXHOD 

76 

REAOXSr-*’)  (DISKPRBd,J),J=l,  ITYPES) 

coxnod 

77 

IF  (EQFI5)  .EQ.  1)  GO  TO  1000 

COXHOO 

78 

29 

CONTINUE 

COXHOD 

79 

WRITE (6, 30) 

COXHOD 

80 

30 

FORMAT! I/m,  ’DISK  15X,'DI9<  ACCESS  PROBABILITIES  * 

COXHOD 

81 

1"BY  TASK  TYPE*) 

COXHOD 

82 

DO  32  I  *  1.  NDISKS 

COXHOD 

83 

WITE(6,3l)  I,  IDISKPRBdiJ),,M>  ITYPES) 

COXHOD 

84 

31 

FORMAT!  mi,  14, 10X,  10!F10.5> ) 

COXHOD 

85 

32 

CONTINUE 

COXHOO 

86 

READ<5,#)  <NUSERS!J>,TCPU(J),D10(J), SIZE,  1UB! J),  J=1,ITYPES) 

COXHOD 

87 

IUSERS  *  0 

INDATA 

33 

ISTATES  »  1 

INDATA 

34 

JMTEI6.35) 

INDATA 

35 

35 

FORMAT! ////10X. "THE  DISTRIBUTION  OF  TASKS  IN  THE  SYSTEM  ", 

COXHOD 

88 

1-BV  TASK  TYPE  V/OX.'TASK  TYPE", 8X, "NUMBER  OF  TASKS", 8X, 

COXHOO 

89 

2"CPU  TIME/tNTERACTI0N",8X.  "DISK  IO'S/INTERACTI0N*,8X, 

COXHOD 

90 

3“ JOB  SIZE  IN  K  BYTES"/) 

COXHOD 

91 

DO  41  I  =  1,  ITYPES 

COXHOD 

92 

WRITE (6, 40)  I,  NUSERSIII,  TCPU(I),  DIOd),  SIZEJOBd) 

COXHOD 

93 

40 

FORHAT(/20X,I3,20X,I3,13X,F10.5,2(20X,F10.5)) 

COXHOD 

94 

IUSERS  *  IUSERS  +  NUSWStl) 

IMMTA 

41 

ISTATES  »  ISTATES  *  (MJSERSd)  ♦  1) 

INDATA 

42 

41 

CONTINUE 

COXHOO 

95 

VRITE(6,42)  IUSERS 

COXHOD 

96 

42 

FORMAT! //17X, "rOTAL", 21X, 13) 

COXHOD 

97 

AVGJOBS  =  0.0 

COXHOD 

98 

DO  43  I  »  1,  ITYPES 

COXHOD 

99 

AWGJOBS  a  PNGJQ8S  ♦  SIZEJOBd)  *  NUSERS(l) 

COXHOD  100 

43 

CONTINUE 

COXHOD  101 

AVGJOBS  a  AVGJOBS  /  DJSERS 

COXHOO  102 

HAVGJOB  a  SIZE®!  (  AVGJOBS 

COXHOD  103 

URITE<6,44)  AVGJOBS, HAVGJOB 

COXHOD  104 

44 

F0RHAT!///l0X,"flWERAGE  JOB  SIZE  *  ",F10.5,/I0X, 

COXHOD  105 

1 ‘AVERAGE  NUMBER  OF  JOBS  IN  PEHORY  =  M3) 

COXHOO  106 

IF!  (IUSERS  .GT.  MAKUSER)  .OR.  (ISTATES  .GT.  MAXSTAT))  GO  TO  1200 

INDATA 

46 

C 

INDATA 

47 

C 

READ  T)€  NODE  CHARACTERISTICS. 

IMMTA 

48 

c 

IMMTA 

49 

HRITE!6,48)  ITYPES 

IMMTA 

50 

48 

FCRMA7(////10X, "NODE  0*RACTERISTICS*//IOX, "NOOEMX,  *N00E*,5X, 

IMMTA 

51 
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l‘SERVIC£*/10X,"NAIE’'i5X,*TYPE,,»7X»,RATE,,6X,*NOOE  SERVICE  RATES  *, 

INDATA  52 

2*F0R  TASK  TYPES  1  TW0UGHM3) 

IIOATA  S3 

DO  65  I  =  1>  NODES 

HOATA  54 

FEAD<5»»)  NODETYP(I),  IDEP(I),  (SERVICE* I, ITYPES) 

COXHOO  107 

IF(  BQF<5)  )  1000,55 

IIOATA  57 

55 

CONTINUE 

IIOATA  56 

IF  (IDEP(I)  .EQ.  2)  READI5,*)  (DEPd.  1, J) , J=1,HAXUSER) 

COXNOD  108 

IF  (  IDEP(I)  .1C.  3)  GO  TO  59 

HOATA  61 

DO  58  J  =  1.  I TYPES 

IIOATA  62 

READ<5.»  (OEP(I,JiK>iK=l,hAXUSER) 

COXHOO  109 

58 

CONTINUE 

INDATA  64 

59 

CONTINUE 

INDATA  65 

WRITE<6,60)  I.  NOOETYP(I).  IDEP(I),  (SERVICE*!, J) , J=1 , ITYPES) 

INDATA  66 

60 

FORMAT  ( /10K,  14, 5X,  14, 8X,  I4-,  6X,8F10. 5, 10(/41K,  EF10.5) ) 

IIOATA  67 

65 

CONTINUE 

IIOATA  68 

RETURN 

IIOATA  122 

1000 

URITE(6> 1005) 

INDATA  123 

1005 

FORMAT ( //10X,  "***l!IEXPECTED  END  OF  INPUT#**1) 

IIOATA  124 

RETURN  ERROR 

IIOATA  125 

1100 

NRITE(6,U05)  NAXHODE,  NAXTYPE,  NODES,  ITYPES 

IIOATA  126 

1105 

FORMAT(//IOX,"MAX1HJH  NODES  *  ",  13/, "MAXIMUM  TASK  TYPES  »  M3 //, 

IIOATA  127 

1"Y0UR  INPUTS  ARE  M3,"  M  M3,"  RESPECT 1VELY.‘> 

IIOATA  128 

RETURN  ERROR 

IIOATA  129 

1200 

WRITE (6, 1205)  MAXUSER,  MAXSTAT,  IUSERS,  ISTATES 

IIOATA  130 

1205 

FORMAT!  //10X,  'MAXIMUM  NUMBER  OF  TASKS  *  M3.//10X, 

COXHOO  110 

1 ‘MAXIMUM  NUMBER  OF  STATES  »  M5.//10X, 

COXNOD  111 

2"Y0UR  DfUTS  ARE  M3,‘  AW  M5,1  RESPECTIVELY. *) 

COXHOO  112 

RETURN  ERROR 

IIOATA  133 

END 

IIOATA  134 

SUBROUTUE  FUNCT(E, SERVICE, IFACTQR,  ISTATE,FNT,NUSERS,NOOETYP, 

FUNCT  2 

INODES,  IUSERS,  ITYPES,  HAXNOOE,  MAXUSER,  NAXUSE1 , 

FUNCT  3 

2MAXTYPE,  MAXSTAT,  IWEX3,  ISIZE3,  ISTATES,  I DEP.DEP),  RETURNS (ERROR) 

FUNCT  4 

DIMENSION  E(MAXNO[£,HAXTYPE),  SERVICE(NAXNODE,MAXTYPE), 

FUNCT  5 

llFACTOR(HAXUSEl),  ISTATE(MAXTYPE),  FNT(HAXSTAT), 

FUNCT  6 

2NUSERS( NAXTYPE),  WKETYP(NAXNODE),  IWEX3(ISIZE3), 

FUNCT  7 

3HEP(HAXN0DE) ,  DBMMAXNOOE,  NAXTYPE,  MAXUSER) 

FUNCT  8 

DOUBLE  PRECISION  I FACTOR 

FUNCT  9 

C 

FUNCT  10 

C 

THIS  SUBROUTINE  CALCULATES  TIE  VALUES  FOR  TIE  FUNCTIONS 

FUNCT  11 

C 

WHICH  ARE  USED  TO  CALCULATE  TIE  PROBABILITY  THAT  THE  COMPUTER 

FUNCT  12 

C 

SYSTEM  BEGIN  MODELED  IS  IN  A  GIVEN  STATE  AT  EBUILIBRIUH. 

FUNCT  13 

C 

FUNCT  14 

IFACTORU)  *  1 

FUNCT  15 

DO  10  I  «  1,  MAXUSER 

FUNCT  16 

IFACTORU+l)  =  IFACTORU)  »  I 

FUNCT  17 

to 

CONTINUE 

FUNCT  18 

DO  60  I  *  1,  NODES 

FUNCT  19 

INDEX  -  1 

FUNCT  20 

DO  20  J  -  1,  ITYPES 

FUNCT  21 

ISTATE(J)  »  0 

FUNCT  22 

20 

CONTINUE 

FUNCT  23 

25 

CONTINUE 

FUNCT  24 

NUMB  >  1 

FUNCT  25 
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DO  30  J  =  l.  ITYPES 

FUNCT 

26 

NUMBER  *  NUMBER  -»  ISTATEU) 

FUNCT 

27 

30 

CONTINUE 

FUNCT 

28 

J  *  2 

FUNCT 

29 

35 

CONTINUE 

FUNCT 

30 

FVALUE1  »  1. 

FUNCT 

31 

FVALUE2  =  1. 

FUNCT 

32 

DO  38  K  *  1,  I TYPES 

FUNCT 

33 

L  sISTATEIK) 

FUNCT 

34 

FVflLUEl  *  FVALUE1  »  <(E(I>K)  *  SBWICE(I.K))*»U 

FUNCT 

35 

FVALUE2  *  FVMJUE2  t  IFACT0R<L+1) 

FUNCT 

36 

38 

CONTINUE 

FUNCT 

37 

IF  (  NCOETYPd)  .Efl.  3  )  FVALUE2  *  l.  /  FVALUE2 

FUNCT 

38 

IF  (  NOOETYPU)  .1C.  3  )  FVALUE2  =  I FACTOR (NUM8ER)  /  FVALUE2 

FUNCT 

39 

FNT(INOEX)  *  FVAUJE1  «  FVM.IE2 

FUNCT 

40 

IF  (  IEEP(I)  .1C.  2  )  GO  TO  41 

FUNCT 

41 

NUH  *  NU1B0T  -  1 

FUNCT 

42 

IF  (  NUN  .ES.  0  >  GO  TO  41 

FUNCT 

43 

FFF  *  l. 

FUNCT 

44 

DO  40  K  =  1.  NUN 

FUNCT 

45 

FFF  «  FFF  *  DEPIM.K) 

FUNCT 

46 

40 

CONTINUE 

FUNCT 

47 

FNT( INDEX)  =  FHT(INBEX)  /  FFF 

FUNCT 

48 

41 

CONTINUE 

FUNCT 

49 

IF  (  IEEP(l)  .1C.  3  )  GO  TO  44 

FUNCT 

50 

FFF  »  1, 

FUNCT 

51 

00  43  K  »  1.  ITYPES 

FUNCT 

52 

IF  (  ISTATE(K)  .EB.  0  )  BO  TO  43 

FUNCT 

53 

MJM  =  ISTATE(K) 

FUNCT 

54 

DO  42  KK  a  1,  HJfl 

FUNCT 

55 

FFF  »  FFF  *  OEPUiK'MO 

FUNCT 

56 

42 

CONTINUE 

FUNCT 

57 

43 

CONTINUE 

FUNCT 

SB 

FNT(INOEX)  *  FNTI INDEX)  /  FFF 

FUNCT 

59 

44 

CONTINUE 

FUNCT 

60 

IMEX  a  INDEX  ♦  1 

FUNCT 

61 

NUWER  a  MME»  f  1 

FUNCT 

62 

ISTATEU)  «  ISTATEil)  ♦  1 

FUNCT 

63 

IF  I  ISTATEU)  .IE.  MJSERS(l)  )  GO  TO  35 

FUNCT 

64 

45 

CONTINUE 

FUNCT 

65 

ISTATEIJ)  «  ISTATEU)  ♦  1 

FUNCT 

66 

IF  (  ISTATEU)  .OT.  NUSERS(J)  )  DO  TO  55 

FUNCT 

67 

K«  J-  1 

FUNCT 

68 

DO  50  L  «  1.K 

RJNCT 

69 

ISTATEIL)  »  0 

FUNCT 

70 

50 

CONTINUE 

FUNCT 

71 

GO  TO  25 

FUNCT 

72 

55 

CONTINUE 

FUNCT 

73 

J  *  J  ♦  1 

FUNCT 

74 

IF  (  J  .LL  ITVPES  )  00  TO  45 

FUNCT 

75 

C 

FUNCT 

76 

C 

TIC  FUNCTION).  VALUES  CALCULATED  FOR  Ml  POSSIBLE  STATES 

FUNCT 

77 

C 

OF  EACH  NODE  ME  SAWED  SO  THAT  TREY  CAN  BE  USED  TO  CALCULATE 

FUNCT 

78 
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c 

THE  NORMALIZATION  CONSTANT. 

FUNCT 

n 

c 

FUNCT 

80 

CALL  URITMS(3.FNT(i)>  {STATES.  I) 

FUNCT 

81 

60 

CONTINUE 

FUNCT 

82 

fCTJRN 

FUNCT 

83 

EM) 

FUNCT 

8* 

FUNCT 

85 

SUBROUTINE  NORMAL (FNT.NORNCDN.NUSERS. 1STATE1 • 

NORMAL 

2 

USTATE2.  ISKIP.IVtfNODE.MAXSTAT, NODES.  ITYPES, 

NORHAL 

3 

2HAXTYPE, INDEX3. ISIZE3.  ESTATES) 

NORMAL 

4 

DIMENSION  FNT(NAXSTAT).  NQRHC0N(HAXSTAT,2), 

NORHAL 

5 

lNUSERSlMAITYPE),  1STATEKHAXTYFE),  ISTATE20WXTYPE). 

NORMAL 

6 

2ISKIP<HAXTYP£) , IM)EX3( ISIZE3) 

NORMAL 

7 

INTEGER  FIWEX 

NORMAL 

8 

REAL  NORHCON 

NORMAL 

9 

C 

NORMAL 

10 

C 

THIS  SUBROUTINE  CALCULATES  THE  NORMALIZATION  CONSTANT 

NORHAL 

11 

C 

FOR  THE  MODEL.  THE  NORMALIZATION  CONSTANT  ASSURES  THAT  THE 

NORMAL 

12 

C 

PROBABILITY  OF  THE  SYSTEM  BEING  IN  All  STATES  SUNS  TO  UNITY. 

NORMAL 

13 

C 

NORHAL 

14 

CALL  R£AQ1S(3i  NORHCON!  1.1MSTATES. I) 

NORMAL 

15 

CALL  URITNS(3.NQRHCGN<1.1).ISTATES.HAXNOOEM) 

NORMAL 

16 

N0RH1  *  1 

NORHAL 

17 

N0RH2  -  2 

NORHAL 

18 

DO  110  I  *  2.  NODES 

NORHAL 

19 

CALL  REA0HS(3.FNT(l).ISTA1ES.l) 

NORMAL 

20 

INDEX  »  1 

NORMAL 

21 

DO  15  J  «  1,  ITYPES 

NORHAL 

22 

ISTATEKJ)  *  0 

NORHAL 

23 

15 

CONTINUE 

NORMAL 

24 

20 

CONTINUE 

NORMAL 

25 

K  «  2 

NORMAL 

26 

DO  25  J  =  l,  ITYPES 

NORHAL 

27 

ISTATE21J)  *  0 

NORMAL 

28 

25 

CONTINUE 

NORMAL 

29 

KK  *  l 

NORMAL 

30 

FINDEX  *  l 

NORMAL 

31 

LINDEX  *  INDEX 

NORMAL 

32 

I7EMP  =  1 

NORMAL 

33 

DO  30  J  *  1,  ITYPES 

NORMAL 

34 

ISKIP(J)  *  ITEMP  ♦  (MJSERSIJ)  -  ISTATE1UH 

NORHAL 

35 

ITEM*  »  ITEHP  *  (NUSERS(J)  ♦  1) 

NORMAL 

36 

30 

CONTINUE 

NORMAL 

37 

ISXIP(l)  =  ISKIP(i)  ♦  t 

NORMAL 

38 

FVALUE  »  0. 

NORMAL 

39 

40 

CONTINUE 

MflflMAI 

WnwWm 

40 

FVALUE  -  FVALUE  <-  NORMCONIFINDEX.NORNl )  *  FNT(LINDEX) 

NORMAL 

41 

IF  (  FIHJEX  .89.  INDEX  )  GO  TO  70 

NORMAL 

42 

ISTATE21KK)  »  ISTATE2000  ♦  l 

NORMAL 

43 

IF  <  ISTATE2CKK)  .ST.  1STATEKKK)  )  GO  TO  50 

NORMAL 

44 

FINDEX  »  FINDEX  ♦  1 

NORMAL 

45 

LIWEX  a  LINOEX  -  1 

NORMAL 

46 

SO  TO  AO 

NORMAL 

47 
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SO 

CONTINUE 

NORMAL 

48 

FINDEX  *  F INDEX  ♦  ISKIP(KK) 

NORMAL 

49 

UNDO  >  LIIDEX  -  ISKIPfKK) 

NORMAL 

50 

KK  «  KK  *  l 

NORMAL 

51 

ISTATE2(KX)  =  ISTATE2(KX)  +  1 

NORMAL 

52 

IF  (  ISTATE2<KK)  .GT.  lSTATEl(KK)  )  GO  TO  50 

NORMAL 

53 

L  ■  KK  -  1 

NORMAL 

54 

DO  60  J  *  1,  L 

NORMAL 

55 

ISTATE2(J)  *  0 

NORMAL 

56 

60 

CONTINUE 

NORMAL 

57 

KK  «  1 

NORMAL 

5B 

GO  TO  40 

NORMAL 

59 

70 

CONTINUE 

NORMAL 

60 

NORKONIINIEXiNOfVtt)  =  FVALUE 

NORMAL 

61 

INQEX  «  INDEX  ♦  1 

NORMAL 

62 

ISTATEKU  *  ISTATE1  (1)  +  1 

NORMAL 

63 

IF  (  ISTBTEKl)  .LE.  NUSERS(l)  )  GO  TO  20 

NORMAL 

64 

80 

CONTINUE 

NORMAL 

65 

ISTATEl(K)  *  ISTATEl(K)  +  1 

NORMAL 

66 

IF  (  ISTATEl(K)  .GT.  NUSERSIK)  )  GO  TO  100 

NORMAL 

67 

L  »  K  -  1 

NORMAL 

68 

00  90  >  1.  L 

NORMAL 

69 

ISTATEl(J)  =  0 

NORMAL 

70 

90 

CONTINUE 

NORMAL 

71 

GO  TO  20 

NORMAL 

72 

100 

CONTINUE 

NORMAL 

73 

K  a  K  ♦  1 

NORMAL 

74 

IF  (  K  .LE.  ITYPES  )  GO  TO  80 

NORMAL 

75 

CALL  URITHS(3>NGRMC0N(  1 » N0RH2)  *  ISTATESi  MAXNOEE+I ) 

NORMAL 

76 

ITENP  *  N0R11 

NORMAL 

77 

N0RN1  «  N0RH2 

NORMAL 

78 

N0RH2  >  ITENP 

NORMAL 

79 

110 

CONTINUE 

NORMAL 

80 

RETURN 

NORMAL 

81 

END 

NORMAL 

82 

SUBROUTINE  MARGIN<FNT,CONNORM,PROBMAR,PROBCUH, 

MARGIN 

2 

INUSERS. ISTATEli ISTATE2. ISKIP.MAXNOOE. NAXSTAT, 

MARGIN 

3 

2MAXTYPE,  NODES,  ITYPES.IN0EX3, 1SIZE3,  ISTATES) 

MARGIN 

4 

DIMENSION  FNT( NAXSTAT),  CONNORNf NAXSTAT, 2), 

MARGIN 

5 

1PR0BHMHNAXSTAT),  PROKUHIHAXSTAT), 

MARGIN 

6 

2NUSERS(MAXTYTE)  •  ISTATEl(NAXTYPE).  ISTATE2(HAXTYPE)> 

MARGIN 

7 

3ISKIP(NAXTYPE),  IWB3(ISIZE3) 

MARGIN 

8 

INTEGER  FINOEX 

MARGIN 

9 

C 

MARGIN 

10 

C 

THIS  9JBR0UT1*  CALCULATES  THE  PROBABILITY  THAT  A  GIVEN 

MARGIN 

11 

C 

NODE  IS  IN  A  GIVB4  STATE.  THIS  PROBABILITY  IS  CALCULATED 

MARGIN 

12 

C 

FOR  ALL  POSSIBLE  STATES  HR  EACH  NODE. 

MARGIN 

13 

C 

MARGIN 

14 

NORM  a  i 

MARGIN 

15 

N0RN2  a  2 

MARGIN 

16 

CALL  AEAM(3iCONNMHliNORHl)iISTATESiMAXNQOE+NQ()ES) 

HARBIN 

17 

CONSTAN  >  t.  /  COMORHI  ISTATESi  N0RH1 ) 

HARDIN 

18 

-n 

i 

■ 

MARGIN 

19 

A-l  1 


UNDO  «  ISTATES 

HARDIN 

20 

PPP  »  0. 

MARGIN 

21 

CALL  REA0NS(3>FNT(l)i ISTATES* NODES) 

MARGIN 

22 

CALL  READMS<3.  CONHOm<  1 ,  NQRH2> .  ISTATES.  HAXNQOE^CnES-i  > 

MARGIN 

23 

CONTINUE 

MARGIN 

24 

ffiOGNARiFlMEX)  *  CONST  AN  •  CCWCGMLIHO.NCRN2)  • 

MARGIN 

25 

lFNT(FlMEX) 

MARGIN 

26 

PRQDCUMIFINOO)  =  CONST AN  •  COWORH<FINOO.NQRH2)  * 

MARGIN 

27 

lFNT(UIOO) 

MARGIN 

28 

PPP  »  PPP  f  PWBHARIFINDQ) 

MARGIN 

29 

FINDO  *  F1HDO  +  1 

MARGIN 

30 

UNDO  *  L  INDEX  -  1 

MARGIN 

31 

IF  (  FINDO  .LE.  ISTATES  )  GO  TO  20 

MARGIN 

32 

IF  (  NODES  .LE.  2  >  GO  TO  150 

MARGIN 

33 

NCYCLES  »  NODES  -  2 

MARGIN 

34 

CALL  WmmPROEHfttU). ISTATES, MAXNOOE*2+NOOES) 

MARGIN 

35 

DO  140  I  >  1.  NCYCLES 

MARGIN 

36 

NSUB  *  NODES  ~  I 

MARGIN 

37 

CALL  READMS(3.FMT(1)*1STATES*NSUB) 

MARGIN 

38 

ITENP  »  N0RN1 

MARGIN 

39 

NORHl  -  NORM 

MARGIN 

40 

NORM  a  ITENP 

MARGIN 

41 

CALL  READNS(3f  CONORM  1  >  NORM)  >  ISTATES.  NAXNOOE+NSUIH ) 

MARGIN 

42 

INDEX  a  i 

MARGIN 

43 

DO  30  J  «  l.ITYPES 

MARGIN 

44 

ISTATEl(J)  »  0 

MARGIN 

♦5 

CONTONE 

MARGIN 

46 

K»2 

MARGIN 

47 

PPP  »  0. 

MARGIN 

48 

CONTINUE 

MARGIN 

49 

DO  50  J  »  1,  I TYPES 

MARGIN 

50 

ISTATE2IJ)  >  ISTATEl(J) 

MARGIN 

51 

CONTINUE 

MARGIN 

52 

FIWO  -  L 

MARGIN 

53 

UNDO  >  INDEX 

MARGIN 

54 

KK  ■  1 

MARGIN 

55 

ITENP  a  1 

MARGIN 

56 

DO  DO  J  «  1,  ITYPES 

MARGIN 

57 

ISKIP(J)  «  ITBP  *  ISTATEl(J) 

MARGIN 

SB 

ITENP  a  ITBP  t  (NUSERSU)  *  1) 

MARGIN 

59 

CONTINUE 

MARGIN 

60 

ISKIP(l)  «  ISKIP(l)  ♦  1 

MARGIN 

61 

FVWJUEt  •  0. 

MARGIN 

62 

FVALUE2  >  0. 

MARGIN 

63 

CONTINUE 

MARGIN 

64 

FVAUUEl  -  FVALUE1  ♦  PRODCUN(UNDO)  • 

MARGIN 

65 

l«OONNMN(FINDOiNORH2)  *  RTT<  INDEX))  / 

MARGIN 

66 

2CONOIMLINDO.NORM)) 

MARGIN 

67 

FVALUE2  a  RJALUE2  f  PROBCON(LINDO)  » 

MARGIN 

68 

1(  (COMIORNdlOOtNQRM)  •  FNT(FINDO))  1 

MARGIN 

69 

2CONCRN(LINDO.NORt1t) ) 

MARGIN 

70 

IF  (  UWO  .ED.  ESTATES  )  GO  TO  100 

MARGIN 

71 

ISTATE20Q0  «  1STKTE2(KX)  +  l 

MARGIN 

72 

A- 12 


I 


90 


90 


100 


110 


120 


130 


HO 

150 


160 


IF  (  ISTATE2(KK)  .GT.  NUSERS(KK)  )  GO  TO  80 

MARGIN  73 

FIMEX  *  FINCEX  +  1 

MARGIN  74 

LINDEX  >  L1NBEX  f  1 

MARGIN  75 

GOTO  70 

MARGIN  76 

COKTINLE 

MARGIN  77 

FINDEX  *  FIMEX  ♦  ISKIPIKX) 

MARGIN  78 

UMEX  *  LIMEX  +  ISXIP(KK) 

MARGIN  79 

KK  *  XX  +  1 

MARGIN  80 

ISrflTE2(XXI  *  KTATEHHO  +  1 

MARGIN  81 

IF  (ISTHiTE2(KK)  .GT.  NUSERS(KX)  )  GO  TO  80 

MARGIN  82 

L  *  KK  -  1 

MARGIN  83 

DO  90  J  *  It  L 

MARGIN  84 

ISTATE2U)  *  ISTATEHJ) 

MARGIN  85 

CONTINUE 

MARGIN  86 

KK  *  1 

MARGIN  87 

GO  TO  70 

MARGIN  88 

CONTINUE 

MARGIN  89 

PRODMARI  INDEX)  =  FVALUE1 

MARGIN  90 

PPP  »  PPP  *■  FVALUEl 

MARGIN  91 

PROBCUMUMEX)  ■  FVALUE2 

MARGIN  92 

IMEX  «  INDEX  +  1 

MARGIN  93 

ISTATEKl)  *  ISTATEl(l)  *  1 

MARGIN  94 

IF  (  ISTATEKl)  .LE.  NUSERS(l)  )  GO  TO  40 

MARGIN  95 

CONTINUE 

MARGIN  96 

ISTAIEl(K)  »  ISTATEl(K)  ♦  1 

MARGIN  97 

IF  (  ISTATEKK)  .GT.  NUSERSIK)  )  GO  TO  130 

MARGIN  98 

L  «  X  -  1 

MARGIN  99 

DO  120  J  «  It  L 

MARGIN  100 

ISTAIEKJ)  »  0 

MARGIN  101 

CONTINUE 

MARGIN  102 

X  ■  2 

MARGIN  103 

GO  TO  40 

MARGIN  104 

CONTINUE 

MARGIN  105 

X  »  X  ♦  1 

MARGIN  106 

IF  (  X  .LE.  ITYFES  )  00  TO  110 

MARGIN  107 

CALL  URITMS(3t  PROBMARt  l )  >  ISTATES.  HAXNQDE*2+NSUB) 

MARGIN  106 

CONTINUE 

MARGIN  109 

CONTINUE 

MARGIN  110 

PPP  »  0. 

MARGIN  111 

DO  160  I  ■  1.  ISTATES 

MARGIN  112 

PPP  »  PPP  ♦  PROBOUH(I) 

KARGIN  113 

CONTINUE 

MARGIN  114 

CALL  MRITHS(3>PR0BCUH(  l ) .  ISTATES»HAXN0DE*2+ 1 ) 

MARGIN  115 

RETURN 

HARDIN  116 

END 

MARGIN  117 

SUBROUTINE  EXPECT(PR06MAR.PflOfiTOT,PROBPAR,NUSERS.ISTATE. 

EXPECT  2 

lNXETYPt  EXWALi  DEPART  i  UTILt  SERVICE.  Ei 

EXPECT  3 

2NAXN0DE.HAXUSEl,NAXTYPE.  NODES.  ITYPES.  IUSERStNAXSTATt  INDEX3. 

EXPECT  4 

3ISIZE3*  ISTATES*  IDEP,  DEP,  HAXUSER*  PRINTON.  NITER) 

COXmO  113 

DIMENSION  PROB4AR(NAXSTAT),PWBTOT(HAXNOOE,HAXUSE1). 

EXPECT  6 

IPROBPMHMUNOOEiNAXPfPEiNAXUSEl)*  NUSERS(HAXTYPE), 

EXPECT  7 

2!STATE(NUTTFE)>  NOETYP(NAXMOE)  *  EXPWt(HAXNODE,HAXTYPE), 

EXPECT  8 

3DEPART(NAXNOOE*NAnVPE)t  UTH(NAXNOOE.HAXTYPE). 

EXPECT  9 

1 


i  A- 13  « 

i  i 


<->  o  o  o  o  o  o 


4SERVICE(NftXH00E,MAXTYPE),  1N0EX3(1SIZE3>,  E(MAXNODE,HAXTYPE), 
5I0EP(HAXN00E) »  DEP(MA)(NOOE,  HAXTYPE.  HAXUSER) 

LOGICAL  PRINTQN 

I»1S  SUBROUTINE  CALCULATES  VARIOUS  PERFORMANCE  (CASURES 
FOR  THE  TlfE-SHARIHS  SYSTEM  BEING  MOOELED.  TIESE  MEASURES 
INCLUDE  DC  MARGINAL  QUEUE  LENGTH  DISTRIBUTION  FOR  EACH  TYPE 
OF  USER  AT  EACH  NODE,  TIE  UTILIAZTION  OF  THE  NODES  BY  EACH 
TYPE  OF  USER,  AND  THE  RESPONSE  THE  FOR  EACH  TYPE  OF  USER. 

DO  10  I  *  l,  MMNOOE 
DO  10  J  *  l.  MAXUSE1 
PROBTOT(ItJ)  *  0.0 
DO  10  K  *  1,  MACTYPE 
PROBPAR(I.K.J)  «  0.0 
UTIL(I.IC)  »  0.0 
DEPARTd.K)  *  0.0 
10  CONTINUE 

IF  (PRINTON)  WRITE (6,1) 

1  FORMAT!  1H1,  ////'10X,  "«**»fHE  OUTPUT  OF  THE  MODEL  FOLLONS#*****) 

DO  70  1  *  1,  NODES 

CALL  RH0HS(3.PR0GMAR<1).ISTATES,NAXN0DE«2«I) 

K  =  1 
INDEX  *  0 
NTOTAL  *  0 
DO  20  J  *  1,  ITYPES 
ISTATE(U)  *  l 
20  CONTINUE 
40  CONTINUE 

NTOTAL  -  NTOTAL  *  1 
INDEX  *  IMEX  ■»  1 
PROG  -  PMDHARdNDEX) 

PROBTOTII, NTOTAL)  *  FROGTOTd, NTOTAL)  4  PRO! 

DO  50  J  *  1,  ITYPES 
M  •  ISTATE(J) 

PRODPAR(I,JtJJ)  -  PRO0PAR(I,J,JJ)  4  PROS 
IF  (  JJ  .EQ.  1H30  TO  50 
PARTIAL  «  TOE  *  (JJ  -  1) 

IF  (  NODETYPd)  .EQ.  3  >  GO  TO  45 
PARTIAL  *  PARTIAL  /  (NTOTAL  -  1) 

UTlLd.J)  »  UTlLd.J)  ♦  PARTIAL 
45  CONTINUE 

PARTIAL  »  PARTIAL  *  11.  /  SERVICEd.J)) 

DEPART (1,J)  «  DEPART(IiJ)  4  PARTIAL 
50  CONTINUE 

IF  (  I  WEI  .EQ.  ISTATES  >  GO  TO  70 
60  CONTINUE 

ISTATE(K)  »  ISTATE(K)  4  1 
IF  1  ISTATE(K)  .GT.  NUSERS(K)  4  1  )  GO  TO  65 
K  >  1 
GO  TO  40 
65  CONTONE 

ISTATE(K)  >  l 


EXPECT  10 
EXPECT  11 
COXHOD  114 
EXPECT  12 
EXPECT  13 
EXPECT  14 
EXPECT  15 
EXPECT  16 
EXPECT  17 
EXPECT  18 
COXHOD  115 
COXHOD  116 
COXHOD  117 
COXHOD  118 
COXHOD  119 
COXHOD  120 
COXHOD  121 
COXHOD  122 
COXHOD  123 
COXHOD  124 
EXPECT  21 
EXPECT  22 
EXPECT  23 
EXPECT  24 
EXPECT  25 
EXPECT  26 
EXPECT  27 
EXPECT  28 
E»en  29 
EVECT  30 
EVECT  31 
EVECT  32 
EVICT  33 
EVECT  34 
-  JBT  35 
EVICT  4 
EVT  ft 

eve*  a 
evcci  a 

EXPECT  « 
EVECT  41 
EVKT  42 
EVECT  43 
EVECT  44 
EVECT  45 
EVECT  46 
EXPECT  47 
EVECT  « 
EVECT  49 
EVECT  a 
EVECT  51 
EXPECT  52 
EXPECT  S3 
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NTOTAL  *  NTOTAL  -  NUS®S(K> 

EXPECT  54 

K  s  K  +  1 

EXPECT  55 

IF  (  K  .LE.  ITYPES  )  GO  TO  60 

EXPECT  56 

70 

CONTINUE 

EXPECT  57 

10SERS1  *  IUSERS  +  1 

EXPECT  50 

IF  (PRINTON)  HRITEI6.74)  NODES 

COXHOO  125 

74 

FORHAT</#'//10X» “TVC  MARGINAL  QUEUE  LENGTH  PftOBABILITlES  *, 

EXPECT  60 

1-AT  EQUILIBRIUHVmOX, -QUEUE  LENGTH', 5X,  ‘QUEUE  LENGTH 

EXPECT  61 

2-PROBABILITY  FOR  NODES  1  THROUGH  M3) 

EXPECT  62 

DO  75  I  *  1,  NODES 

EXPECT  63 

EXPVALd,l>  =0. 

EXPECT  64 

75 

CONTINUE 

EXPECT  65 

DO  110  J  *  1,  1USERS1 

EXPECT  66 

JJ  *  J  -  1 

EXPECT  67 

DO  90  I  *  1,  NODES 

EXPECT  68 

EXPVAL(M)  *  EXPVAL(I,1)  ♦  JJ  »  PROBTOT(I.J) 

EXPECT  69 

90 

CONTINUE 

EXPECT  70 

IF  (FRINTON)  HR1TE(6,80)  JJ,  (FROBTOTII.J),  1=1, NODES) 

COXHOO  126 

80 

FORMAT! /14X, 13,77, 10F10.5, 10(/24X, 10F10.5) ) 

EXPECT  72 

no 

CONTINUE 

EXPECT  73 

IF  (PRINTGN)  WnE(6,ni) 

COXNOD  127 

111 

FORMAT ( / // /10X,  ‘THE  EXPECTED  NUMBER  OF  TASKS  AT  EACH  N0DEV//1OX, 

EXPECT  75 

l'NOOEMOX, ‘EXPECTED  VALUE") 

EXPECT  76 

DO  113  I  *  1.  NODES 

EXPECT  77 

IF  (FRINTON)  UR1TE(6,112)  I,  EXPVAL(l.l) 

COXNOD  128 

112 

F0RMAT(/10X,I4,10X,F10.5) 

EXPECT  79 

113 

CONTINUE 

EXPECT  80 

IF  (PRINTON)  URITE(6,114) 

COXNOD  129 

114 

FORMAT! /(7/7L0X, -THE  MARGINAL  QUEUE  LENGTH  PROBABILITIES  ", 

EXPECT  82 

l’AT  EQUILIBRIUM  CLASSIFIED  BY  TASK  TYPE") 

EXPECT  83 

DO  160  J  ■  1.  ITYPES 

EXPECT  84 

IF  (PRINTON)  URITE(6, 115)  J,  NODES 

COXNOD  130 

115 

FORMAT! //7/10X,  "THE  MARGINAL  QUEUE  LENGTH  PROBABILITIES  ", 

EXPECT  86 

1-AT  EQUILIBRIUM  FOR  TASK  TYPE',  I3///10X, ‘QUEUE  LENGTH*, 5X, 

EXPECT  87 

2-QUEUE  LENGTH  PROBABILITIES  FOR  NODES  1  THWXJOH  M3) 

EXPECT  88 

IUSERS1  *  NUSERS(J)  ♦  1 

EXPECT  89 

DO  116  1  »  1,  NODES 

EXPECT  90 

EXPVAL(I.J)  *  0. 

EXPECT  91 

116 

CONTINUE 

EXPECT  92 

DO  150  K  «  1,  IUSERS1 

EXPECT  93 

KK«K-  l 

EXPECT  94 

DO  120  1  *  1,  NODES 

EXPECT  95 

EXPVAL(I,J)  «  EXPVAL(I.J)  +  KK  «  PROBPARd,J,K) 

EXPECT  96 

120 

CONTINUE 

EXPECT  97 

IF  (PRINTON)  UR1TEC6.130)  KK,  (PROBPAR(I,J,K),  1-1, NOTES) 

COXHOO  131 

130 

FORMAT  ( /14X,  13, 6X,  IOF10.5, 10(/24X,  10F10.5) ) 

EXPECT  99 

150 

CONTINUE 

EXPECT  100 

160 

CONTINUE 

EXPECT  101 

HRITE(6,162)  ITYPES 

EXPECT  102 

162 

F0RMAT(///('10X,-THE  EXPECTED  NUMBER  OF  TASKS  AT  EACH  NODE  V 

EXPECT  103 

1 ‘CLASSIFIED  BY  TASK  TYPE* //dOI, ‘NODE",  10X, -EXPECTED  NUMBER  -, 

EXPECT  104 

2-OF  TASK  TYPES  1  THOUGH  M3,*  FOUM)  AT  NODE*) 

COXHOO  132 

00  166  I  >  1,  NODES 

EXPECT  106 
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WRITER,  164)  I,  (EXPVAL(I.d)  „M,  ITYPES) 

EXPECT  107 

164 

FCRKATdlOX, 14, 10X. 10F10.5, 10( /24X, 1QF10.5) ) 

EXPECT  108 

166 

CONTINUE 

EXPECT  109 

IF  (PRINTON)  WRITE16, 170)  1TYPES 

COXNOD  133 

170 

FORHAT(////1M,HNODE  UTILIZATION  BY  EACH  TASK  TYPEV//10X, 

EXPECT  111 

l’NME'i  10X,  "UTLIZATION  OF  NODE  BY  TASK  TYPE  t  THROUGH  M3) 

EXPECT  112 

DO  180  I  3  1,  NODES 

EXPECT  113 

IF  (  MOOETVP(I)  .EQ.  3  )  GO  TO  180 

EXPECT  114 

IF  (PRINTGN)  UR1TE16, 175)  l,  (UTIUI*d)„M,  ITYPES) 

COXNOD  134 

175 

F0RHAT(/10X,  14, 10X.10F10.5, 10I/24X,  10F10.5)) 

EXPECT  116 

180 

CONTINUE 

EXPECT  117 

DO  190  I  3  1,  NOTES 

EXPECT  118 

DO  185  J  3  2,  ITYPES 

EXPECT  119 

UTILd.l)  3  UTILd.l)  ♦  UTlLd.d) 

EXPECT  120 

185 

CONTINUE 

EXPECT  121 

190 

CONTINUE 

EXPECT  122 

IF  (PRINT®!)  WUTE(6,195) 

COXNOD  135 

195 

FORNAT(////10X, ’TOTAL  NODE  UTILIZATION’ ///10X, 

EXPECT  124 

fNOOEMOI ,' " UTILIZATION* ) 

EXPECT  125 

DO  200  I  3  1,  MHES 

EXPECT  126 

IF  1  NOOETYPd)  .ED.  3  )  GO  TO  200 

EXPECT  127 

IF  (PRINTON)  URITEI6, 198)  I.  UTILd.l) 

COXNOD  136 

198 

F0RNAT(/10X,  14, 10X.F10. 5) 

EXPECT  129 

200 

CONTINUE 

EXPECT  130 

IF  (PRINTON)  HRITE(6,210)  ITYPES 

COXNOD  137 

210 

FORWnT///lOK,’lHE  MEAN  Tilt  EACH  TASK  TYPE  SPENDS 

EXPECT  132 

1’AT  THE  VWI0U3  NODES  FOR  EACH  lNTBttCTIONV//10X,*MflEMOX, 

EXPECT  133 

21CAN  Tilt  FOR  TASK  TYPE  1  THROUGH  M3) 

EXPECT  134 

DO  230  I  3  1,  NODES 

EXPECT  135 

IF  (IDffd)  .ED.  1)  00  TO  205 

EXPECT  136 

IF  (PRINTON)  URITE(6,204) 

COXNOD  138 

204 

FORNAT(/10I.'THE  MODEL  DOES  NOT  CALCULATE  THIS  TINE  FOR  ",/10K, 

EXPECT  138 

1 ’NODES  MITH  STATE  DEPENDENT  SERVICE  RATES.') 

EXPECT  139 

GO  TO  290 

EXPECT  140 

205 

CONTINUE 

EXPECT  141 

DO  220  J  3  1,  ITYPES 

EXPECT  142 

UTlLd.d)  3  0. 

EXPECT  143 

IF  <  (OEPWTd.d)  .ED.  0.)  .OR.  (EU,d)  .ED.  0.)  )  GO  TO  220 

EXPECT  144 

UTlLd.d)  3  EXPVAL(I.J)  /  DEPART(I.d)  *  Ed,J> 

EXPECT  145 

220 

CONTINUE 

EXPECT  146 

IF  (PRINTON)  MRITE(6>225)  I,  (UTIL( I, J),>1, ITYPES) 

COXNOD  139 

225 

F0RNAT(/10X,  14, 10X,  1CF10.5, 10I/24X,  I0F10.5)) 

EXPECT  148 

230 

CONTINUE 

EXPECT  149 

IF  (PRINTON)  URITE(6,240) 

COXNOD  140 

240 

FORNATd////10X,’TT€  YtAH  RESPONSE  TINE  FOR  T)E  VARIOUS  ", 

EXPECT  151 

ITASK  TYPES’) 

EXPECT  152 

DO  290  I  3  1,  NODES 

EXPECT  153 

IF  (  NOOETYP(l)  ,NE.  3  )  GO  TO  280 

EXPECT  154 

IF  (PRINTON)  HRITE(6,250)  I 

COXNOD  141 

250 

FORMAT! //7/lOX,’T)€  MEAN  RESPONSE  Tilt  FOR  TASK  ", 

EXPECT  156 

1 ’REPRESENTED  AT  NODE  ",  I3//10X,'TASK  TYPEMOX, 

EXPECT  157 

2’RESPONSE  Tilt’) 

EXPECT  158 

DO  270  J  «  l,  ITYPES 

EXPECT  159 
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TINE  »  0. 

IF  (  DEPART ( I. J)  .EQ.  0.)  GO  TO  255 

TIK  *  INUSERSIJ)  -  EXPVALt I t J) )  /  DEPART C I . J > 

255  CONTIME 

IF  (FRINTUN)  URITE(6,260)  J,  TINE 
260  FQRHAT(/13X,I3,14X,F10.5) 

270  CONTINUE 
280  CONTINUE 
RETURN 
C 

C  THAT'S  ALL  FOLKS!! 

C 

EW 

SUBROUTINE  CHEM75(PR(».NDI3t(S»  DISKPRB. TCPU,  DIO,  NAVGJOB. 

1EXPVAL,  SERVICE. DELTAi PEXPVAL, IUSERS.HAXNODE.NAXTYPE. 
2ITYPES.NOOES, NITER, E.X.UAREA) 

DIMENSION  PR08(NAXN0OE,HAXN0DE),DISKPRB(HAXN0DE,HAXTYPE), 
lDIO(NAXTVPE) , EKPVflL(NAXNODE> HAXTYPE) , SERVICE (NAXNODE, HAXTYPE) « 
2£(MAXNO0£,HAXTYPE> , X(NAXNOOE) , TCPU< HAXTYPE) , 

2UAREA! NAXNODE) 

REAL  DELTA, PEXPVAL 

INTEGER  NDISKS, NAXNODE, HAXTYPE, IUSERS 
C 

C  THIS  SUBROUTINE  CALCULATES  THE  DELTA  BETWEEN  THE  LAST 

C  ITERATIONS  ANSWER  FOR  THE  NUICER  OF  JOBS  IN  THE  "THINK  STATE" 

C  AND  THE  CURRENT  ITERATIONS  ANSWER  CONTAINED  IN 
C  EXPVfiL(N0ISKS+3,l).  ALSO,  USING  THE  SUBROUTINE  CONFUTE  IT 
C  RECOMPUTES  THE  PROBABILITY  TRANSITION  »TRIX,  PROS,  USING  TIC 
C  FORMULAE  CONTAINED  IN  CHETS  ARTICLE,  “QUEUEING  NETWORK  MODEL 
C  OF  INTERACTIVE  COHWIMG  SYSTEMS",  IN  THE 

C  PROCEEDINGS  CF  THE  IEEE,  VOL.  63,  NO.  6,  JUtE  75. 

C  - 

C 

UR1TE(6,67)  NITER 

67  F0RMAT(1HI,////10X, "TRANSITION  PROBABILITIES  DESCRIBING  ", 
llttVOENT  OF  TASKS  AMONG  THE  MODEL'S  NODES  FOR  ITERATION  *,I5) 
DO  110  K  «  l,  ITYPES 

WRITE (b, 68)  K,  NODES 

68  F0RNKK///10X, "TRANSITION  PROBABILITIES  FOR  TASK  TYPE  ", 
113/mX,  "DEPARTURE", 7X,  "PROBABILITY  OF  TASK  HOIBENT  FWKV12X, 
2*NQ0E".9X,  "DEPARTURE  HIDE  TO  NODES  1  THROUGH  M3) 

CALL  C0NPUTE(K,PRO8,  NDISKS,  D1SXPRB,TCPU,  DIO. 

1EXPVAL, SERVICE, IUSERS, HAXNOOE, HAXTYPE, NODES, NAVGJOB) 

IF  (NITER  .NE.  0)  00  TO  70 
PROB(l.I)  *  PROB(l.l)  4  PR0B(N0ISKS+2,1) 

PR0B(»ISKS«-2,1)  «  0.0 
70  CONTONE 

DO  75  I  •  l,  NODES 

HRITE(6i74)  I,  (PROD!  J,  I ),J«1, NODES) 

74  F0RNAT(/12X, 14, 10X, 10FI0.5, 101 / 24X, 10F10.5) ) 

75  CONTINUE 

NUMPROB  ■  NAXNODE  4  NAXNODE 


EXPECT  160 
EXPECT  161 
EXPECT  162 
EXPECT  163 
COXHOD  142 
EXPECT  165 
EXPECT  166 
EXPECT  167 
EXPECT  168 
EXPECT  169 
EXPECT  170 
EXPECT  171 
EXPECT  172 
COXHOD  143 
COXHOD  144 
COXHOD  145 
COXHOD  146 
COXHOD  147 
COXHOD  148 
COXHOD  149 
COXHOD  150 
COXHOD  151 
COXHOD  152 
COXHOD  153 
COXHOD  154 
COXHOD  155 
COXHOD  156 
COXHOD  157 
COXHOD  158 
COXHOD  159 
COXHOD  160 
COXHOD  161 
COXHOO  162 
COXHOD  163 
COXHOD  164 
COXHOD  165 
COXHOD  166 
COXHOD  167 
COXHOD  16S 
COXHOD  169 
COXHOO  170 
COXHOD  171 
COXHOD  172 
COXHOD  173 
COXHOO  174 
COXHOD  175 
COXHOD  176 
COXHOD  177 
COXHOD  178 
COXHOD  179 
COXHOO  180 
COXHOD  181 
COXHOD  1S2 
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CALL  URimSt3.Pm8(l.l).NUWR0B.l) 

COXMOD  183 

DO  M  I  =  1.  NQIES 

COXHOD  184 

PROBU.I)  =  PRCB(I,I)-l. 

COXHOD  185 

PROBU.I)  *  0. 

COXHOD  186 

K(I)  *  0. 

COXHOD  187 

80 

CONTINUE 

COXHOD  188 

PROB(l.l)  -  L. 

COXHOD  189 

x<»  *  i. 

COXHOD  190 

C 

COXHOD  191 

C 

THE  RELATIVE  TRANSITION  RATE  BETWEEN  NODES  IS  CALCULATED 

COXHOD  192 

C 

SOLVING  THE  SET  OF  LIICAR  EQUATIONS  SPECIFIED  BY  PROB  AND  X. 

COXHOD  193 

C 

THE  SWOJTIIC  LECT1F  SOLVES  THIS  SET  OF  LOCAR  EQUATIONS. 

COXHOD  194 

C 

LEQT1F  IS  A  SUBROUTINE  FROM  THE  IHSL  PACKAGE  SUPPORTED  AT 

COXHOD  195 

C 

URI0HT-PAT1ERSCN  AFB,  OHIO. 

COXHOD  196 

c 

COXHOD  197 

CALL  LEQTlFIPROB.^NODES.nAXMSE.X.O.UAREA.IER) 

COXHOD  198 

CALL  RB»MS(3,PRCBU,1),NUNPR0B,1) 

COXHOD  199 

DO  100  I  *  1.  NODES 

COXHOD  200 

GAMA  *  0. 

COXHOD  201 

DO  90  J  a  It  NODES 

COXHOD  202 

GAMA  «  GAMA  +  (PROB(I.J)  »  X(J>) 

COXHOD  203 

90 

CONTINUE 

COXHOD  204 

IF  (  ABS(GAMA  -  X(I)  )  ,LE.  .1  )  GO  TO  100 

COXHOD  205 

URITE(6,95)  I.  (X(L),L»l, NODES) 

COXHOD  206 

95 

FGRHATdOXi'THE  RELATIVE  SERVICE  RATE  BETWEEN  NODES  CALCULATED  ", 

COXHOD  207 

1‘FRON  THE  COEFFICIENT  ARRAY  PROS  IS  IN  ERROR  FOR  NODE  ", 

COXHOD  208 

2I3//5X, "THE  RELATIVE  SERVICE  RATES  ARE  ",/iOX, 

COXHOD  209 

310(/10(F10.5,2X))) 

COXHOD  210 

GO  TO  115 

COXHOD  211 

100 

CONTINUE 

COXHOD  212 

DO  105  I>  1,  NODES 

COXHOD  213 

E(I,K)  *  X(I) 

COXHOD  214 

105 

CONTINUE 

COXHOO  215 

110 

CONTINUE 

COXHOO  216 

115 

CONTINUE 

COXHOD  217 

RETURN 

COXHOD  218 

END 

COXHOD  219 

SUBROUTINE  COWUTEIJCLASS,  PROB,  NDISKS,  DISKPRB,  TCPU,  DIO, 

COXHOD  220 

1EXPVAL,  SERVICE, IUSERS, NAXNODE, NAXTYPE, NOOES,  NAVGJOB) 

COXHOD  221 

DIMENSION  PRCBOIAXNGOEi NAXTYPE), DISKPRB(HAXNH)EiHAXTYPE), 

COXHOD  222 

lDIO(HAXTYP€),EXPVAL(NAXNODEil1At(TYPE). 

COXHOD  223 

2SSWICE(NAWNDE,  NAXTYPE), TCRJ(HAXTYPE) 

COXHOO  224 

INTEGER  XLASS, NDISKS,  IUSS1S,  NAXNODE,  NAXTYPE,  ITYPES 

COXHOO  225 

C 

COXHOD  226 

C 

THIS  SUBROUTINE  CONFUTES  THE  PROBABILITY  TRANSITION  MATRIX 

COXHOO  227 

C 

FOR  THE  MODEL  BASED  ON  THE  INPUTS  OF  CHEN'S  MODEL  (NDISKS,  DISKPRB, 

COXHOO  228 

C 

TCPU,  AND  DIO)  AM)  HIS  FORMULAE. 

COXHOD  229 

C 

COXHOO  230 

C 

INITIALIZE  PROBABILITY  TRANSITION  MATRIX 

COXHOD  231 

DO  10  I  »  1,  NODES 

COXHOO  232 

DO  10  J  *  1,  NODES 

COXHOD  233 

IF  (J  .EQ.  1)  00  TO  5 

COXHOD  234 

PRDB(J,I)  >  0.0 

COXHOO  235 

i 


GO  TO  10 

5  PROBUiI)  =  1.0 

10  CONTINUE 

PftOB(NDISK$+3»  1 )  =  $0WIC£(1.JCLASS)  /  TCPU(JCLASS) 

BO  20  I  *  t,  NQISKS 

PROBUM, 1)  =  DIO(  JCLASS)  *  PflOBf NDISKSt-3, 1  >  *  DISKPRB(IrJOASS) 
20  CONTINUE 

IP  ( (EXPVAL(NDISKS*3,  l >  .<£.  IUSERS  -  fMVGJOB)  .AND. 

1( IUSERS  .GE.  NWGJCB)) 

2FfiOB(NDISKS*-2> 101SKS+3)  =  (IUSERS  -  NAWGJQBI  /  EXPVAL(«ISXS+3.1> 
IF  (  IUSERS  .LE.  NAVGJOB  )  PR0B!NDISKS+2,NDISKS+3)  *  0.0 


IF  <  EXPVAL(NDISXS+3,1)  .LE.  (IUSERS  -  NAVGJOB)  ) 
1  PRCB(NDlSKS+2, NDISKS+3)  -  1.0 


PR0B(1,»ISKS*3)  =  1.0  -  PRC8(NDlSKS+2,NDISKS+3) 

CH0ICE1  *  0.0 
CH0ICE2  =  0.0 

ACTM  =  EXPVALd,  1)  +  EXPVAL(W)ISKS+2,1) 

DO  40  I  =  1.  HDI9CS 

ACTIVE  =  ACTIVE  ♦  EffWL(M.l) 

CHOICEl  =  CHOICEl  <-  PfiOB(I+lil) 

40  CONTINUE 

CHOICEl  =  CHOICEl  *  PR0B(NDISKS+3,1) 

CHOICEl  =  1.0  -  CHOICEl 

IF  ((ACTIVE  .LT.  NAVGJOB)  .OR.  (IUSERS  .LE.  NAVGJOB))  00  TO  50 
CH0ICE2  «  PR0B(NDISXSf3,l)  t  (ACTIVE  -  NAVGJOB)  / 

1  (IUSERS  -  NAVGJOB) 

FR0B(NDISKS+2, 1)  =  AHIN1 (CHOICEl ■ CHQICE2 ) 

GO  TO  40 

50  PR0B<N0ISKS+2,t)  =  0.0 

60  FftOBdd)  *  1.0  -  PR0B(N0ISKS+2,t)  -  d.O  -  CHOICEl) 

RETURN 

EW 

«EOR 

*E0F 

} 


COXHOD  234 
COXHOO  237 
COXNOD  238 
COXNOD  239 
COXHOD  240 
COXHOO  241 
COXHOD  242 
COXHOO  243 
COXHOD  244 
COXHOD  245 
COXHOD  244 
COXHOD  247 
COXHOD  248 
COXHOD  249 
COXHOD  250 
COXHOO  251 
COXHOD  252 
COXNOD  253 
COXHOO  254 
COXHOD  255 
COXHOD  256 
COXHOD  257 
COXHOD  258 
COXHOD  259 
COXHOO  240 
COXHOO  261 
COXHOD  262 
COXHOD  263 
COXHOD  264 
COXHOD  265 
COXHOD  266 
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COXHOD  269 
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